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Numerous correlated electron systems exhibit a strongly scale-dependent behavior. Upon lower- 
ing the energy scale, collective phenomena, bound states, and new effective degrees of freedom 
emerge. Typical examples include (i) competing magnetic, charge, and pairing instabilities in 
two-dimensional electron systems, (ii) the interplay of electronic excitations and order parameter 
fluctuations near thermal and quantum phase transitions in metals, (iii) correlation effects such as 
Luttinger liquid behavior and the Kondo effect showing up in linear and non-equilibrium transport 
through quantum wires and quantum dots. The functional renormalization group is a flexible and 
unbiased tool for dealing with such scale-dependent behavior. Its starting point is an exact func- 
tional flow equation, which yields the gradual evolution from a microscopic model action to the 
final effective action as a function of a continuously decreasing energy scale. Expanding in powers 
of the fields one obtains an exact hierarchy of flow equations for vertex functions. Truncations of 
this hierarchy have led to powerful new approximation schemes. This review is a comprehensive 
introduction to the functional renormalization group method for interacting Fermi systems. We 
present a self-contained derivation of the exact flow equations and describe frequently used trun- 
cation schemes. Reviewing selected applications we then show how approximations based on the 
functional renormalization group can be fruitfully used to improve our understanding of correlated 
fermion systems. 
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FIG. 1 (Color online) Important energy scales in high- 
temperature superconductors of the cuprate family. Magnetic 
interactions and superconductivity are generated from the ki- 
netic energy (hopping) and the Coulomb repulsion. 



scales are plagued by infrared divergences, and are there- 
fore often inapplicable even at weak coupling, especially 
in low dimensions. 

It is thus natural to treat degrees of freedom with dif- 
ferent energy scales successively, descending step by step 
from higher to lower scales. This is the main idea behind 
the renormalization group (RG). 



I. INTRODUCTION 
A. Motivation 

The Coulomb interaction between electrons in solids 
leads to a virtually unlimited variety of phenomena, 
such as magnetic correlations and magnetic order, high- 
temperature superconductivity, metal-insulator transi- 
tions, phase separation and stripes, and the formation of 
exotic quantum liquid phases. The latter include Lut- 
tingcr liquids, quantum critical points, and fractional 
quantum Hall states. 

Interacting electron systems usually exhibit very dis- 
tinct behavior on different energy scales. Composite ob- 
jects and collective phenomena emerge at scales far be- 
low the bare energy scales of the microscopic Hamilto- 
nian. For example, in cuprate high-temperature super- 
conductors one bridges three orders of magnitude from 
the highest scale, the bare Coulomb interaction, via the 
intermediate scale of short-range magnetic correlations, 
down to the lowest scale of d-wave superconductivity and 
other ordering phenomena (sec Fig.[T]). This diversity of 
scales is a major obstacle to a straightforward numerical 
solution of microscopic models, since the most interest- 
ing phenomena emerge only at low temperatures and in 
systems with a large size. It is also hard to deal with by 
conventional many-body methods, if one tries to treat all 
scales at once and within the same approximation, for 
example by summing a subclass of Feynman diagrams. 
Pcrturbativc approaches which do not separate different 



B. RG for Interacting Fermi systems 

Renormalization group methods have a long tradition 
in the theory of interacting Fermi systems. Already in 
the 1970s, various versions of the RG have been used to 
deal with infrar ed singularities arising in one-dimensional 
Fermi systems (jSolvoml Il979l ). Naturally, the RG was 
also applied to (mostly bosonic) effective field theories 
describing critical phenomena at continuous classical or 
quantum phase transitions in interacting Fermi systems 
(iFradkinL flQQlt ISachdevl . ll999D . 

Renormalization group approaches dealing with inter- 
acting fermions in arbitrary dimensions d have been de- 
veloped much later. Due to the extended (not point-like) 
geometry of the Fermi surface singularity in dimensions 
d > 1, the renormalization group flow cannot be reduced 
to a finite number of running couplings. However, the 
main reason for the delayed development of a compre- 
hensive RG approach for interacting Fermi systems in 
higher dimensions was probably not this difficulty, but 
rather a lack of motivation. The few infrared singulari- 
ties appearing in three-dimensional Fermi systems could 
usually be han dled by simple r e summations of pertur- 
bation theory (lAbrikosov et al \ 119631: iNozieresl 119641) . 
Triggered by the issue of non-Fermi liquid behavior in 
two-dimensional systems, and the related discussion on 
the validity of perturbation theory, systematic RG ap- 
proaches to interacting Fermi systems in arbitrary di- 
mensions have been developed by various groups in the 
early 1990s. 

Aiming at a mathematical control of interacting 
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Fermi systems, iFeldman and Trubowitd (Il990l 1991 ). 
and independently iBenfatto and Gallavottil (|l990al lll. 
have formulated a riR orous fermionic version of Wilson's 
momentum-shell RG ( Wilson and Kogutl . Il974l) . Impor- 
tant rigorous results have indeed been obtained in one- 
dimensional ([Benfatto eii"al j_ 1994 ) and two-dime n sional 
([Bcnfatto et al!. '2006': 'Piscrtori and Rivasseaul . l2000l: 
iFold man et a/., 2003, 2004, 1992, 1996.) systems. The 
essential message from these results is that no hitherto 
unknown instabilities or non-perturbative effects occur in 
Fermi systems with sufficiently weak short-range interac- 
tions, at least in the absence of special features such as 
Van Hove singularities at the Fermi level. 

The Wilsonian RG for interacting Fermi systems was 
populari z ed among (non-i nathematical) ph ysicists by 
IShankaj (|l99ll Il994[ ) and iPolchinskil (|l993[) . who pre- 
sented some of the main ideas in a pedagogical style. In 
particular, they provided an intuitive RG perspective of 
Fermi liquid theory. Subtleties associated with the sin- 
gularities of the interaction vertex for forward scatter- 
ing were clariiied a b it later ( Chitov and Senechallll995t 
iMetzner et aZ.I . [T998[ ) . A Hamiltonian-based RG interpre- 
tation , of Fermi liquid theory was presented by iHewsonI 
(I1994D . who discussed not only translation invariant sys- 
tems but also models for magnetic impurities in metals. 

As an alternative to the Wilsonian RG one may also 
use flow equations for Hamiltonians based on infinites- 
imal unitary transformations, which make the Hamilto- 
nian successively more diagonal (IWegneil . [l994h . This ap- 
proach has been used succ essfully for qua ntum impurity 
models and other systems ( Kehreinl |2006| ). A weak cou- 
pling truncation of the flow equations has been applied 
to iden tify instabili t ies of the two-dimensional Hubbard 
model (|Grote et aZ.! . |2002[) . 

There is much current interest in RG methods for 
correlated fermions in non-e q uilibrium. The pe rturba- 



tive RG (jMitra et all 
numerical RG 



son s 



20061: iRosch et a/.l. l200lh. Wil 



Anders and Schillerl . |2005|). as 



well as Wegner's flow equation approach ( Kehreinl . 
I2006D were extended to non-eq uilibrium, a nd re a l-time 
RG methods were develo ped (jSchoelleii l2000l 120091: 
ISchoeller and Konid . [2OO0I) . 



C. Functional renormalization group 

The Wilsonian RG is not only useful for a deeper 
and partially even rigorous understanding of interact- 
ing fermion systems. A speciflc version of Wilson's RG 
known as exact or functional RG turned out to provide 
a valuable framework for computational purposes. Ap- 
proximations derived from exact functional flow equa- 
tions have played an increasingly important role in the 
last decade. These developments are the central topic of 
this review. 

Exact flow equations describe the evolution of a gen- 
erating functional for all many-particle Green or vertex 
functions as a function of a flow parameter A, usually 



an infrared cutoff. They can be derived relatively easily 
from a functional integral representation of the generat- 
ing functional. Exact flow equations have been known 
si nce the early years of the RG , starting with the work 
of IWegner and Houghto^ (|l973l ). IPolchinskil (|l984[) em- 
ployed an exact flow equation to formulate a relatively 
simple proof of renormalizability of the $^-thcory in four 
dimensions. Renormalizability proofs can be further sim- 
plified by using a Wick-ordere d variant of Polchinski's 
equation ( Wieczerkowskil [l988() . 



For computational purposes the exact fiow equation 
for the effective action, first derived in the context of 
bosonic field theories bv lWetterichI (|1993[ ) turned out to 
be most convenient. The effective action r^[0] is the 
generating functional for one-particle irreducible vertex 
functions. The latter are obtained by taking deriva- 
tives with respect to the source field 0. The flow pa- 
rameter A describes a rcgularization of the underlying 
bare action, which regularizes infrared divergencies in 
perturbation theory. The rcgularization is removed at 
the end of the flow, say for A — >■ 0. The initial regu- 
lator (for A = Aq) can be chosen such that r'^"[(/)] is 
given by the bare action. The flow of then pro- 

vides a smooth interpolation between the bare action 
of the system and the final effective action F [</>], from 
which any desired information can be extracted. This 
fiow is d etermined by an exact functional diffe r ential 
equation (lEUwanger and Wetterichl . Il994l: iMorrisl Il994l: 
IWetterichl . 1993). Expanding in the fields one obtains 
a hierarchy of flow equations for the one-particle irre- 
ducible vertex functions. The advantage of that hier- 
archy compared to others, obtained, for example, from 
Polchinski's equation, is that self-energy feedback is in- 
cluded automatically and no one-particle reducible terms 
appear. 

The expression functional RG stems from the feature 
that the exact flow equations describe the flow of a func- 
tional or (equivalently) of a hierarchy of functions. An 
important difference compared to Wilson's original for- 
mulation is that a complete set of source flelds is kept in 
the flowing generating functionals, not only those corre- 
sponding to scales below A. Hence, the full information 
on the properties of the system remains accessible, not 
only the low energy or long wavelength behavior. 

Exact flow equations can be solved exactly only in 
special cases, where the underlying model can also be 
solved exactly (and more easily) by other means. "'^ How- 
ever, the functional RG is a valuable source for devising 
powerful new approximation schemes, which can be ob- 
tained by truncating the hierarchy and/or by a simplified 
parametrization of the Green or vertex functions. These 
approximations have several distinctive advantages: i) 



^ An instructive example is provided by the exact solution of the 
Tomonaga-Luttinse r model via functional RG flow equations 
jSchiitz aLl . l2005t) . 
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they have a renormahzation group structure built in, that 
is, scales are handled successively and infrared singular- 
ities are thus treated properly; ii) they can be applied 
directly to microscopic models, not only to effective field 
theories which capture only some asymptotic behavior; 
iii) they are physically transparent, for example one can 
see directly how and why new correlations form upon 
lowering the scale; iv) one can use different approxima- 
tions at different scales. Small steps from a scale A to a 
slightly smaller scale A' are much easier to control than 
an integration over all degrees of freedom in one shot, 
and one can take advantage of the flexibility provided by 
the choice of a suitable flow parameter. 

Approximations derived from exact flow equations 
have been applied in m any areas of quantum field theory 
and statistical physics ( Berges et all [2OO2I ) . In the con- 



text of interacting Fermi systems, functional RG methods 
were first used for an unbiased s tability analysis of the 
two-dimensional Hubbard model (IHalboth and Metzner , 
2000al : lHonerkamD et aZl . l200ll : [Zanchi and Schulzj . ll998 . 
2OOOD . Since then, approximations derived within the 



functional RG framework have been applied to numer- 
ous interacting fermion systems. 



D. Scope of the review 

This review provides a thorough introduction to the 
functional RG in the context of interacting Fermi sys- 
tems. It should serve as a manual and reference for many- 
body theorists who wish to apply approximations based 
on the functional RG to their own problem of interest. 
We will first describe the functional RG framework and 
derive in particular the exact flow equations, which are 
the starting point for approximations. We will discuss 
general aspects related to the flow equations such as the 
choice of cutoffs, power counting and truncations. Links 
to the use of flow equations in the mathematical liter- 
ature will be pointed out along the way. We will then 
review some of the most interesting applications of trun- 
cated functional RG equations (see table of content) . Our 
aim is not to deliver an exhaustive overview of all appli- 
cations, but rather to show via selected applications how 
the functional RG method can be fruitfully used. 

The functional RG wa s recently extended to Fermi sys- 
tems out of equilibrium ([Gezzi et aZI. 2007 1 Jakobs, 2003"; 



tems out 01 equnibrium (itjezzi et at.l.lzuuaiJa K:obs. zuuj 
I Jakobs a/1 l2007al . l2010allbl : iKarrasch et aU . JOlOaaa) 
In the derivation of the flow equations in Sec. II we 
restrict ourselves to the equilibrium formalism. Func- 
tional RG flow equations for non-equilibrium Keldysh 
Gree n and vertex functions can b e derived in close anal- 
ogy (iGezzi et all l2007t iJakobsl . l2003l I2OIOI: iKarraschl 
I2OIOD . The necessary extensions are briefly mentioned 
when discussing the application of this method to finite 
bias steady state transport through correlated quantum 
wires and quantum dots in Sec. VI. 

A number of reviews with a focus on the functional 
RG are already available. Mathematically rigorous de- 



velopments until the e nd of the last mill enium were sum- 
marized in a book bv ISalmhofed ( 1999t ). a large portion 
of which is dedicated to interacting Fermi systems. Ex- 
amples of approximations derived from the exact flow 
equation for the effective action with many applications 
in quantum field theory and sta tistical physic s were pre- 
sented in the review article by iBerges et al. I (I2002D . A 
detailed introduction to the functional RG in a text- 
book style supplemented by selected applications (includ- 
ing Fermi sys t ems) can be found in the recent book by 
iKopietz et all (|2010t) . 



II. FUNCTIONAL FLOW EQUATIONS 

In this section we present the general functional 
RG framework. The reader should be familiar with 
the functional integral formalism for quantum many- 
body systems, as d e scribe d in classic textbooks such as 
iNegele and Orlaiidl ( 1987t ). After introducing the gen- 
erating functionals for Green and vertex functions in 
Sec. II. A, we derive the exact functional flow equations in 
Sec. II. B. The flow equation (|35|) for the effective action 
is the central equation of this review. Expanding in 
the fields we derive the hierarchy of fiow equations for 
vertex functions in Sec. II. C, which is the starting point 
for approximations. Possible choices of flow parameters 
arc reviewed in Sec. II. D. The general structure of the RG 
hierarchy and power counting are discussed in Sec. II. E, 
with various references to the closely related mathemat- 
ical literature. Sec. II. F is dedicated to flow equations 
for observables and correlation functions. Coupled flow 
equations for fermions and bosons, which are useful for 
studies of spontaneous symmetry breaking and quantum 
criticality, are derived in Sec. II. G. 



A. Generating functionals 

We consider a system of interacting fermions which can 
be described by Grassmann fields V'j V'l ^^id an action of 
the form 



= -(Vi,GoV) + F[^,V;] 



(1) 



where is an arbitrary many-body interaction, 

and Go is the propagator of the non-interacting sys- 
tem. The bracket (.,.) is a shorthand notation 
for the sum ^^^i.^) ^^0^"^)^^)^ where {G'^^il^){x) = 
X^a;' Gq {x^x')'il}{x'). The Grassmann field index x col- 
lects the quantum numbers of a suitable single-particle 
basis set and imaginary time or frequency. In case of 
continuous variables, the sum over x includes the ap- 
propriate integrals. Prefactors such as temperature or 
volume factors depend on the representation (e.g. real or 
momentum space) and are therefore not written in this 
general part. A two-particle interaction has the general 
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form 



VbPyi^] 7 X! Vix'l, Xi, X2) ij{x[)ij{x'2)'llj{x2)ij{xi) . 



El ,^2 
^1 '^2 



(2) 

In particular, for spin-i fcrmions with a single-particle 
basis labeled by momentum k and spin orientation a, 
one has x = (fco,k, cr), where kg is the fermionic Mat- 
subara frequency. If the bare part of the action is trans- 
lation and spin- rotation invariant, the bare propagator 
has the diagonal and spin-independent form Go(x,.t') = 
4ofe^,<5kk'<5o-o-'Go(fco,k) with 



Go(fco,k) = 



1 



(3) 



where ^k = fk ~ A* is the single-particle energy relative to 
the chemical potential. 

Connected Green fu nctions can be obtained f rom the 
generating functional ( Negele and Orlaiidl . Il987() 



g[T],f]] = -Iny VtpViP e-^^''''''''^ e(*'''^)+('^'''' , (4) 

where / VipVijj . . . — J Y[x dip (x) dip (x) .... Completing 
squares yields the identity 

VipViPe^'f'-'°o'^)e(fi,i')H4':n) ^ Zq e^^*'''^"") , (5) 



where Zq = J VipVip e^^'^o '^'^ is the partition function 
of the non- interacting system. Hence C/ [77, 77] = — In Zg -I- 
(77, Gq??) in the non-interacting case = 0. For 

vanishing source fields, ^J[0,0] = — InZ, where 



/ 



(6) 



is the partition function of the interacting system. The 
connected m-particle Green functions are given by 



^ (3^1 , ... J Xm ]Xi, . . . , Xj^) — 
-{lp{xi) . . .1p{Xm)'ip{x'^) . ..1p{x[)), 



dri{xi) . . . dr]{xm)dri{x'^) . . . dr]{x[) 



, (7) 



r},r}—0 



where (. . .)c is the connected average of the product of 
Grassmann variables between the brackets. The one- 
particle Green function G^^-* is the propagator of the 
interacting system, which we will usually denote with- 
out the superscript by G. Expanding G[ri, fj] in the fields 
yields a formal power scries with the connected Green 
functions as coefficients, 



GlVi v] = - In Z + (?/, Gr/) + 



^ y 

(2!)2 ^, , 



G'^'^'> {xi, X2; x[, X2) 'n{xi)T]{x2)r]{x2)T]{x[) 



(8) 



Renormalization group equations are most conve- 
niently formulated for the Legendre transform of Q[ri,f]\, 
the socalled effective action 



T[ip, ip] = {-q, ip) + {ip, 77) -I- g[T], 77] , 



(9) 



with ip — —dQ I drj and -0 = dQ / drj , which generates one- 
parti cle irreducible vertex functions (jNegele and Orlandl . 
Il987l) 



^ ('''1 1 • ■ • 1 x\ 



dil){x\) . . . dip{x',j^)dip{xm) ■ ■ ■ dip{xi) 



(10) 



In the non-interacting case one obtains r['(/','!/)] = 
— InZo — {iP,Gq^iIj). The Legendre correspondence be- 
tween the functionals G and T yields relations between 
the connected Green functions G'^™' and the vertex func- 
tions r'^'"'. In particular, 

p(2) ^ (^-1 = G-i _ s , (11) 
where E is the self-energy. The connected two-particle 



Green function is related to the two-particle vertex by 
G^'^\xi,X2;x[,X2) = V G{xi,y[)G{x2,y2) 



E 

yi,y2,yi,y'2 



xrW(2/'i,yi;yi,y2)G(yi,x'i)G(y2,4) 



(12) 



while the three-particle Green function G'^^'' = 
(^3p(6)(53 _^(-3p(4)g<p(4)g.3 involves r(4) and T^^). More 

generally, the connected r/i-particlc Green functions are 
obtained by adding all possible trees that can be formed 
with vertex functions of eq ual or lower order and G-lines 
(jNegele and Orlandl . fl987l) . 

The effective action obeys the reciprocity relations 



dip 



V 



(13) 



The second functional derivatives o f G and F with re- 
spect to the fields are also reciprocal ( Negele and Orlaiidl . 
I1987D . We define the matrices of second derivatives at fi- 
nite fields 
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and 



dfj{x)dr}{x') 

ay 

dri{x)d7]{x') 



dr]{x)dri{x' 
drj{x)drj{x') 



__( mxWx')) (v(x)v(x')) 

(^(x)^(x')> (^(x)^(x')) 



(14) 



2p g2p 



d'ili{x')di)[x) ajp{x')d'4>{x) 



ddr[ip,ip]{x',x) aar[^,v](x',x) 



(15) 



where the matrix elements in the second matrix of the 
last equation are just a more conventient notation for 
those in the first matrix. The reciprocity relation for the 
second derivatives reads 



with the functional Laplacian 



AG„ = (a„Goa,)=^^Go(x,x')^. (20) 



(16) 



Note that anomalous components are involved as long 
as the source fields are finite. Only at 77 = ^ = and 
■0 = -0 = 0, and in the absence of U{1) charge symmetry 
breaking one has the simple relation F^^^ = (G^^^) 

Another useful generating functional is the effective 



Another useiui generating 
interaction i Salmhoferl . Il999f ) 



V[x,x] = -ln|^ J VipV^jje^'^-^o''''^e-^'^'''+^''''+^^ 

(17) 

A simple substitution of variables yields the relation 

V[x,x]^Ghfl]+lTiZo-ifj,Gov) , (18) 

where x = GoV and x = GqV- Here Gg is the transposed 
bare propagator, that is, Gq{x,x') = Go{x',x). Hence, 
functional derivatives of V[x, x] with respect to x ^^nd x 
generate connected Green functions with bare propaga- 
tors amputated from external legs in the corresponding 
Feynman diagrams. The term In Zq — {fj, GqT]) cancels the 
non-interacting part of ^[f?,?/] such that V[x, x] = for 
= 0. The effective interaction V can also be ex- 
pressed via functional derivatives, instead of a functional 
integral: 



It is sometimes convenient (see Sec. II. G) to combine 
the fields tp and ^ in a Nambu-type field 



ip{x) 
ip{x) 



and similarly for the source fields rj and 77, 



H{x) 



f]{x) 
-fi{x) 



(21) 



(22) 



The minus sign in the definition of H makes sure that the 
source term [fj, ■0) + (V'l appearing in the definition of 
and also in the Legendre transform relating Q and F, 
can be written concisely as (_ff,^). In Nambu notation, 
the matrices of second derivatives of Q and F have the 
compact form 



and 



respectively. 



G(2)[i7] 



dH{x)dH{x') 



d^{x')d'^{x) 



(23) 



(24) 



1 

1 



X e 



in,i'+x}+{v^ip+x) 



ri,ri—0 



^ ^-Vldf,M„] g(»),Gor,)g(f;,x) + (r,,x) 



7],7] — 



-V[dr,M^] Jd^,God^)Jf„x) + in,x) 



„Ago p-Vlx,x] 



ri,ri—0 



(19) 



B. Exact fermionic flow equations 

In this section we derive exact flow equations de- 
scribing the evolution of the generating functionals de- 
fined above, as a function of a flow parameter A which 
parametrizes a modification of the bare propagator Go- 
Usually A is an infrared cutoff or another scale depen- 
dence. For example, in a translation invariant system 
one may impose a momentum cutoff, modifying Go to 



G^(fco,k) 



61^ (k) 

iko - 



(25) 
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where 6'^(k) is a function that vanishes for |^k| *C A 
and tends to one for ^ A. In this way the infrared 
singularity of the propagator at fco = and = (corre- 
sponding to the non-interacting Fermi surface in k-space) 
is cut off at the scale A. A simple choice for 6''^(k), which 
was often used in numerical solutions of truncated flow 
equations, is 



d\k) = e(iai - A) , 



(26) 



where 9 is the step function. With this choice momenta 
close to the Fermi surface are strictly excluded, as il- 
lustrated in Fig. [5] for a two-dimensional lattice fcrmion 
system. Alternatively, one may also use a smooth cutoff 



V obviously plays a dual role: It is the generating func- 
tional for (amputated) Green functions of a system with 
a cutoff A, and at the same time the inter action for the 
remaining low en ergy degrees of freedom ( Morrisl Il994l : 
ISalmhofeAlT999l) . 

The effective interaction satisfie s the following exact 
renormalization group equation ( Brvdges and Wrightl . 
[1988; Salmhofel 119991) 



dx 



- tr G, 



dxdxJ 



where Gg 



■j^Gq and tr denotes the trace trA 



(30) 



J2x A{x,x). Its derivation is simple: 




FIG. 2 (Color online) Momentum space region around the 
Fermi surface excluded by a sharp momentum cutoff for 
fermions with a tight-binding dispersion on a two-dimensional 
square lattice (lattice constant = one). 

function. In the absence of translation invariance it is 
more convenient to use a frequency cutoff instead of a 
momentum cutoff. The cutoff excludes "soft modes" be- 
low the scale A from the functional integral. Instead of a 
cutoff one may also choose other flow parameters such as 
temperature. The various possibilities will be discussed 
more extensively in Sec. III.DI For the derivation of the 
flow equations it does not matter how Gq depends on A. 

The bare action constructed with Gg (instead of Gg) 
will be denoted by S^['ip,ilj], and the generating func- 
tionals introduced in Sec. III.Al bv ^^[77,77], V^[x, x], and 
r^['ip,'ip], respectively. The original functionals Q, V and 
F are recovered in the limit A — 0. 

In the presence of a cutoff, Eq. (|19p becomes 



'0 e 



(27) 



At the highest energy scale Ag one has Gg" = 0, and 
thus V"^" — V. Hence, interpolates smoothly between 
the bare interaction V and the generating functional V. 
Introducing the soft mode propagator 



GA y^A 
~ '-'O ^ <-fo 



which has support on scales below A, wc can write 



e e 



(28) 



(29) 



dA 



yA U _yA 

dA 



-e' -rrle e 
dA 



— e Aqa e 



= right-hand side of Eq. i^U^ . 

In the second step we have used Eq. (P7)) . With the initial 
condition 



(31) 



the RG equation determines the flow of uniquely for 
all A < Ag . The initial value Ag must be chosen such that 
Gq° vanishes. For a sharp momentum cutoff, Ag can be 
chosen as the maximal value of |^k|; for a frequency cutoff 

Ag = 00. 

An expansion of the functional ^"^[XjX] in the renor- 
malization group equation ([5(11) in powe rs of x and y leads 
to the fcrmionic analog of Polchinski 's (iPolchinskiLlTgsl 
flow equations for amputated connected m-particle Green 
functions T/(2™)A. 

From the flow equation for V"^, Eq. ([3D]), and the re- 
lation ([T8|) applied to and Q^, one obtains an exact 
flow equation for Q^: 



fdg^ 



,0, 



df] 



-tr 



(32) 



where Qq = (Gg )^^ and the dot denotes a A-derivative. 
This flow equation can also be derived more directly, by 
applying a A-dcrivative to the functional integral repre- 
sentation of Q^. 

The flow equations for G'^™''^ and T/(2m)A generate, 
among others, also one-particle reducible terms, which 
require some special care. In this respect the flow equa- 
tions for one-particle irreducible vertex functions F^^™^^, 
obtained from the scale-dependent effective action. 



(33) 



are easier to handle. Note that 77 and fj are A- 
dependent functions of -ip and Vi a-s they arc deter- 
mined by the A-dcpendent equations tp ~ —dQ^/dfj and 
■i/j = dQ^/drj. Since the A-dcpendcncc does not change 
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the structure of the action as a function of the fields, all 
standard relations between the connected Green func- 
tions G'^™-' and the vertex functions F^^™^ carry over to 
the ones for G't^™)^ and F^^™)^. 

The A-derivative of F"^ can be written as ^F^ [?/', V"] = 
imt. ^) + (V^, sv^^) + -k ^""N^' '7^]' ^here the deriva- 
tive in front of Q'^ acts also on the A-dcpendence of 
■q^ and 77^. Due to the relations dQ^/drj = 'tp and 
dQ^/dfj = —ijj, most terms cancel and one obtains 



,'7 



f^A fixed 



(34) 



Inserting the flow equation ([5^ for and using the 
reciprocity relations (|13p and (fT6|) . one obtains the exact 
functional flow equation for the effective action 



^F^[V',V"]=-(V5,Qo^)-^tr 

(35) 

Here r'^^^ [■(/;, V'] is the matrix of second functional deriva- 
tives defined in Eq. (fTSl) . and 



-Q^}* 



diag(Q^,-Q;}*) 



(36) 



where Q^*(a;,a;') = Q^{x',x). 

Alternative definitions of the effective action F^, dif- 
fering by interaction-independent terms, have also been 
used. One variant is to normalize the functional integral 
defining ai V = dividing by Zq. This yields an 
additional contribution In Zq to and to its Legendre 
transform F"^. In the flow equation for F"^ this leads to 
an additional term tr((5^GQ ), which is field independent 
and therefore does not coup l e to t he other contributions 
Salmhofer and Honerkamp . 2001 ). Another variant is 



Berges et a/J . l2002t lEUwanger and Wetterichl . [l993 ) 



F^[^,^]=FA[V',Vi] + (V",i?M 



(37) 



where = Qq — Qq. The additional quadratic term 
cancels the first (trivial) term in the flow equation ([55)) 
for F^, and one obtains the equivalent flow equation 



^r^[^,^].-itr 



.{2)A 



(38) 

where R"^ = diag (i?^, --i?^*) . The functional F^ and 
its analog for bosonic fiel ds is known as effec tive average 



action in the literature ( Berges et al\ . |2002() . Both F^ 



and F^ tend to the same effective action F in the limit 
A 0, where i?^ vanishes. At the initial scale Aq, one 

has F^" [V', 



5 [■(/', ■0] J while 



= 5[V;,^]-(^,i?^» . (39) 



the final effective action F, while F^ interpolates between 
the regularized ha.Te action and F. On the other hand, 
the functional F^ has the advantage that its second func- 
tional derivative directly yields the inverse propagator 
{G^y^ without the need to add R^. 

In Appendix [A] we present yet another version of ex- 
act flow equations, based on a Wick ordered effective 
interaction. That version also contains one-particle re- 
ducible contributions, but it has the distinct advantage 
that the vertices are connected by propagators with an 
energy scale at o r below A. This facilitates a systematic 
power counting ( Salmhofeil Il999t ). and also a numerical 
evaluation of flow equations, since the integration regions 
shrink upon lowering A. 

It is instructive to compare the functional RG flow 
equations with the traditional W ilsonian momentum 
shell RG (IWilson and KogutL Il974^. wh ich was applied 
to Fer mi systems by Shankad ( 199ll . 1994| ) and lPolchinskil 
(I1993D . ,n the commonly used version of Wilson's RG, 
the flow of the effective action is computed only for soft 
fields, that is, for fields with energy or momentum vari- 
ables below the scale A, while in the functional RG the 
effective action with unrestricted source fields is com- 
puted. This allows for a direct calculation of correlations 
functions with arbitrary external variables such as mo- 
menta or Matsubara frequencies. Furthermore, in the 
traditional implementations of Wilson's RG the integra- 
tion of degrees of freedom is combined with a rescal- 
ing of momenta and fields, which is chosen such that 
the momentum space and certain terms in the quadratic 
part of the action remain invariant during the fiow. 
This facilitates the classification of interactions as rel- 
evant, marginal or irrelevant, and helps to identify fixed 
points of the flow. The functional RG flow equations 
derived above do not involve any rescaling. Rcscaling 
momentum space in a shell around the Fermi surface re- 
quires a non- linear transformation in dimensions d > 1, 
which s poils the simple linea. r forrn of momentum conser- 
vation (IKopietz and Buschd . I2OOII: iMetzner efall . Il998l: 
IShankaij . Il994 ). and is therefore of questionable value. 
Power-counting can be done also without rescaling, as 
shown in Sec. II. E. Rescaling of the flelds can be im- 
plemented easily by a simple substi t ution of variables 
(|Kopietz and Buschd . I2OOII: IShankail . Il994l ). However, 
in many applications of the functional RG, quantitative 
results including power-laws with anomalous scaling di- 
mensions are obtained simply by direct calculation of the 
(unsealed) physical quantities. 



C. Expansion in the fields 

1. Hierarchy of flow equations 



Hence, F;^ has the attractive feature that it interpolates 
smoothly between the (unregularized) bare action S and 



The functional flow equation for the effective action 
can be expanded in powers of the fields. To this end we 
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expand the effective action as 

oo 



wfiere ^(^'"^^[V', ip] is liomogeneous of degree 2m in tlie 
fields, 



(40) 



rn=0 



(-ir 



m 



n2 



r^^"'^^{x[,. ..,x'^;Xi,..., Xm) H^'l) ■ . . i^{x'^)i^{Xrn) ■ ■ ■ '^{xi 



(41) 



r 



for m > 1. Tlie field- independent constant „4'-*'-'^ yields 
the grand canonical potential: 



(42) 



Here we have restored the explicit temperature factor, 
since it is independent of the representation of the fields. 
To expand the inverse of r*^^)^ on the right hand side of 
the flow equation, we isolate the field- independent part 
of r(2)A as 



r(2)A[^,V5] = (G^)-i-s^[V',^], 



(43) 



where 



= diag(G'\ -G^*) (44) 



is the full propagator, and (cf. Eq. (fT5|)) 



(45) 

Note that ■0] contains all contributions to 

r(2)A[^^ ^] which are at least quadratic in the fields. We 



1 - G^S 



AvA^ 



can now expand (r*^^)^) 
geometric series. Inserted in ([35]) . this yields 

1 

2'^ 



qA 



as a 



-tr(o^}G^j -(v;,goV') + 



where 



diag(S'^,-S'^') 



-G^Q^G^ 



(46) 



(47) 



Using the Dyson equation [G^)"^ — Qq — E^, the so- 
called single-scale propagator can also be written as 
A-derivative of the propagator at fixed self-energy, 



qA _ x-<A| 

~ dA '^'^ fixed 



(48) 



The expansion of the flow equation in powers of V'l V' 
is now straightforward and leads to a hierarchy of flow 



equations for E"^, the two-particle vertex F*^^)^, and the 
higher-order vertices r'^^^, F'^®)'^, etc. The first three 
equations in this hierarchy are shown diagrammatically 
in Fig. [3] Note that only one-particle irreducible one- 
loop diagrams contribute, and internal lines are dressed 
by self-energy corrections. The hierarchy does not close 
at any finite order, since the flow of each vertex F*^^™)"^ 
receives a contribution from a tadpole diagram involving 
p(2m-i-2)A^ and TO-particle vertices with arbitrary m are 
generated by the flow, irrespective of their presence in 
the bare action. 



d A 



d p(4)A , 

dA 



d p(6)A 

dA 



+ 

(4)/\" — ^ p(4)A 
(4)A 



r^""^ qA r 



p(6)A 



+ 

(4)A 



rW)A qA 



,(4)A oA 

pA r(8)A 



FIG. 3 (Color online) Diagrammatic representation of the 
flow equations for the self-energy E"^, the two-particle vertex 
r^*''^, and the three-particle vertex r'®-"^ in the one-particle 
irreducible version of the functional RG. Lines with a dash 
correspond to the single scale propagator S^, the other lines 
to the full propagator G^. 



Let us derive explicitly the first two flow equations 
from the hierarchy. Comparing coefficients of quadratic 
contributions (proportional to ipip) to the exact flow 
equation yields 



2)A 



-(V^,QoV')-tr(5^aayl('*)^) • (49) 



Inserting Eq. (|4T|), and using F^^^"^ = Qo - E^, one ob- 
tains the flow equation for the self-energy, 

^E^(x', x) = Yl '^^(2^' y') 2/'; y) ■ m 

y,y' 
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Comparing coefScicnts of quartic contributions (propor- 
tional to {tjjilj)'^) yields 



dA 2 V 

- itr (^S^ddA'-''^^G^'ddA^^^^ + S^'ddA^^^^G^BdA^^^^) - tr (^S^ddA'-'^^^ 
Inserting Eq. (|41l) , one obtains the flow equation for the two-particle vertex, 



(51) 



Vl.Vl V2,V2 



X {r(4)A(a;;,x'2;yi,y2)r(4)^(y'i,y^;xi,x2) 

„ r^^)^(x[,y2;xi,yi)T'^^^^{y[,X2;y2,X2) + {yi ^ y2,y'i ^ y'2) 

+ [rWA(4,y^;xi,yi)r(4)A(y;,x;;y2,X2) + (^1 ^ 2/2 , ^ J/^ )] } 

- ^5"^(2/,j/')r(*'^'^(xi,a;2,2/';2;i,X2,y) . 



(52) 



r 



Note that there are several distinct contributions involv- 
ing two two-particle vertices, corresponding to the fa- 
miliar particle-particle, direct particle-hole, and crossed 
particle-hole channel, respectively, as shown diagram- 
matically in Fig. 21 Similarly, one can obtain the flow 



The flow equation p5p and the ensuing equations 
for the vertex functions can be easily generalized to 
cases with J7(l)-symmetry breaking by allowing for off- 
diagonal elements in the matrices Qg , and S^. 




pp 



ph 



ph- 



FIG. 4 Contributions to the flow of the two-particle vertex 
with particle-particle and particle-hole channels written ex- 
plicitly, without the contribution from F^^'"^. 

equation for F*^^) and all higher vertices. 

Since F [■(/;,?/;] at = = is essentially (up to a 
factor T) the grand canonical potential VL, the flow equa- 
tion evaluated at vanishing fields, yields also a flow 
equation for the grand canonical potential: 



d 

dA' 



(53) 



2. Truncations 

The exact hierarchy of flow equations for the vertex 
functions can be solved only for systems which can also 
be solved more directly, that is, without using flow equa- 
tions. Usually truncations are unavoidable. A natural 
truncation is to neglect the flow of all vertices F*^^™^"^ 
beyond a certain order niQ. We call this the level-mQ 
truncation. The structure of the resulting equations and 
general properties of their solution will be discussed in 
Sec. II. E. Note that the level-mo truncation contains all 
perturbative contributions to order mg in the bare two- 
particle interaction. 

In practice, in applications to physically interesting 
systems, vertices F^^™^^ with m > 3 have so far been 
neglected, and the contributions from F'^^"^ to the flow 
of F^^^A are usually restricted to self-energy corrections 
(see below) or discarded completely. In particular, the 
analysis of competing instabilities (see Sec. Ill) is based 
entirely on a level-2 truncation given by the flow equa- 
tion (|52p for the two-particle vertex, with F^^^^ replaced 
by zero, where in addition the self-energy feedback is 
neglected. This seemingly simple approximation cap- 
tures the complex interplay of fluctuations in the particle- 
particle and particle-hole channel, which leads to inter- 
esting effects such as the generation of d-wave super- 
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conductivity from antiferromagnetic fluctuations. In the 
quantum transport plienomena reviewed in Sec. VI, the 
self-energy as given by the flow equation ([50| plays a 
crucial role. Some of the phenomena described there are 
already obtained by a level- 1 approximation where the 
flowing two-particle vertex in Eq. (|50p is approximated by 
the bare one. That truncation might look like a Hartrce- 
Fock approximation, but it is in fact very different, and it 
works well in cases where Hartrce-Fock fails completely. 

The truncated flow equations are still rather compli- 
cated. They involve the flow of functions, not just a lim- 
ited number of running couplings. For example, the ef- 
fective two-particle interaction in a translation invariant 
system is a function of three independent momentum and 
energy variables. Hence, a simplified parametrization of 
effective interactions is necessary even for a numerical so- 
lution. A useful strategy is to neglect dependences which 
become irrelevant in the low-energy limit, that is, whose 
contributions to the fiow scale to zero. 

Contributions to the effective action are called "rele- 
vant", "marginal", and "irrelevant", if their importance 
increases, stays fixed, or decreases, respectively, upon 
lowering the scale A. This classification can be ob- 
tained from power counting. To this end, one tradi- 
tionally considers a renormalization group transforma- 
tion where one rescales momenta and fields after the inte- 
gration over fields in a momentum shell of width dA such 
that a certain quadratic part of the action (the Gaussia n 
fixed point) remains invariant ( Wilson and Kogutlll974l) . 
From the behavior of the other terms of the action under 
this transformation one can assess directly whether they 
increase, remain invariant, or decrease compared to the 
quadratic part. 

For Fermi systems in dimensions d > I the conven- 
tional RG transformation is not applicable, since the re- 
duction of momentum space by the mode elimination 
cannot be compensate d by a linear rescaling of momenta 
(jShankad . ll99lL [ 19941 ). However, one can perform the 
power counting more directly by estimating the scale de- 
pendences of Feynman diagrams on the right hand side 
of the fiow equations. As described in Sec. II. E. 3, this 
can be done rigorously and to all orders. At the crud- 
est level the power counting is independent of dimen- 
sionality and corresponds to what one would get from 
the above-mentioned RG transformat ion applied to one- 
dimensional systems ( Shankarl . Il994[) . that is: (i) the 
self-energy has a relevant piece describing a Fermi sur- 
face shift, while linear dependences on frequency and mo- 
mentum perpendicular to the Fermi surface are marginal; 

(ii) a regular two-particle interaction is marginal; its de- 
pendences on frequencies and momenta perpendicular 
to the Fermi surface are irrelevant, such that one can 
parametrize it by its static value on the Fermi surface; 

(iii) regular m-particle interactions with to > 3 are ir- 
relevant. This basic classification does not depend on 
dimensionality because the bare propagator Gq, Eq. ([3]), 
is singular on a (d— l)-dimensional surface, such that the 
codimcnsion of the singularity in the {d+1) dimensional 



space spanned by momentum and frequency is always 
two. 

One should, however, not jump to the conclusion that 
the TO > 3 terms can simply be discarded from the RG 
hierarchy in general. This is because effective interac- 
tions with TO > 3 may diverge for small A even in case 
of finite two-particle interactions. For example, the first 
contribution to the fiow of F*^^)^ in Fig. 3 generates a 
three-particle interaction of order A~^ if the external mo- 
menta add up to zero at each vertex. When inserted into 
the equation for F^*^^ in Fig. 3, this may give rise to 
a marginal term of third order in F^^^^. For d ~ 1, 
this term is indeed marginal. However, if d > 2 and 
the Fermi surface is curved, this and other contribu- 
tions are suppressed below the basic power counting es- 
timate due to geometrica. l ly reduced integratio n volumes 
( Feldman and Trubowitd . Il990l: IShankad . [l99l . This im- 
proved power counting is described in App. B.3. It can 
also be used to give a precise, scale-dependent meaning 
to nesting of the Fermi surface. 

A less obvious effect is that this improvement becomes 
uniform, that is, independent of t he external momenta , 
in graphs with ov erlapping loops ( Feldman et aU Il996t 
[SaShofe3,[i998i), so that their contribution gets further 
suppressed (see also App. B.3). It is this effect which im- 
plies that the derivative of the self-energy is not marginal, 
but irrelevant for curved Fermi surfaces in dimension 
d > 2. Moreover, it justifies truncated flows beyond the 
weak-coupling regime, as follows. Consider again the first 
contribution to the flow of F^^^^, shown also in Fig. [5] 
(a). When this term is inserted in the equation for F'^^^, 
the two lines can be joined in two ways, shown in Fig. 
[5] (b) and (c). The graph in (b) gets no extra small 
factor, but the graph in (c) has overlapping loops and 
for positively curved Fermi surfaces in o? = 2, its con- 
tribution gets suppressed by an additional small factor 
~ log ^ at scales below a scale A/, which depends 
only on the geometry of t he constant energy sur faces of 
the initial dispersion Ck ( Feldman et all Il998bl ). This 
suppression holds uniformly for all values of the external 
momenta. For c? = 3, a similar bound holds, without 
the logarit hm. Similar (in . general, weaker) estimates are 
shown in ( Feldman et aZ.I . Il996l ) for general non- nested 



regular Fermi surfaces in d > 2 and for Fermi surfaces 
with V an Hove singularities in ( Feldman and Salmhoferl . 
I2008allbl ). The contribution from graph (c) remains small 
compared to the second-order term if [F^"*)^! log ^ is 

small. Note that this condition does not require |r'^^''^| it- 
self to be small: curvature effects justify dropping these 
terms beyond the weak-coupling regime, provided that 
the above condition is satisfied. This will be used in 
Sec. mil The detailed argument and a discussion of the 
consequences for th e functional RG flow, are given in 
Sections 1 and 5 of ( Salmhofer and Honerkarnpl . l200lt ). 

In the theory of interacting Fermi systems, one is not 
only interested in low-energy fixed points and scaling, but 
also in the behavior at intermediate scales, and formally 
irrelevant terms may play an important role. There are 
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FIG. 5 (a) Third order graph contributmg to T'^''^. (b) Tad- 
pole contraction, (c) Contraction to form a graph with over- 
lapping loops. 



cases where one would like to know the full temperature, 
momentum, or frequency dependence of physical quanti- 
ties, because a low-energy expansion contains insufficient 
information. One of the advantages of the functional RG 
framework is that such dependences can be computed 
directly. 

In many situations, a comparison to standard rcsum- 
mations of the perturbation expansion is desirable, and 
it is also an interesting question to what extent such 
resummations can be reproduced by truncations of the 
functional RG flow equations. A very im portant obser- 
vation regarding this was made by iKatanin. (.200 j ). who 
showed that a partial inclusion of the six-point vertex in 
the flow allows to recover approximations of the type 
Hartree-plus-ladder summations (in cases where these 
approximations are a good starting point). This also al- 
lows to continue fermionic flows into symmetry-broken 
phases (see Sec. lIVp . If we drop the eight-point vertex 
from the equation for F^^^^, it is determined by a Fcyn- 
man graph containing three four-point vertices, depicted 
in Fig. [S^. When backsubstitutcd into the equation for 
the four-point vertex, two external legs get contracted in 
all possible ways. We have just discussed that the contri- 
bution from the graph in (c) is suppressed by improved 
power counting. When two legs of a single-four-point 
vertex are contracted to form a tadpole (see Fig. [S})), 
the value of the thus obtained subgraph is S^, by the 
flow equation for the self-energy. Thus a factor G^S^G^ 
appears in the integral for the value of the graph. By 
Dyson's equation. 



G^ = G^t^G 



(54) 



so this can be combined with the second-order contri- 
bution to replace by G"^. If all other effects of 
the six-point function, corresponding to graphs of the 
type shown in Fig. [5]: are dropped, the equation for 
the four-point function gets changed to one where the 
product G^{k)S^{k') + S^{k)G^{k') is replaced with 
G^(fc)G^(fc')+G^(fc)G^(fc') ^ ^ (G^(fc)G^(fc'))- If one 
now restricts further to a single channel in the four-point 
equation, it becomes explicitly solvable by a ladder sum- 
mation in that channel. Backsubstitution in the equa- 
tion for the self-energy gives the c orresponding H artree- 
type term. This is explained in ( Kataninl . |2004[ ). and, 
also in its extension to flows with symmetry breaking, in 
(ISalmhofer et aLl . l2004D . 



D. Flow parameters 

In the derivation of the exact functional flow equation, 
the scale dependence of the bare propagator Gq was not 
specified. The derivation holds for any choice of Gq, 
provided all functions involved are indeed diflterentiable 
with respect to A, and provided that the resulting flow 
equation is well defined. These conditions arc not trivial; 
in fact, badly chosen flow parameters may lead to diver- 
gences on the right hand side of the flow equation. On 
the other hand, one can exploit the flexibility provided 
by the choice of the A-dcpcndcnce to ones own advan- 
tage. Besides regularity issues, the scale dependence of 
Gq is only constrained by the initial condition 



G^" = 



and the final condition 



G, 



A^O 



— Go 



(55) 



(56) 



The functional F = F"^"*" reached at the end of the ex- 
act flow is independent of the choice of Gq. However, 
in most practical calculations, where approximations are 
unavoidable, a judicious choice of Gg is mandatory. Im- 
portant aspects related to the choice of Gq are: regular- 
ization of infrared singularities, minimization of trunca- 
tion errors, respecting symmetries, technical convenience. 
In the following we will review the most frequently used 
cutoff schemes along with their merits and drawbacks. 



1. Momentum and frequency cutoffs 

For the sake of a concise discussion, let us focus on 
translation and spin-rotation invariant one-band systems, 
such that the bare propagator Go can be written as a sim- 
ple function of frequency and momentum as in Eq. ([3|). 
The scale dependence can then be introduced by multi- 
plying Go with a suitable cutoff function 6^ , 



G^(fco,k)==0^(fco,k) Go(A:o,k) , 



(57) 



with 6^" ~ and 6^^*^ ~ 1. To regularize the infrared 
divergence of Go at zero frequency and for momenta on 
the Fermi surface (^k = 0), the cutoff function 6^{kQ,'k) 
has to vanish sufficiently quickly for /cq ^ 0, at 

fixed A > 0. The most frequently used cutoff functions 
are either pure momentum cutoffs of the form 0'^(k) = 
i?(|Ck|/A) or frequency cutoffs 9^{ko) = 7?(|A:o|/A), where 
'd{x) is a function that vanishes for a; ^ 1 and tends to 
one for a; ^ 1. Mixed momentum and frequency cutoffs 
of the form 6*^(^0, k) = d[(k'^ + Ck)/^^] are preferred 
in the mathematical literature, as they facilitate power 
counting and rigorous estimates. 

A technical advantage of momentum cutoffs compared 
to frequency cutoffs is that Matsubara sums on the right 
hand side of the flow equations can often be performed 
analytically. Furthermore, a momentum cutoff does not 
spoil the analytic structure of propagators and vertex 
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functions in the complex frequency plane. However, 
there are also serious drawbacks, which are specific to 
Fermi systems. Once self-energy effects are taken into 
account, the Fermi surface is usually deformed in the 
course of the flow, such that the momentum cutoff has 
to be continuously adapted to the new Fermi surface, 
which complicates the flow equations considerably. Sec- 
ond, particle-hole excitations with a small momentum 
transfer q are suppressed by the momentum cutoff for 
l^k+q ~ Ckl < 2A. As a consequence, the limit of van- 
ishing momentum transfer q — ^ in interaction vertices 
and res ponse functions does n ot commute with the limit 
A -j> (jMetzner et all Il998l ). In other words, forward 
scattering interactions and the response to homogeneous 
fields can be obtained only by taking the limit q — » at 
the en d of the flow, at A = ( Honerkamp and Salmhofen . 
l2001aD . This is a serious drawback in stability analyses 
(see Sec. IIII[) . where one compares the increase of the 
effective interaction in different momentum channels (in- 
cluding forward scattering), or different susceptibilities, 
upon lowering A until a divergence occurs in at least one 
channel at a finite scale Ac > 0. 

A frequency cutoff has the advantage that it does not 
interfere with Fermi surface shifts, and that particle-hole 
processes with small momentum transfers are capture d 
smoothly by the flow ( Husemann and Salmhofeil . |2009() . 
It can al so be used in systems w ithout translation in- 
variance ( Andergassen et all . 1200^ . where a momentum 
cutoff is less useful since the propagator is not diagonal 
in momentum space. However, a frequency cutoff affects 
the analytic properties of propagators and vertex func- 
tions in the complex frequency plane. Depending on the 
sort of truncation used, this may pose a serious prob- 
lem if one likes to continue results to real frequency. For 
a frequency cutoff the initial cutoff is Ag = oo. Since 
the contributions to the self-energy flow are of order 
A~^ at large A, one has to retain the convergence fac- 
on the r i ght h and side of the flow equation 



tor 6'*^°° on the r i ght h and side of the flow cquatior 
( Andergassen et aZI . |2004J ). analogously to the convcr 



gence f actor in the perturbation expansion of the self- 
energy ( Negele and Orland . 19871): for a rigorous justifi- 



cation, see Pedra and Salmhoferl (|2008[ ). 

For a bare propagator Go(fco, k) = (ifco — ^k)~^ and a 
multiphcativc cutoff as in Eq. ([57)). the full propagator 



has the form 



G^(fco,k) 



zfco-ek-0^(fco,k)I]A(fco,k) ' 



(58) 



and the single-scale propagator 
Eq. (gT]), reads 



G^g^G^, see 



(^fco-Ck)^Ag^(fco,k) 
[zfco-Ck-e^(fco,k)SA(fco,k) 



2 • 



(59) 



For a sharp cutoff function such as 9^{ko) = Q{\ko\ — 
A), the single-scale propagator seems ill-defined, since 
i9a0(|^oI — A) = —^(1^0 1 — so that a delta peak in 
the numerator of Eq. (j59p coincides with a discontinuity 



(due to the step function) in the denominator. However, 
this ambiguity can be easily removed by viewing the step 
function &{x) as a limit of increasingly sharp regularized 
step functions 8e(x), where t he discontinuity is smeared 
over a width e (jMorrid . [l993) . With S,{x) = d^Oeix), a 
simple substitution of variables yields 



(60) 



for any continuous function /. Note that the right hand 
side is unique, that is, it docs not depend on the shape of 
the smeared step function Qe{x) for e > 0. For a sharp 
frequency cutoff 9^{ko) = O(|fco| — A), for example, the 
single-scale propagator thus simplifies to 



S'\ko,k) = — 



S{\ko\-A) 



fco-a-SA(fco,k) 



(61) 



as long as it does not appear in products where other 
factors are also discontinuous at |fco| = A. Otherwise, for 
example in products of the form S''^(fco, k)[G^(fco, k)]™, 
one has to apply Eq. (pO)) to the entire product. 

A sharp cutoff has the obvious technical advantage 
that the integration over the cutoff variable (fco or ^k) can 
be carried out analytically, thanks to the delta-function 
in the numerator of S^. On the other hand, a sharp 
cutoff generates discontinuities in the momentum or fre- 
quency dependences of the vertex functions, correspond- 
i ng to a pron ounced non- locality of the effective action 
(|Morrisl . [l99l . which is often not amenable to a simple 
parametrization. At finite temperature the flow equa- 
tions arc ill-defined for a sharp frequency cutoff, since 
the Matsubara frequencies are discrete: fco — {2n + l)7rT 
with integer n. Continuous cutoff functions at T > are 
conveniently chosen such that the A-derivative is non- 
zero only in a frequency range of width 27rT, since then 
only two frequencies contribute to the Mat subara sum 
on th e right hand side of the flow equation ( Enss et~al\ . 
120051 ). 

There are useful cutoff schemes which are formulated 
more naturally by adding a regulator function to the 
inverse propagator (instead of multiplying): 

Q^(fco,k)= [G^(fco,k)]"' =Qo(fco,k)+i?^(fco,k) , 

(62) 

with i?"^" = oo and _R ^^" = 0. In particular, regulator 
functions of the form (iLitiml . I2OO1I) 



R\k) - -Z^[sgn(Ck)A - ek]e(A - 1^1) 



(63) 



or its frequency dependent analogue, i?^(fco) = 
iZ^[sgn{ko)A — ko]<d{A — \ko\), have some distinct ad- 
vantages. The prcfactor is initially one and is then 
determined by a momentum (or frequency) derivative of 
the flowing self-energy T,^{ko, k). The Litim cutoff satis- 
fies a criterion of "optimal" re gularization o f the infrared 
singularity of the propagator (jLitiml . I2OOII ) . For simple 
truncations it also leads to a very convenient form of the 
integrands, facilitating the integrations. 
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It is easy to choose the cutoff function in a way that 
does not affect the global symmetries of the system, such 
as global charge conservation or global spin rotatation 
invariance. However, local conservation laws are typically 
spoiled. The corresponding Ward identities are modified 
by cutoff dependen t additional terms, which vanish only 
in the limit A — > ( Ens^ [20051 ) . It is very hard to devise 
truncations which satisfy the modified Ward identities 
at each scale A, and hence truncated fl ows often violate 
Ward identities also in the limit A (|Kataninl |2004[) . 
In these cases it is better to compute only independent 
quantities from the flow, and determine the remaining 
quantities, which are fixed by local conservation laws, 
via the Ward identity. 



2. Temperature and interaction flows 

For fermion systems the infrared singularity of the 
bare propagator can also be regularized by tempera- 
ture, instead of a cutoff, since the fermionic Matsub- 
ara frequencies stay away from zero at a distance ttT. 
A flow equation with temperature as a flow parame- 
ter can be obtained from the general flow equation de- 
rived in Sec. III.Bl if one manages to shift all temper- 
ature dependences of the bare action to the quadratic 
part. This is indeed possible by a s imple r escaling of the 
fields ( Honerkamp and Salmhofeil l2001a[) . Let us con- 



sider a translation invariant system of spin-i fermions 
for definiteness, where the fields depend on a momen- 
tum k, a spin index a, and a Matsubara frequency 
ujn = (2n -I- l)7rT. Rescaling the fields as ip'^{Ti,'k) ~ 
T^/^^„{uJn: k) and ^'„{n, k) = T3/4^^(tj„, k) removes aU 
explicit T-factors from the (quartic) interaction in the 
bare action. The temperature dependence is thereby 
shifted entirely to the quadratic part of the action, given 
by the inverse bare propagator for the rescaled fields. 



(64) 



The effective action P-'" [?/»', -0'] for the rescaled fields 
obeys the exact flow equation Eq. ([35]) , with temperature 
as the flow parameter. The unsealed vertex functions 
p(2m) recovered from the vertex functions r'^™)-^ by 
multiplying with 7^3'"/2_ The temperature flow has sev- 
eral advantageous features. First, it generates directly 
a temperature scan of the computed quantities. In cut- 
off schemes one has to run a full flow for each temper- 
ature separately. Second, the temperature flow includes 
particle-hole excitations with small momentum transfers 
uniformly at each scale. Third, local symmetries and the 
corresponding Ward identities arc respected at each step 
at least for the exact flow, which makes the still difficult 
issue of Ward identities in truncated flows at least more 
transparent. 

A particularly simple choice of a flow parameter is pro- 
vided by a uniform factor A scaling the bare propagator 



(iHonerkamp et a/.l . l2004 ). 

Gq = XGq , 



(65) 



with Aq = 0, and A — > 1 at the end of the flow. By a sim- 
ple rescaling of the fields one can see that this is equiv- 
alent to multiplying the bare quartic interaction with a 
factor A^, which means that the interaction is scaled up 
continuously from to its full strength in the course of the 
flow. Hence the name "interaction fiow" for this scheme. 
In the absence of self-energy feedback the interaction fiow 
has the technical advantage that the propagator has the 
same form at each scale, such that certain loop integrals 
need to be done only once. However, the global scaling of 
the propagator does not regularize the infrared singular- 
ities, such that one easily runs into infrared divergences. 
Nevertheless, for suitable problems and simple trunca- 
tions the interaction flow has been shown to yield results 
simila r to flows with a cutoff, an d with less computational 
effort (jHonerkamp et aLl . l2004) . 



E. General properties of the RG equations 

In this section we discuss the general structure of the 
RG hierarchy of equations and provide power counting 
bounds for its solution. These bounds are simple, but 
mathematically exact, and they provide a strict sense to 
the notion of relevant and irrelevant terms. We shall also 
briefly discuss improved power counting bounds, which 
provide sharper estimates for bulk Fermi systems in d > 2 
and exhibit the role of Fermi surface geometry. 

The generating functionals were introduced to obtain 
the Green functions and vertex functions of the model by 
differentiation, cf. ([7]) and (|10p . In the framework of the 
RG as an iterated convolution, they acquire an indepen- 
dent importance. Indeed, in many situations in bosonic 
field theory, an expansion i n the fields is avoid ed in fa- 
vor of a gradient expansion (|Berges et a/.l . l2002D or other 
types of parametrization (see also Sec. V), and the flow 
may lead to a non-analytic function of the fields. Func- 
tions of Grassmann variables are, however, defined only 
by power series expansions in these variables, so in this 
case the meaning of the RG is strictly that of the infinite 
hierarchy. This is only a seeming disadvantage because 
by the anticommutation properties of Grassmann vari- 
ables, the fully regularized functionals (as they appear 
in the flow equations) have co nvergent expans ions in the 
fields (lAbdesselam and Rivasseau, ,1998; Feldman et al\. 



1998a . 



Lesniewski 



20021 119861: iGawedzki and KupiainenL I198E 
19871: ISalmhofer and Wieczerkowskil I200C 



In contrast, the expansion in the fields of bosonic func- 
tionals is almost always divergent, even in the regular- 
ized theory. Convergent expansions then take the form 
of cluster expansions that distinguish bet ween regions of 
small and large fields (see, for example, iBalaban et al\ 
( 2010l ) and references therein). 
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1. Inductive structure of the RG hierarchy 

The functional appearing in p6)) has an 

expansion similar to (|40)) . namely x) = 

X]m>i ^](a;', a;), where Sf^™)^ is homogeneous 



of degree 2m in the fields, hence has a representation with 
coefficient functions S^^™)^ similar to (HI]). By definition 
(i45|) of S, the S(2™)A are determined by r(2™+2)A^ foj. 
example 



^ (x]^,..., Xi , . . . , X77i) ^ ; Xi , , . . . , Xtti, x) . (66) 

Here the indices refer to the matrix structure of (j45p . The other matrix elements arc given by similar expressions. 
We use this to expand (|46p in homogeneous parts in and and compare coefficients. This gives 



5^ E tr S^I;(2™°)A[^^^]-Qg^S(2"«)^[V',V^] . (67) 



P>2 mo,.-.,mp>l \ g=l 

TTiQ + .-. + mp = ni 



In the sum over p, each of ttiq, . . . ,mp is at least one 
because S only contains field-dependent terms, and mo-l- 
... -I- nip = m must hold since ^(2m)A i^omogeneous of 
degree 2m in the fields ['0, V"]- These two conditions imply 
that p < m and that mq < m — p for all < q < p, so 
that for every given m, the sum only runs over finitely 
many terms. Since the coefficient in ^(2?ti)A p(2m)A 
^(2m)A ^ p(2m+2)A^ Comparing coefficients of powers of 
ip and -0 in (|67l) gives a hierarchy of differential equations 
for the r(2™)A^ labelled by m. We rewrite §7} as 

_d_p(2m)A _ ^Ap(2)n+2)A ^ ^A ^p(4)A^ p(2m)A 
m 

+ ^£^(r«2")A) . (68) 

p=2 

The three summands on the right hand side are obtained 
from the three terms in (|67p . and it is understood that 
both sides are functions of 2m variables xi, . . . , X2m- The 
action of the operator y,^ on r(2™+2)A jg linear, as is 
that of /^^(rf*)^) on r(2™)A^ while is nonhnear in the 
lower-m vertex functions r^^™)^ = T^^'>^, r(2™-2)A^ 
Specifically, the action of Ti^ is given by a tadpole-type 
contraction and summation, the action of /C^(r(^)'^) is 
given by the evaluation of a one-loop diagram formed 
from r(2™)A and the four-point function T^'')^, and 
is given by a sum over one-loop diagrams involving p+1 
vertex functions, each of which has ruq < m. Thus Ti^, 
/(;A^p(4)A^ and £p also depend on A and on the self- 
energy via the propagators and G^. The H^- 
term couples the higher vertex function r(2ni+2)A ^^^^^ 
the equation for r'2™)A, Thus the hierarchy does not 
close among finitely many m, and therefore truncations 
need to be employed to obtain solutions. 



2. Truncated hierarchies and their iterative solution 

If for some mo > 1, the initial vertex functions r'^2m)Ao 
vanish for all m > mo + 1, one may employ the approxi- 
mation of setting r(2'")A ^ r(2m)Ao foj. all m>mo + l. 
That is, all vertices with m > mg 4- 1 remain zero, and 
the (mo + l)-particle vertex is kept fixed at its initial 
value. This level-mo truncation reduces the infinite hier- 
archy to a system of finitely many differential equations 
for (r(2™)'^)„<„i„. The vertex r(2"o+2)Ao enters in the 
equation for r'^2™o)A^ Specifically, in the level-1 trun- 
cation, the two-particle vertex F^^^^ is fixed to its bare 
value r'^'')'^", and the self-energy is the solution of (|50p . 
The level-2 truncation is given by ([5^ . with T^^^^ fixed to 
its initial value F^^^^" (which may vanish), together with 
([501. The term /C^(F(4)A)r(4)A jg quadratic in r(4)A. 

In the level-mo truncation of the hierarchy, with 
mo > 2, and at given T,^ and (r(2"')A)^,^^^^^ ^q. §^ 
for f(2™o)A becomes a linear inhomogcneous differential 
equation for r'^2mo)A^ which can be solved by an operator 
version of the standard method of variation of the con- 
stant: when all sums and integrals corresponding to the 
traces in (j67p are written out, it takes the form of a lin- 
ear integro-differential equation which is, viewed more 
abstractly, a linear ordinary differential equation in a 
suitable space of functions, to which standard techniques 
apply. Together with the initial condition f(2"')Ao^ ^^i^jg 
determines f(2™)A uniquely in terms of (r'^2"» ''^)m'<mo- 
Backsubstitution of this solution into the 'H'^-term for 
the equation for f(2™o-2)a ^^i^^u yields an equation for 
p(2mo-2)A^ which can be solved, in terms of the not 
yet determined lower vertex functions (F(2™ ^^)m'<mo-i- 
Proceeding downwards in m in this way, one can formally 
solve the truncated hierarchy, with the final equation de- 
termining . We write "formally" here because after at 
most two steps of this iteration, the differential equations 
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become nonlinear, so that existence of the solution is typ- 
ically known only for short flow times, and because the 
question of blowup of solutions is rather nontrivial. In- 
deed, we shall see below that blowup generically occurs 
in RG equations if relevant terms have not been taken 
into account. This phenomenon is related to the infrared 
divergences of unrcnormalizcd perturbation theory. The 
major advantage of thc RG method is that the growing 
terms can be identified and studied long before they get 
singular, and then removed by taking into account ap- 
propriately chosen relevant parts in the flowing action. 

Increasing mg to improve the accuracy is then a natu- 
ral strategy for approximation of the true solution; how- 
ever, explicit and numerical calculations can be done only 
for small mo, because the number of variables increases 
rapidly with m. Nevertheless, one can get useful infor- 
mation in the form of bounds for the maximal possible 
value of the vertex functions (or other norms that mea- 
sure their size). This is done in the following section. 

3. Running coupling expansion and power counting 

We turn to the standard situation of a model with 
two-body interactions, where the initial interaction of 
the fermion system is quartic, i.e. r^^™)^" = for all 
m > 3. We also assume that this interaction is short- 
range so that its Fourier transform is bounded (e.g. an 
unscreened Coulomb interaction is long-range). In a per- 
turbative expansion in powers of the initial four-point 
interaction = r*^^)^", the vertices are given by sums 
over irreducible graphs. An irreducible Feynman graph 
formed with r four-legged vertices can have at most 2r 
external legs, so that in order r in that expansion, all 
vertex functions with m > r vanish. 

As we shall now explain, one can solve the RG hierar- 
chy in terms of a similar expansion in the scale- dependent 
four-point function = F'^'*'''^, again by integrating 
the RG hierarchy downwards in scale, but keeping the 
2m-point functions for m > 2 only to a fixed order in 
V^. The equation for itself then becomes an integro- 
differential equation with a power r nonlinearity on the 
right hand side (the equation for remains unchanged) . 
This leads in a natural way to power counting estimates 
for the higher 2m-point functions in terms of the maximal 
value of the four-point vertex that occurs in the flow. 

We denote the O {{V^Y) contribution to F^^")^ by 
-p(^2m)A^ Its scale derivative equals 

_(i_p(2m)A _ ^Ap(2m+2)A _|_ ^Aj--^A-j p(2m)A 

The primed sum runs over all sequences (mo, . . . nrip) and 
all sequences (ro,...,rp) with m^ > 1 and > 1 for 
all 1 < q < p, mo + . . . + mp = m + p, and ro + . . . + 
Tp — r. The solution of the RG hierarchy for an initial 

quartic interaction has the property that F^^™^^ = for 



all m > r. Therefore, for m = r, thc T-L term drops 
out of (|68p . and all remaining terms contain only or 
terms of order at most r — 1 in V^. Thus, given these 
lower-order F's, FJ. can be obtained by integration. 
Then the right hand side of the equation for m = r — 1 
is determined, so F^^'^"^''^ can be determined, and so 
on. Successive backsubstitution then leads to a system of 
equations where ^Tr^™''^ gets contributions from a sum 
of graphs with r vertices of type , where Ao > A' > A, 
and propagators and , and all the intermediate 
scales A' etc. are integrated. Thus the equation becomes 
nonlocal in the flow parameter A but the right hand side 
is known once and J:^ are known. V is given by 
a degree r nonlinear equation, with a similar graphical 
background as discussed above, and by the standard 
self-energy equation ([50)1 . While more restricted than 
the level-m truncation, the running coupling scheme also 
captures effects that cannot be seen in any fixed order of 
bare perturbation theory, such as screening or asymptotic 
freedom of certain coupling functions. 

We use this inductive structure to derive basic power 
counting bounds for thc vertex functions in terms of 
thc flowing four-point function, for a d-dimensional bulk 
fermion system (d > 1). For simplicity, we focus 
on spin-i fermions with translation-invariant action, so 
that we can use Xi = {ki,ai) = (fco,,, k^, (7^), as dis- 
cussed at the beginning of Section III.Al Wc also as- 
sume that the symmetries of the action remain unbro- 
ken. These specific assumptions are for presentation 
only; power counting can be done without them. By 
translation invariance, F^^™''^ ((^i: ^1)1 ■ • ■ ik2m, o'2m)) = 
S{J2ih) rl^™''^(fc, ct), where the delta function forces 
conservation of the spatial momentum k (up to recip- 
rocal lattice vectors) and conservation of the frequency 
variable fcp, and we have introduced the abbreviations 
a = (cTi, . ..a2m) and k= (/ci, . . . , k2m-i)- For A = Ao, 
thc function rf^™-*^ is smooth and bounded because the 
initial interaction is short-range, and this stays so during 
the flow above critical scales. 

Consider the maximal size of the vertex functions, 
||F^'"'^^|| = sup |ff'"^'^(fc,a)|. Then, for m > 3, 

|]p(2m)A|| < ^{2m) ^^r-m+lj^r ^ ^jq^ 

where 7^^™'' is independent of A and /?, 

fA= sup ||y^|| (71) 

A<^<Ao 

is thc maximal value of the four-point coupling on all 
scales between A and Aq, and 

5A = maxE / 5fc|^a,a'(fc)l , (72) 

a' 

with / dk . . . = TJ2ko / {27T)'' ■ ■ ■ general power 

counting, we do not need to assume that the propagator 



17 



is diagonal in the spin indices a, a', so also carries 
these indices) . The dependence of s\ on A is determined 
by the shape of the Fermi surface. As shown in Appendix 
IB. 11 sa is of order one for regular Fermi surfaces. If 
the Fermi surface contains Van Hove points, sa grows 
logarithmically in A for A —> 0. 

At first sight, one may worry about the factor A^~™, 
which diverges for m > 3 as A — > 0. For the maximum 
value of the vertex functions, it is indeed true - and eas- 
ily verified in examples - that there are always particular 
values of the external momenta where these vertex func- 
tions become very large in /3 = 1/T (and diverge at zero 
temperature). However, this happens only on a "small" 
set of momenta. For a general m-point function, it is very 
involved to determine this set, but this is not necessary 
for power counting. One can use the norm instead, i.e. 



consider \\r?"''>\ ^ j:^Jdh...dk^ \rr"'\La) 



Using generalizations of 
if /a remains finite 



TOl). one can then show that 



(73) 



with constants c^^™'' t hat are independ ent of A, /3 and 
the system size L (see lSalmhofe^ (I1999D . Section 4.4.3). 
This implies that, even in the limit /? -> oo, the 2m- 
point vertices can become singular only on a set of zero 
Lebcsgue measure in momentum space. In general, this 
set can be rather complicated, but, loosely speaking, it 
will have codimension at least one. 

It is one of the appealing features of the flow-equation 
RG that exact statements like (|70p can be proven in a 
few lines, see Appendix [BJ The argument given there 
also implies an at-most logarithmic growth of the coeffi- 
cients in the equation for /a itself. The self-energy then 
comes out of order /a, provided that rcnormalization is 
done correctly, see Section FILE. 41 At small /a, the size 
of the vertices Y^^"^")^ with m > 3 is thus determined 
by /a. The terms with m > 3 are the RG-irrelevant 
ones, m = 2 is marginal, and to = 1 is relevant. This 
classification is explained in detail in Appendix IB] In a 
Taylor expansion of the vertex functions in the Matsub- 
ara frequency around zero, and in momentum around the 
Fermi surface, additional small factors arise, which cancel 
the small denominators of the propagators; at the same 
time, the vertex function is replaced by a differentiated 
one. Hence, the flow obtained by projecting the frequen- 
cies to zero and the momenta to the Fermi surface gives 
the dominant contribution for small A. This is expected 
from a simple counting of bare scaling dimensions, and 
can be established more rigorously by power-counting ar- 
guments similar to those used above and in Appendix IB] 
In the case of a curved Fermi surface in d > 2, /a 
indeed stays small in a weakly inte racting system abo ve 
a BCS-like temperature (see, e.g., ISalmhofeil ( 1998bf )). 
indicating the absence of symmetry-breaking. At zero 
temperature. f\ grows as A decreases, and the four-point 
function has singularities at points corresponding to nest- 
ing vectors of the Fermi surface; for details, see Appendix 



IB. 31 This growth of f\ with decreasing scale A is usu- 
ally called the "flow to strong coupling" in RG studies, 
and is described in more detail in Section IIIII A sin- 
gularity of the two-particle vertex in momentum space 
means that the interaction becomes long-range in posi- 
tion space. This is associated with the formation of criti- 
cal fluctuations, and in case of spontaneous breaking of a 
continuous symmetry, with the appearance of Goldstone 
bosons (see Section irv|) . 



4. Self-energy and Fermi surface shift 

The self-energy is important for all effective one- 
particle properties of the system, and it can cause dras- 
tic effects, as compared to the non-interacting fermions. 
Accordingly, in the RG flow, the self-energy is a rele- 
vant term. In absence of symmetry breaking, it modifles 
the inverse propagator to iko — ^ 5](fco, k). The long- 
distance behaviour of the fermionic propagator is deter- 
mined by the behaviour of this function around its zero 
set. A Taylor expansion around ko = gives 



iko-^k - S(fco,k) = 



+ p(fco,k) (74) 



with Z^^ = 1 + i(aoS)(0,k), Zi^^ek = Ck + S(0,k) and 
a Taylor remainder p. If p vanishes faster than linearly 
in fco a-s fco ~^ 0, we thus obtain an effective descrip- 
tion in terms of quasiparticles with dispersion relation ek, 
hence "interacting" Fermi surface {k : ^k + S(0,k) = 0}, 
Fermi velocity Vck, and quasiparticle weight Z^. If Z^ 
is bounded and nonvanishing for all k, the long-distance 
decay of the fermion propagator in position space is the 
same as for the free theory. 

The question whether S is smooth enough for the 
above to hold is nontrivial. In the one-dimensional Lut- 
tingcr model, 9oS(0, kp) diverges, and the small-fco be- 
haviour of the self-energy, E(fco,^F) ^ |^o|" with a < 1 
depending on the interaction strength, implies that the 
self-energy effects dominate at small ko and the decay 
in position space becomes more rapid, so that the occu- 
pation number density n(k) becomes a continuous func - 
tion of k even at zero temperature ( GiamarchH . |2004[ ). 
An even more drastic change is spontaneous symmetry 
breaking, where the propagator cannot be written any 
more in the simple form given above (see Section |IV[) . 

In the RG flow, S is replaced by the A-dependent 
self-energy ~ (Gq )^^ — {G^)~^ . An important phe- 
nomenon in this context is the shift in the Fermi surface 
entailed by T,^. In terms of power counting, this shift 
is the most relevant term. Cutting off the propagator 
around the free Fermi surface then fails to regularize the 
propagator, which leads to spurious singularities in the 
RG flow. A convenient method to avoid this is to intro- 
duce a counterterm. In the context of the bare pertur- 
bation expan sion, the count erterm method was already 
described in (|Nozieresl[l96l . The main idea is to antic- 
ipate the form that the propagator takes at the end of 
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the flow and to rearrange the flow such that this form, 
not the bare one is used as the starting point for the 
RG analysis, hence the Fermi surface is fixed to that of 
the interacting system in the flow. The difference be- 
tween the two dispersion functions appears as a (finite) 
counterterm. To obtain a one-to-one relation between 
the model given by the Hamiltonian and the one with 
fixed interacting Fermi surface, one has to solve a self - 
consistency equation. iFeldman and Trubowit j (Il990l ) 
used the counterterm method for the r adius shift of a 
circular Ferm i surf a ce in a RG flow. iFeldman et al\ 
(|l996l . Il998bl Il999l I2OOOI ) generalized this to the case 
of non-circular curved Fermi surfaces, solved the self- 
consistency equation, and showed that remains fi- 
nite for d > 2 to all orders. The corresponding fixed- 
po int problem for the Fermi sur face was also considered 
bv iLedowski and Kopietd ( 20031). The role of Van Hove 
singularities was analyzed by Feldmaii and Salmhofeij 
|2(308a,,b). Feldma n and Trubowitz ( 199l[) also used the 
counterterm method to derive the equation for the su- 
perconducting gap from an RG flow; for further work 
in that direction, see also Section IIVI In alternative to 
counterterms, one can try to avoid a momentum spa ce 
cutoff altogether (iHone rkam p and Salmhoferl . l2001allbl : 
iHusemann and Salmhof cr. 200^ or to use an ada ptive 
scheme (see the appendix in, Honcrk amp et al\ ( 200l[ ) and 
iBenfatto et ail (l2006^ : ISalmhofcr. (,20071) '). 



a charge invariant bilinear is of the form 



F. Flow equations for observables and correlation functions 

All observables of the fermionic system are given by 
polynomials in the fields, so they can be calculated 
from the connected Green functions G*-^™^^, hence by 
the above-mentioned tree relations also from the irre- 
ducible vertex functions r'^™^"^. It is nevertheless conve- 
nient, and due to the limitations of approximations often 
mandatory, to calculate the flow of observables and their 
correlation (or response) functions by separate flow equa- 
tions, which we derive and discuss now. 

For simplicity we restrict the presentation to the par- 
ticularly important class of observables that are corre- 
lations of fermionic bilinears. Charge-invariant bilinears 
are of the form 



y,y' 

Charge-non-invariant bilinears are of the form 



(75) 



y,y' 



B{x) = Y] ( V(2/) B{x; y, y')i>{y') + i>{y) B{x; y, y')i:{y': 

(76) 

The functions B and B determine the spatial and spin 
structure of these bilinears. For translation-invariant sys- 
tems, we can choose a momentum representation where 
X = (fco, k) and y = {po, p, a), as explained above Eq. ([3]). 
With the notations p = {po,p) and J dp = Tj^p^ J dp, 



B{k) 



dp Ipaip) K^a' {P, k) -0<T' (p+k). 



(77) 



The frequency fco is an integer multiple of 27rT. The case 
k = and b — ai, where denotes the i'^ Pauli ma- 
trix, corresponds to a uniform spin density. The case 
fco = 0, k = (tt, TT, . . . , tt) and h = Ui corresponds to 
a staggered spin density. Similarly, the same choices of 
fc but with ba,cr' = Sa-.a' corrcspond to charge densities. 
Other choices of k can be used to test tendencies towards 
noncommensurate magnetic or charge ordering. Cooper 
pair fields correspond to the non-charge invariant combi- 
nations 



B{k) 



dp [ Ipa- (P) ^a,cr' {p, k) {^P + k) 



+ ^/v(p) A^',a(-p + fc,fc)V'<T'(-p + fc)] . (78) 



Again, the simplest choice is uniform singlet pairing, 
where fc = and Aa,cr'{p,k) = A(p)£cr,CT'- In this case 
A(p) is the gap function. Triplet pairing, extended 
Cooper pairs, and spatially nonuniform gaps, are de- 
scribed by suitable generalizations. 

A convenient way of generating correlation functions 
of the bilinears B is to couple the B{x) to external source 
fields J{x), i.e. to add a term (J, S) = J2x J[x)B{x) to 
the action. The external field J is not an integration 
variable, so it can be regarded as a (functional) parameter 
on which Q depends. Writing Q — Q{J,r],fj), we then 
have 



{B{x)B{y))^{B{x)){B{y)) = 



dJ{x)dJ{y) 



.7=0 

TJ,T?=0 



(79) 

In presence of J, the effective action F = r( J, i/j, -0), as 
well as all other quantities appearing in the fermionic 
Legendre transform (|9]), depend on J as well. Since re- 
lations such as ^(J, t/),-!/)) = ri{J,ip,'ip) remain valid for 
any J, straightforward differentiation yields 



dJ{x)dJ{y) 



,7 = 
T7 ,77=0 



g'r(j,^,v;) 

dJ{x)dJ{y) 



.7=0 



(80) 



Graphically, this relation is intuitive in that the bilinears 
always couple to the (effective) vertices by two lines, and 
the fermionic vertices are all even, so that the graphs 
that contribute are automatically irreducible. 

Again, because J plays the role of a parameter, the flow 
equation psp is unchanged. Flow equations for the re- 
sponse functions arc then obtained simply by expanding 
F in the fields J and comparing coefficients. This again 
leads to a hierarchy of equations for the vertex functions 
p(2m,n) fj^g^f^ have 2m fermionic and n external bosonic 
lines. The J-independent term corresponds to the stan- 
dard fermionic hierarchy for the r^^™'°^ = F^^'"^ which 
therefore remains unchanged. When counting powers 
in the fermionic fields, each J corresponds to a bilin- 
ear, so that the truncation F^^'"' = for m > toq 
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for the fermionic vertices corresponds to a truncation 
p(2m,n) = for m + n > mg. The flow equations re- 
maining after a truncation for mg = 3 are shown dia- 
grammatically in Figure [51 



o = --<0 

= --<m: 



--0-- 



FIG. 6 The truncation of the hierarchy for the response func- 
tion tliat corresponds to keeping only the irreducible two- 
particle vertex in the fermionic hierarchy. 

Clearly, the two-point correlation {B{x) B{y)) — 
{B{x)) {B{y)) of any fermionic bilinear B only involves 
the fermionic two- and four-point functions, hence could 
simply be calculated from the knowledge of F^^^ and F'"*^ 
by the reciprocity relation and . Eq. ([50)1 shows that 
the route via external fields in the one-particle irreducible 
equations is strictly equivalent to this if the hierarchy is 
treated exactly. When making truncations to the hier- 
archy and other approximations, the two are no longer 
the same. Anomalous scaling dimensions of fermionic 
bilinears (or other composite objects) are captured eas- 
ily by separate flow equations for these quantities, while 
they are hard to obtain from F^^^ and F'^^\ if the latter 
are computed from a truncated flow equation. An in- 
structive example is given by the calculation of the den- 
sity profile near a s t atic i mpurity in a Luttinger liquid in 
lAndergassen et all (|2004D . 



G. Flow equations for coupled boson-fermlon systems 

The focus of this review is on fermion systems. How- 
ever, even if the bare action involves only fermionic fields, 
bosonic degrees of freedom are frequently generated as 
fermion composites and order parameter fields. For ex- 
ample. Cooper pairs and the order parameter in a su- 
perconductor are bosons. Technically, bosonic fields are 
introduced in an originally purely fermionic theory by a 
Hubbard-Strato novich decou pling of an interaction be- 
tween fcrmions (lPoT)ovl . ll987D . Often the fermionic fields 
are subsequently intregrated out, such that an effective 
action involving only bosons remains. Otherwise one has 
to deal with a coupled theory of fcrmions and bosons. 
In this section we generalize the flow equations derived 
in Sec. III.BI to interacting boson-fermion systems. Flow 
equations for coupled boson-fermions systems have been 
derived b y various groups, with slight differences in the 
notation ( Berges et al\ . |200"3 iKopietz et al\ . 



We first introduce some notation for bosons, and write 
down the bosonic analogues of some of the most impor- 
tant equations from Sec. IIFBI Bosonic particles are de- 
scribed by complex fields 0. It is convenient to combine 
(j) and its complex conjugate 0* in a bosonic Nambu field 



ct>*{x) 



(81) 



The generating fu nctional for connected Gre en functions 
can be written as ( Negele and Orlandl . Il987l) 



g[H] = - In y"l?$e-^[*le(^*'*) , (82) 
where iS[$] is the bare action, and 



H{x) 



h{x) 
h*{x) 



(83) 



the source field. Connected Green functions are obtained 
as functional derivatives 

^2;! , . . . , Xm ,Xi, . . . , a;„j) = 

{(f>{xi) . . . (p{Xm)(p* ix'„^) . ..(p*{x\))c = 



(84) 



H=0 



dh*{xi) ...dh* {x„,)dh{x'„-^) . . . dh{x[) 

The effective action is defined as Legendre transform 

F[$] = (i/*,$)+g[i?], (85) 

with $ = —dQ/dH*. Functional derivatives of F[$] yield 
the bosonic m-particlc vertex functions 



F^ ^ (xi , . . . , Xjm 2;]^ , . . . , a;„j ) 



$=0 



90* (xi ) . . . dcjy* {xra)d<P{x'^ )...d(jy{x'^) 
The matrices of second derivatives at finite fields 



(86) 



dH*{x)dH{x') 
{4>{x)4>*{x')) (0(x)0(x')) 

{^*{x)cj>*{x')) (r(x)0(x')) 



(87) 



and 



r(2)[$] 



a^F 



9$*(x)a$(x') 



obey the reciprocity relation r(^)[$] = (G^^) [iJ])~^. 

Endowing the bare propagator Gq with a cutoff or an- 
other scale dependence, one can derive exact flow equa- 
tions for the generating functionals in complete analogy 
to the fermionic case. In particular, the flow equation for 
the effective action F^ [$] has the form 



1 



.{2)A 



(89) 
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where = diag(g^,Q^*) with = {G^)-\ Note 
that the first term on the right hand side can also be 
written as (0*, Qq4>). The above flow equation is equiva- 
lent to the freq uently used flow equ ation for the effective 
average action (jBerges et all |2002() 



r J? 



with ^ Qo - Qo, which reads (jWetterichl . \l99. 



4-r^[$] = -tr 

dA ^ 2 



,(2)A 



(90) 



(91) 



Order parameters are often associated with real (not 
complex) bosonic fields. In that (simpler) case the 
above equations are still valid if one replaces the com- 
plex Nambu fields $ and H by the real fields (j> and h. 

A generalization to coupled fermion-boson systems is 
now straightforward. Bosonic and fermionic fields are 
conventiently collected in a "super- field" 



(92) 



where $ and 5* are the bosonic and fermionic Nambu 
fields defined above (see Sec. II. A). The conjugate super- 
field is given by 



(93) 



The generating functional for connected Green functions 
involving both bosons and fermions reads 



(94) 

where is the bare action, and Hi, and Hf are 

the Nambu source fields for bosons and fermions, respec- 
tively. Functional derivatives with respect to the source 
fields generate connected Green functions with an arbi- 
trary number of bosonic and fermionic fields, the only 
general constraint being that the number of fermion fields 
is always even. 

The effective action r[$,^'] is given by the Legendre 
transform 

r[$,*] = (if^$) + (i//,*) + g(ifb,i/^) , (95) 

where $ = -dG/dH* and * = -dG/dHf. The source 
fields may also be collected in a super-field 



n = 



Hf 



(96) 



The Legendre transform can then be written more con- 
cisely as r[s] = {n,E) + G{n). 

The matrix of second functional derivatives of G at 
finite fields 

d^G 



dn{x)dn{x') 

($(a;)$*(x')) 
(*(x)$*(x')) 



($(a;)*(x')) 



involves also mixed boson-fermion propagators, which 
vanish only for H ~ 0. The matrix of second derivatives 
of the effective action 



r(2)[s] 



dE{x)E{x') 

( 



9$*(a;)9$(x') d^*{x)d^{x') 



(98) 



\ d^{x)d^{x') d^[x)d'^{x') ) 



is related to G'^^ \H\ by the reciprocity relation F^^^ [S] = 
(G(2)[H])-i. 

A flow of the generating functionals is generated by 
modifing the bare propagators for bosons and fermions, 
Gbo and G/o, such that they depend on some scale pa- 
rameter A. We denote the scale dependent bare propaga- 
tors by G^Q and G^q, and their inverse by Q^o ^^'^ Qfo- 
The generalization of the exact flow equations for the ef- 
fective action in purely bosonic or fermionic systems to 
coupled boson-fermion systems reads 



-(cI.*,QA<i>)--(^,QA^M/) 



1 
2 

-Sty 
2 



where 







fo 



(99) 



(100) 



The supertrace Str incorporates a minus sign in the 
fermionic sector. T he flow equat i on (|9 9l) is equivalent 
to the flow equation ( Berges et all |2002|) 



^rA[.,v.]4str 



R^ 



(2)A 



$, + R^ 



(101) 



with R'^ = Qq — Qo, for the effective average action 

r^[<i>,vI/]^rA[$,vI/]_i(<I,*,RA$) + i(^r,RAvI,). (102) 

The expansion of the exact functional flow equation 
dnni) proceeds in complete analogy to the purely fermionic 
case. Inserting 



.(2)Ar 



[$,*] = (G^)-i-SA[$,^-], 



(103) 



with 



qA ^ /'r(2)A[$^^] 



*=*=0 



diag(G^G^^ 



(97) into the functional flow equation one obtains 



(104) 
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-Sti 
2 



(105) 



with the single-scale propagator 



diag(S^S 



A .A) 



(106) 



The expansion in powers of the fields is now straight- 
forward and leads to a hierarchy of flow equations for all 
vertex functions. The first few terms are shown diagram- 
matically in Fig. [71 



dA-^f o w 



dA^" 



O 

d W / 
dA 



+ -^P — ^ + + j5( 



dA 



FIG. 7 Diagrammatic representation of the flow equations 
for the (fermionic and bosonic) self-energies and some of 
the interaction vertices in a coupled boson-fermion theory. 
Solid lines denote fermionic, dashed lines bosonic propaga- 
tors. Propagators with a dash are single-scale propagators. 



The flow equations derived above are also valid in 
case of C/(l)-synimetry breaking, if one allows for off- 

diag onal elements in the matrices Q^^, Q% , 

etc. (|Berges et aZ.1 . 120021: [Schiitz and Kopietzl l2006l) 

Coupled flow equations for fermions and bosonic 
Hubbard-Stratonovich flelds are particularly convenient 
to treat fluctuations associated with spontaneous symme- 
try breaking (see Sec. IV) and quantum criticality (see 
Sec. V), but they may also be used to study Luttinger 
liquids a nd other symmetric sta t es in interacting Fermi 
systern s (iBartosch et all 2009at iLedowski and Kopietzj . 
l2007ai rbl: ISchiitz et all homi . 



III. COMPETING INSTABILITIES 

In this section we describe how one can apply the level- 
2 truncation of the fermionic RG, mainly without self- 
energy corrections, to two-dimensional fermion systems. 



and study the interplay of ordering tendencies. In re- 
summations of perturbation theory, their manifestation 
are singularities in the four-point function and in certain 
susceptibilities. In the RG, the precursor to a singularity 
is the growth of some parts of the vertex function (often 
termed "flow to strong coupling"). Since singularities in 
the vertex function change the power counting drasti- 
cally, this truncated flow then has to be stopped before 
a singularity happens, at a scale A* > 0, where one can 
read off the dominant interactions and infer a tentative 
phase diagram (in this, susceptibilities are used to com- 
pare the strength of different ordering tendencies and to 
determine A*). As discussed in Sec. III. C. 21 curvature 
effects of the Fermi surface imply that the truncations 
discussed here can be used also when the interaction is 
no longer small, provided that the power counting im- 
provement factor times the interaction strength remains 
small. To obtain a true phase diagram, however, one 
needs to integrate over all degrees of freedom, also those 
with scales below A*. This has been achieved in some 
cases (see Section HV)) . but much remains to be done. 

This first step of monitoring the fiow to strong cou- 
pling above A* , as described in this section, is important 
for the following reasons. (1) It allows to determine the 
effective interaction just above transition scales from the 
given microscopic model without any additional a pri- 
ori assumptions about the nature of symmetry-breaking, 
and thereby provides an initial condition for the integra- 
tion at scales below A*. (2) It exhibits how the interplay 
of the scale-dependent scattering processes on different 
parts of the Fermi surface gradually builds up the effec- 
tive interaction. (3) It has by now become a versatile 
tool for analyzing models with an elaborate microscopic 
structure, such as multiple bands. 



A. Hubbard model and A'^-patch RG schemes 

The Hubbard model and its extensions have become 
standard in correlated fermion systems: on the square 
lattice as a candidate model for high-temperature su- 
perconducting cuprates (And e rson . ,1997: Fulde, 1991i), 
in a multiband-gener alization, for the ne wly discovered 
iron superconductors ( Mivake et all 2010l). on triangular - 
lattices for orga nic crystals ( Kino and Fukuvainal Il996l : 
iMcKenzie ^ 19971). on the honeycomb lattic e for graphene 
(iHerbutl . 120061: iLopez-Sancho all l2009l) . The Hamil- 
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tonian for the simplest one-band Hubbard model reads 



(107) 



analyses with momentum- dependent vertices were 
done using the Polchinski (IZanchi and Schulz . 1997 , 
19981. l2000l i. the Wick ordered (iHalboth and Metznei 



2000a 



b ), and the o ne-pa rticle irreducible flow equations 



(|Honerkamp et all l200l[) . all with a momentum space 



where U-j = is the hopping amplitude between sites regulator. To include ferr omagnetism, the tempera 



« and j and U is the Hubbard on-site repulsion. We 
consider here mainly the case with only nearest-neighbor 
hopping t and next-to-ncarcst neighbor hopping t' on a 
square lattice. Additional hopping terms can be added 
if a more detailed description of the band structure is 
required, and other interaction terms may be added. The 
chemical potential fi and t and t' determine the band 
structure = — 2i(cos + cos ky) — At' cos k^ cos ky — /i, 
and hence the shape of the Fermi surface. 

Resummations of perturbation theory in U suggest 
singularities in different channels, arising fro m Fermi 
surfa ce nesting and Van Hove singularities (jSchuld . 
I1987D . hence competing effects, which are best treated 
by RG methods. After two-patch studies, which pro- 
vided a very crude approximation to the momentum 
depen d ence of the fo ur -point vertex (iDzvaloshinskii , 



19871: iFurukawa et all 



Lederer et 



199a 



■ 19981 : [Gonzalez et all 
19871: ISchuU Il987l) . more careful 



ture flo w was introduced by 'Ho nerkamp and Salmhofeil 
(|2001al lbl) and jl oncrkamp (20 pi|), and further devel- 
oped by iKatanin and Kampj (|2003| ). The results of 
these studies at Van Hove filling were confirmed using 
a refined parametrization of the wavevector dependence 
( Husemann and Salmhofeil |2009() . The decoupling of the 
various ordering tendencies in the limit of small U very 
close to the instability and t he influence of non-loca l in- 
teractions were discussed bv lBinz et al ] (|2002l . l200,'^ . 



In the general RG setup of SectionHH the fermion fields 
now carry a spin index s and a multiindex K consisting 
of Matsubara frequencies w, wavevectors k, and possibly 
a band index b. To avoid bias, the action is required to 
reta in all symmetries of the initial action. This implies 
fsee | Honerkamp et all ( 200ll ): [Salmhofer and HonerkampI 
(l200l|)) that 



(108) 



for a spin- rotation invariant system. By lattice- and time-translation invariance, if 4 is fixed by Jv 1,7^2 and in the 
one-band model (in multiband models, the fourth band index 64 still remains free). We therefore abbreviate notation 
to V^{Ki, K2, K:}). In the truncation F^^'''^ = 0, the flow equations for the self-energy and for the coupling function 
become 



^^\K) = - j dK' [2V\K,K',K)-V\K,K',K')\ S\K') , = T^p + T^H,d + T^H^cr 

with the particle-particle term Tpp and the direct and crossed particle- hole terms Tp^j^ and 7^^^,^: 



(109) 



dKV^{Ki,K2,K) L^{K,-K + Ki+K2)V^{K,-K + Ki+K2,K3) , (110) 



dK 



^2V^{Ki,K,K3) L^{K,K + Ki- K3)V''{K + Ki - K3,K2,K) 



+V^{Ki,K, K + Ki- K3) L^{K, K + Ki^ K3) V^{K + Ki - K3, K2, K) 



+V^{Ki,K, K3) L^{K, K + Ki- K3) V^{K2,K + Ki - K3, K) 
T^H,cr{Ki,K2; K3, Ki) = I dK V^{Ki,K + K2- K3, K) L^{K, K + K2 - K3) V^{K, i^a, i^a) 



(111) 
(112) 

Here L^{K,K') ^ S^{K)G^{K') + G^{K)S^{K') is the product of single-scale propagators and fuU propagators 
with momentum assignments corresponding to the diagr ams in Fig. [51 



For the Hubbard Hamiltonian (|107p . the initial condi- bare coupling, since a contribution to F*^^)"^ is generated 

tion is {Ki, K2, K3) = U. Other interactions can be at third order in the two-particle interaction, which leads 

dealt with by modifying this initial condition. The trun- to third order contributions to the flow of (see Sec. II). 

cation F'^''^ = is justifled only for a sufficiently small In most studies the self-energy feedback into the fiow of 
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was also neglected, since it also affects the flow only 
at third order in V^. 

The coupling function V^{Ki,K2,K3) depends on 
three wavevectors and three Matsubara frequencies, so 
that the RG equation for a two-dimensional system is a 
differential equation in a 9-dimensional space. As dis- 
cussed in Section III.El its most singular part sits at 
zero Matsubara frequency. Hence one may neglect the 
frequency dependence. Then defines an effective 
Hamiltonian. Similarly, the k-dependence is most im- 
portant in the angular direction along the Fermi sur- 
face. This dependence can then be taken into account 
by a discretization, i.e. by devising patches in the 
Brillouin z one in which t he co upling function is kept 
constant. iFeldman et al. I (I1992I) showed that using N 
patche s leads to a natural TV- vector mo del in two dimen- 
sions. IZanchi and Schuld (Il998l |2000D were the first to 
use it in studies of the Hubbard model. 

Usually one forms elongated patches that extend 
roughly perpendicular to the Fermi surface but are rather 
narrow parallel to the Fermi surface (sec Fig. |9]). The 
coupling function is then computed for wavevectors ki 
to k3 at the Fermi surface in the center of the patches. 
We label the patches by = 1,...A^. The function 
is thus approximated by 0{N^) interpatch couplings 
V^{ki, K2, K3). Even if ki,k2 and ks are on the Fermi 
surface, k4 can be anywhere. In the calculation of 
the loop integrals it is however necessary to assign a 
patch number K4 to k4, which amounts to an approx- 
imation of projecting k4 on the Fermi surface. Note 
that this projected A^-patch discretized coupling function 
V^{ki, K2, K3) then has fewer symmetries; for instance 
V^{ki, K2, K3) 7^ V^{k2, Ki, K4) in general, as in the lat- 
ter object ka is not necessarily on the Fermi surface. For 
sufficiently large N, this discretization captures the an- 
gular variation of the coupling function along the Fermi 
surface with good precision. 

The results obtained within this approximation, 
described in the following, have been found to 

be robust when t he d e pendence on frequen cies uji 

(iHonerkamp et all 120071: iKlironomos and Tsail . l2006l ) 
and t he component of k^ transvers a l to the Fermi sur- 
face fHalboth and Mctzncrj, l2000at iHoncrkamp^. '200]1 
IHonerkamp et al.. . .2004) are included. [Katanin (20091) 
performed a flow to third order in the scale-dependent 
four-point- vertex (see Section IH.E.Sp . with the fre- 
quency dependence in the sa. me approximation as 
IHonerkamp and Salmhofeil ( 20031 ). 



a. /v'2 . 



B. Results for the two-dimensional Hubbard model 

Starting from the initial condition given by the Hub- 
bard model, the flow is run from Aq down to a charac- 
teristic scale A, , where the largest coupling reaches some 
multiple a of the bandwidth. The choice of a varies 
widely in the literature; the discussion here is based on 
the comparably cautious choice a = 2 or 3, as well as 



e::, 1 





FIG. 8 Top row: The coupling function V^{Ki, K2, K3) with 
the spin convention, and the diagrams entering in the flow 
equation for the self-energy (middle and right diagram). Mid- 
dle and bottom row: The diagrams for the flow of the couphng 
function. The internal lines are either full propagators or 
single-scale propagators S^. 




FIG. 9 (Color online) A'^-patch discretization of the Brillouin 
zone for the one-band Hubbard model on the 2D square lat- 
tice. The colored region is a patch in which the coupling 
function is approximated as a constant. 



on the consistency check that the results do not change 
drastically as a is changed. The characteristic scale A* 
corresponds to a temperature T*. If T is clearly above 
T*, the fiow can be integrated to scale zero without any 
instabilities. is only an upper bound for the tempera- 
ture where ordering can set in because of order parameter 
fiuctuations at scales below A*. In two dimensions they 
are so strong that long-range order that breaks continu- 
ous symmetries does not occur at any T > 0, thus "or- 
dering" is to mean either short-range order with a very 
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large correlation length, or ordering in a related system 
with a small coupling in the third direction, as is present 
in most materials. 



1. Antiferromagnetism and Superconductivity 

The results discussed here are obtained with a slightly 
smeared-out step-function as cutoff on k (no cutoff on 
the frequencies) and by dropping the self-energy. 
Antiferromagnetism. For t' = and ^ = 0, the band is 
half-filled and the Fermi surface a perfect square. Ev- 
ery vector connecting parallel sides of the Fermi surface 
is a nesting vector, and Vfk = at (7r,0) and (0,7r). 
This strongly enhances particle-hole terms at wavevector 
Q = (tTjTt). a random-phase approximation summation 
of these bubbles results in a divergent static spin sus- 
ceptibility at Q for any [/ > at sufficiently low T, 
indicating the formation of an antiferromagnetic (AF) 
spin-den sity wave (SD W), in accordance with mean-field 
studies (|Fuldel . Il99l[ ). The basic RG resuhs at low T 
are shown for U = 2t in Fig. (TU] The labelling of the 
= 32 patches along the Fermi surface can be read off 
Fig. [TO] a). Fig. [TO] b) shows as a function of the 
patch indices ki and K2, at A^ ^ 0.16t and with K3 = 1 
(i.e. ka near (— tt, 0)). Strongly enhanced repulsive in- 
teractions appear as a vertical line at K2 = 24 (i.e. for 
k2 — ks = Q), almost ki -independent, and horizon- 
tal line at ki = 24 (corresponding to ki — ka = Q) with 
only a weak dependence on K2, roughly half as large as 
the vertical feature. In an extrapolation where the regu- 
lar profiles are narrowed down to delta functions with an 
appropriate prcfactor J, V^{ki, K2, K3) = ■j(2(5k2-k3,Q + 
(^ki-ka.q), corresponding to a mean-field AF spin in- 



teraction Hamiltonian Jj2{ij)'^^^'^^' ^^^Si • S^, with 



Si = ^c^crCi . The effective Hamiltonian consisting of 
the low-scale hopping term and this interaction exhibits 
AF long-range order at suf ficiently low T. An analysis of 
the flow of suscep t ibiliti es! Halboth and Metzneil . l2000al : 
iHonerkamp et al. 1, I2OOII) as described in Sec. III.FI con- 
firms this picture. 

The extrapolation to a mean-field Hamiltonian is a 
drastic oversimplification, in which the spin fluctuations 
are lost, but they are retained in the obtained by 
the RG flow. As the leading instability is clearly ex- 
posed by this analysis, one can also resort to a bosonized 
description t hat t reats the collective infrared physics 
(jBaier et al\ . \2004 l. 

d-wave Cooper pairing. For t' = — 0.3t and /.( = 
— 1.2t; the Fermi surface still contains the saddle points 
(tt, 0) and (0, tt) but is curved away from these points 
fFig. ITOb)). Now Cooper pair scattering dominates, well 
visible in Fig. ITOH ) on the diagonal lines ki -I- k2 = ( 
I'^i — '*2| = N/2 in terms of patch indices). It is attractive 
when the incoming pair ki, — ki is near the same saddle 
point (±7r, 0) as the outgoing pair ka,— ks, and repul- 
sive when incoming and outgoing pairs are at different 
saddle points. This is the symmetry of the formfactor 
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FIG. 10 (Color online) A^-patch functional RG data obtained 
with the momentum-shell functional RG for the repulsive 
Hubbard model on the 2D square lattice. Upper plots: /i = 0, 
t' ^ and initial U = 2t, lower plots: /i = 1.2t, t' = -0.3t, 
U — 3t. To the left: Fermi surfaces for the two cases and the 
= 32 discretization points for the two incoming ki , k2 and 
the 1st outgoing wavevector ks. To the right: the coupling 
function V'^* (ki, K2, K3) with K3 = 1 and k\ and K2 moving 
around the Fermi surface. The colorbars on the right indicate 
the values of the interactions. 



d(k) do^coskx — cosfcj,) for (i^2_y2 -Cooper pairing. In 
an extrapolation as above, y^(ki,k2,k3) gives rise to 
the mean-field Hamiltonian 



d{k)d{k') cl, .^c\, ^^c-w.ick,t ■ 

k,k' 



which has a d-wave singlet-paired ground state. This 
rf-wavc pairing instability was found in a number 
of studies using different functional RG schemes 



(iHalboth and Metzner. 2000aib: 


Honerkamo, 2001 


Honerkamo and Salmhofed. l2001allb 


; HonerkamD et al.. 


2001 


; Tsai and Marstonl. 2001; Zanchi and Schuld. 1998 


200c 


), in a rather large parameter region. This con- 



stitutes convincing evidence that the weakly coupled 
Hubbard model possesses a d-wave superconducting 
ground state. 

Interplay of AF and SC. In Fig.fTOtl). the sign structure 
of the d-wave term goes together, and fits perfectly with, 
enhanced repulsive interactions near ki = 8 and K2 = 24, 
which are the remnants of the SDW feature in Fig. ITOb ) . 
Their larger width is due to the Fermi surface curvature. 
As A is decreased, these SDW features appear first, due 
to approximate nesting at high scales, and then create 
an attractive component in the d2,2_y 2 -pairing channel, 
which then grows as A is lowered further, while the SDW 
is cut off by Fermi surface curvature, as discussed also 
in Appendix IB. 31 When the SDW-enhancing terms are 
removed by hand from the right hand side of the RG 
equation, the d-wave terms are suppressed as well. Thus 
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the d-wave pairing interaction is induced by AF-spin fluc- 
tuations that appear on higher scales. 

At fixed U, t and t' , there is a sizable interval of /i 
for which the Fermi surface remains close to the sad- 
dle points. Since both AF-SDW and d-wave SC are 
driven by repulsive scattering between (7r,0) and (0, tt), 
both grow and reinforce one another. In the saddle 
point regime, it becomes impossible to single out one 
over the other in the truncation used here. By anal- 
ogy with the quasi-one-dimensional ladder systems, it 
has been argue d that in this regime, the Fermi surface 
gets truncated (iFurukawa et all Il998l : iHonerkamp et all 
I2OOII : iLauchli et all \20oK 



2. Ferromagnetism vs. Superconductivity 

At the Van Hove filling, ferromagnetic (FM) tenden- 
cies are enhanced by the logarithmic divergence of the 
density of states, and the Stoner criterion for the bare 
interaction suggests an FM ordered state at arbitrarily 
small U. However, the Van Hove singularities also make 
the 0{U'^) Cooper pair scattering log^-divergent, hence 
put the two terms into direct competition. 

As discussed in Section IH.D.ll the momentum- 
shell cutoff artificially suppresses FM. For this rea- 
son, the T-flow (see Section | II.p.2| was invented 



( Honerkamp and Salmhofeil . l2001al lb(). and we discuss re- 
sults obtained by T-flow here. The main difference to the 
AF/SC scenario discussed above is that at zero transfer 
momentum, scattering processes driving FM must have 
the opposite sign from those driving singlet SC, hence 
mutually suppress one another. This simple picture is 
confirmed by the RG with momentum-dependent ver- 
tices, in a study where t' and /i are varied at fixed U and 
such that the Fermi surface always contains the saddle 
points: near to t' = — gets strongly suppressed, 
hinting at a quantum critical point between the dSC and 
FM phases (lower left plot in Fig. [TT|) . These results were 
later confirmed by a two-p article self-consistent approach 
( Hankvevvch et a^ . l2003[ ) and in the so-called fi-scheme, 
which employs a soft infrared regulator on the Matsub- 
ara frequencies ( Husemann and Salmhofeiil2009l ): see the 
lower right plot in Fig. [TTJ In the latter study, the 
A'^-patch scheme was replaced by a parametrization of 
the vertex functions in terms of exchange bosons. The 
much higher value of A* in the transitional regime near 
t' = — 1/3 is believed to be due to a form factor that was 
not fully resolved there. 



3. Charge instabilities 

The effective interaction develops a pronounced mo- 
mentum dependence also in the charge sector. In the for- 
ward scattering channel, this amounts to the formation 
of non-uniform contributions to the Landau interaction. 
If strong enough, the latter can lead to a Pomeranchuk 



instability ( Pomeranchukl Il959l ). that is, a symmetry- 
breaking deformation of the Fermi surface. 

In particular, the antiferromagnetic peak drives 
the combination of couplings V'^^(ki, K2, K3) = 
2y^(Ki, K2, K3) — F^(k2, Ki, K3) at certain Q = — ki. 
Near to Q w and Q (tt, tt), 

Vti^u ^2, ^z) ~ -/<i(ki)/d(k2)Ki(k3 - ki) , (113) 

where /d(k) has the same symmetries as d{h) ~ cosfc^ — 
cos ky , but is more strongly peaked near the saddle 
points. For Q = (tt, tt) the corresponding mean-field 
state is the d-dens ity wave state, wh i ch br eaks time- 
reversal invariance (jChakravartv et al I 120011) and gaps 
the single-particle states, except at nodal points on the 
Brillouin zone diagonal. For forward scattering, Q = 0, 
the mean-field state only breaks the lattice rotational 
symmetry of the electronic dispersion and hence of the 
Fermi surface. This tendency to form a nematic state 
( Fradkin et Ml, I2010I) via a d-wave Pomeranchuk in- 
stability driven by forward scatt ering interactions was 
discov ered using functional RG ( Halboth and Metzneil . 
l2000b[ ). Although the Pomeranchuk instability is 
not leading in t he fl ow for the Hubbard model 
(jHonerkamp et all |2002[) . a nematic state can coexist 
with the superconducting s t ate (|Neumavr and Metzneil . 
I2OO3I : lYamase and Metzneil . 12007^ " and it may get less 
suppressed by fiuctuations since it breaks no contin- 
uous symmetry. The d-wave Pomeranchuk instability 
has been investigated possible source of nematic- 
ity of the electronic state in relation with experiments on 
various correlated elec tr on sys te ms (iHonerkampl 12005 
Metlitski and Sachdev . 2010cl: Okamoto et all 



Yamasd . 120091 : lYamase and Metzneil . l2006l ) 
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4. Flows with self-energy effects 

We briefly summarize functional RG studies where the 
self-energy has been included. If a frequency- independent 
vertex function is directly inserted in the right hand 
side of Eq. ([50]) . then is real and independent of the 
frequency, hence only changes the dis persion. This was 
taken into account in the appendix of IHonerkamp et all 
( 200l[ ). where t he adaptive scale decomposition method 
later detailed in ISalmhofed (l200lD was used. To keep the 
density fixed, /x is adjusted as a function of A. Since 
the interaction grows in the fiow, it is a nontrivial check 
of the validity of the truncation that the feedback from 
the interaction does not lift the low-kinetic-energy modes 
to high energies, which would drastically shift the Fermi 
su rface and lead to spurious divergences. The first study 
by IHonerkamp et all ( 200 1[) showed that the Fermi sur- 



face tends to become flatter as A decreases, but that it 
indeed shifts very little before the flow is stopped at A* . 
Thus including the real part of the self-energy does not 
lead to any essential changes in the AF/SC scenario de- 
scribed above. However, correlations that only feed on 
the immediate vicinity of the saddle points, like FM, are 
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FIG. 11 (Color online) Leading instabilities as found by A''- 
patch functional RG in the t-t'-Hubbard model. Left upper 
plot: T, vs. jj, for band filling larger than unity, at t' — —0.3t 
and U = 3t. There is a high-energy-scale AF SDW instability 
with a weaker d.^2_„2-wave pairing instability when the AF- 
SDW is cut off. Data from iHonerkamd (120011 ). Right upper 
plot: Data for the same t' and U on the 'hole-doped' side wit h 
band fillings smaller than one, from lHonerkamp et al\ (|200ll ). 
Now there is a broad crossover 'saddle point regime' between 
the nesting-driven AF-SDW instability and the d^2_y2-wave 
pairing regime. Lower left plot: T* vs. t' at the Van Hove 
filling where the Fermi surface contains the points (vr, 0) and 
(0, tt ). For large t' one finds a ferromagn etic instability. Data 
from lHonerkamp and SalmhofeJ (|2001al lbl) obtained with the 
T-fiow. Right lower plot: $1, vs. t' at Van Hove filling, 
now obtained with the simplifie d vertex paranietrization of 
iHusemann and Salmhoferl (|2009l ) and with a soft frequency 
regulator il. 



affected more strongly, and a full analysis of the cou- 
pled flow of self-energy and vertex directly at the saddle 
points remains an open problem, in spite of partial results 
(fFcldman and Salmhofer„ ,2008lj) . 

The imaginary part and the frequency-dependence of 
the self-energy can be approximated by inserting the in- 
tegrated f low of the in t eracti on vertex in the self-energy 
equation ( Honerkamd . l200ll ). This effectively includes 
two-loop frequency-dependence effects, and captures the 
T^-dependence of the quasiparticle scattering rate in a 
Fermi-liquid situation and the exponent of the vanishing 
quasiparticle weight in the L uttinger liquid up to second 
order in the bare couplings ( Honerkamp and Salmhoferl . 
I2OO3I ). For the 2D Hubbard model, the quasiparti- 
cl e lifetime and ren o rmalization factor was calculate d 
in IHonerkamd ( 200lh : lHonerkamp and Salmhoferl ( 20031) . 
exhibiting a strongly k-dependent quasiparticle degrada- 
tion_as_/Vj is ap proached. This trend was also found 
by IZanchil (120011 ) in a slightly different approximation 
for the self-energ y, and is also r obust in a more elab- 
orate treatment ' (|Kataninl . l2009t) . where the six-point 
vertex was included partially. The anisotropy of the 
quasiparticle lifetime was found to have a non-Fermi- 



liquid temperature dependence and to correlate with 
the s trength of the gener ated d-wave pairing interac- 
tion ( Ossadnik et a/.l . l2008[ ). similar to what is observed 
experimentally in overdoped cuprates. More refined 
studies of the frequency-dependence revealed, however, 
that a simple parametrizati on in terms of a qu a sipar- 
ticle weight is insufficien t ( Katanin and Kampj . 12004 
iRohe and Metzned . [20051) . It was shown that near A*, 
the small-|a;|-behavior of S^(a;,k) leads to a split-up of 
the quasiparticle peak. All these findings are consistent 
with an anisotropic break-up of the Fermi surface that 
one would like to co nnect with the phenomenology of 
thehigh- Tc cuprates ( Honerkamp et aLl . 120011 iLee et all . 
I2006D . but a quantitative comparison is difficult due to 
the strongly coupled nature of the cuprates. 



C. Pnictide superconductors 

The functional RG has been very useful in the study 
of t he newly discovered iron pnictide superconduc- 
tors (iHirschfeld and Scalapinol l2010t llshida et al. 1. 120091: 



n 



20081 ). Here the functional RG may work 



even better, as the pnictides are less strongly corre- 
lated than the high-Tc cuprates. This can already be in- 
ferred from the experimental phase diagram, where one 
only finds metallic antiferromagnetic phases (if at all), 
but never Mott insulating antiferromagnetism. Theoret- 
ical works that try to assess the iron c?-orbital onsite- 
interaction strengths find values that put the materi- 
als into the range of weak to moderate correlations 
(|Anisimov et all l2009t [Mivake et al\ . |2010[ ). Regarding 
the electronic structure, the pnictides are more complex 
than the cuprates. At least three of the five iron d- 
orbitals have non-negl i gible weight near th e Fermi level 
(jPaghofer et all l2010l : iMazin all |2008[ ). Therefore, 
even if one is only interested in the vicinity of the Fermi 
surface, the multi-band character has to be kept. The 
Fermi surface (see Fig. [T^j b)) is divided into two hole 
pockets, centered around the origin of the Brillouin zone 
at k = 0, and two electron pockets around k = (tt, 0) 
and k = (0, tt) in the unfolded zone corresponding to 
the small unit cell with one iron atom (or k = (tt, tt) in 
the folded zone corresponding to the large unit cell with 
two i ron atoms). As po inted out early ( Kuroki et al\ . 
120081 : [M azin et al. 1, Hooi), there is approximate nesting 
of electron- and hole pockets which enhances particle- 
hole susceptibilities with the wavevector connecting these 
pockets. In addi tion, depending o n the parameters and 
approximations ( Ikeda et al\ . \2Q\Bt ). there can be a third 
hole pocket at (7r,7r) in the unfolded zone. 

The first iV-patch studies o f the pnictides were per- 
formed by IWang et al\ (|2009al [2010., 2009bl ) for a five- 
band model. These authors obtained a sign-changing 
s-wave pairing instability driven by AF fluctuations as 
the dominant pairing instability. Further they found 
strongly anisotropic gaps around the electron pockets, 
with possibility of node formation. The basic struc- 
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ture of the phase diagram with the sign-changing pairing 
gap between electron- and hole-pockets can be under- 
stood already fro i n sim plified few-patch RG approaches 
( Chubukov et Ml, I2008D . This would however predict 
isotropic gaps around these pockets ( Piatt et~cil\ . l2009t) . 



To understand the gap anisotropy one has to take into ac- 
count the multi-orbital nature of the electronic spectrum 
in the ir on pnictides , as wa s done already in the initial 
studies (|Wang et all . l2009ai l2010l l2009bl ). In order to 
understand this point, let us start with a single-particle 
Hamiltonian in wavevector-Fe-d-orbital space 



H 



(114) 



where the matrices h(k)oo' take into account intra- and 
inter-orbital terms for orbital index o = o' ox o ^ d 
respectively. s is the spin quantum number. The 
energy bands are obtained by a unitary transforma- 
tion from orbital to band operators (index 6), Ck,b,s = 
U6o(k)ck,o,s. The standard choice for the interaction 
between the electrons is to introduce orbital-dependent 
intra- and inter-orbital onsite repulsions, plus Hund's 
rule and pair hopping terms. While these local terms 
lead to k-independent interactions in the orbital basis, 
parametrized by a tensor Voi^o2,o3,o4, after the transfor- 
mation to bands one arrives at a k-depcndent interaction 
function 



14l,62,63,64(ki,k2,k3,k4) = ^ Kl,o2,o3,o4 Mf)l,ol(ki)lt62,o2(k2)ltb3,o3(k3)lt64,o4(k4) • (115) 

ol,o2,o3,o4 



The combination of UboS behind the inte raction tensor 
is sometimes called the 'o rbital make-up' ( Graser et all 
l2009l : iMaier et all |2009[) . These prefactors cause a 
marked k-structure already in the initial interaction 
which is then renormalized during the functional RG 
flow. It turns out that this orbital make-up has an 
essential influence on the competition between differ- 
ent channels in the flow and is responsible for the 
gap anisotr opies found i n the multi-band funct ional RG 
studies by IWang et al\ (|2009al l2010l l2009bD and in 
subsequent functional RG studies ( Piatt etal\ . l2011at 
iThomale et all l2011a[ ) . A typical result for the predicted 
pairing gaps is shown in Fig. [T^la). Note that according 
to the functional RG analysis, the pairing sta t e shou ld 
be strongly doping-depcndent (iThomale et all l2011allbl . 
[200911 . 
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FIG. 12 (Color online) Functional renormalization group re- 
sults for the iron pnictide compound LaFeAsO at moderate 
hole doping, a) Superconducting form factor as the outcome 
of functional RG, plotted versus the position on the hole pock- 
ets at r and M and electron pockets at X, numbered as 
depicted in b). The competing fluctuations manifest them- 
selves in diverging ordering susceptibilities at low RG scales as 
shown in c), including in particular spin density wave (SDW), 
superconductivity (SC), Pome ranchuk (PI) and char ge den- 
sity wave (CDW) instabilities (jThomale et all\20liali ). 



Summarizing this brief section, the iron superconduc- 
tors pose an interesting problem where the functional RG 
has been instrumental in obtaining the main ordering 
tendencies in good agreement with current experiments. 
For future research, one goal should be to make the 
functional RG a useful bridge between ab-initio descrip- 
tions providing the effective model at intermediate energy 
scales and the many-body effects seen in the experiments 
at low scales. In particular it will be interesting to relate 
experimentally observed materials trends in, e.g., the gap 
structure or the energy scales of the different systems, 
to changes in the microscopic Hamiltonian taken from 
ab-initio descriptions. Furthermore, the functional RG 
studies may have to be extended to include the disper- 
sion orthogonal to the iron-pnictide planes, as this would 
yield additional possibilities for nodes in the gap func- 
tion ( Hirschfc ld and Scalapino j_ 20101 : llshida et al. 1. 120091: 
iNormani . i2008l : iPlatt et all l2011bf) . 



D. Other systems 

Besides the above-described two larger fields of appli- 
cation, the functional RG truncations described in this 
Section have also been employed in a number of other 
models in strongly correlated electron physics. Here we 
briefly list some of these activities. 

In relation to possible unconventional superconductiv- 
ity in organic crystals and layered cobaltates, Hubbard- 
type models on t he triangular lattice hav e been stud- 
ied (iHonerkamd . 120031: iTsai and Marstonl . l200l[ ). At 
large U, the spin exchange interaction between the 
sites of the triangular lattice is geometrically frus- 
trated, leading to a much weaker appearance of an- 
tiferromag netism and a possibl e non-magnetic insulat- 
ing phase (iMorita et al . 2002: Sahebsara and Senechall 
l2008l : lYoshioka et all l2009( ). At weak coupling and for 
nearest-neighbor hopping, Fermi surface nesting is ab- 
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sent, so that near to or at half band fiUing, only low-scale 
Kohn-Luttinger-like superconducting instabilities occur 
out of an innocuous Fermi liquid. However, there appears 
to be a strong dependence on details of the microscopic 
modelling. 

To study interaction effects in graphene, the A^-patch 
functional RG has been applied to the extended Hubbard 
model on the honeycomb lattice. In nominally undopcd 
graphene, the Fermi surface becomes a set of Dirac points 
where the density of states vanishes, and no instabilities 
arc found for sufficiently small interactions. If the in- 
teraction strength exceeds a certain value, various insta- 
bilities driven by particle-hole fluctu ations between the 
two Dirac points ( Honerkamd . |2008( ) are found. Inter- 
estingly, for larger second-nearest neighbor interactions, 
there is the possibility of an instabi l ity to wards a quan- 
tum spin Hall phase ( Raghu et all |2008() . However, a 
spin-liquid phase for intermediate strength of the Hub- 
bard onsitc repulsion that was recently found in quan- 
tum Monte Carlo calculations ( Meng et a/.l . |2010| ) is not 
reflected in the functional RG results on this level of ap- 
proximation. When the Fermi level is moved away from 
the Dirac points, the functional RG again detects pair- 
ing instabilities. In the case of dominant nearest neighbor 
repulsion, the l eading pairing tend ency is in the /-wave 
triplet channel ( Honerkamd . l2008l ) . 

The unbiasedness of the functional RG, and the ac- 
cess it gives to k- and w-dependence of vertex functions, 
is also of great use in (quasi-)one-dimensional models. 
The half-filled extended Hubbard model in one dimen- 
sion has been studied in the search for bond-order-wave 
phases, which could indeed be found with a refined patch- 
ing of the k-de pendence o f the in teraction away from the 
Fermi points |Tam et all |2006[) . For quasi-lD models 
with a small transverse hopping in a second direction 
the change from a gapless Luttinger liquid in a strictly 
one-dimensional situation to Fermi liquid instabilities to- 
ward ordering can be monitored as a function of th e 
transverse hopping ( Honerkamp and Salmhofeii |2003[) . 
The F ermi surface in coupled metall ic c hains was stud- 
i ed by iLedowski and Kopiet j ( 2007al lbl): iLedowski et al\ 
(|2005l ). The possibility of triplet pairing driven by den- 
sity wa ve fluctuatio n s has been explored in such situ- 
ations ( Nickel et al I I2005L I2006D . In these quasi-one- 
dimensional systems, including the frequency dependence 
of the interaction vertex is numerically more feasible 
than in two dimensions. This has been used to study 
the interplay of phonon-media ted and direct ele ctron- 
electr on intera ctions for chains ( Tam et all , l2007a^ , lad- 
ders ( Tam et a l . 20073) and systems with sm all trans- 
verse hopping ( Bakrim and Bourbonnaij . l2010l ). 

Many-fermion lattice Hamiltonians can also be real- 
ized with ultracold atoms in optical lattices, opening up 
new directions. For example, mixtures of mor e than 
two hyperflne states ( Honerkamp and Hofstetteil . |2004| ) 
and boson-mediated pairing on two-dimens i onal lattices 
()Klironomos and Tsail 120071: iMathev et all . l2Q06l l2007l ) 
have been investigated using fermionic iV-patch methods. 



Another promising development is the application 
of the functional RG to quantum spin systems 
(iReuther and Wolfld . [2010h . Here, an auxiliary-fermion 
representation is used for the spins in generalized Heisen- 
berg models, and the functional RG can be formulated in 
terms of these fermions. As important difference to sys- 
tems of itinerant electrons, in the quantum spin system 
the kinetic energy for the pseudo-fermions is zero and 
the interactions only depend on one spatial or wavevec- 
tor variable. This allows one to keep the full frequency 
dependence of the self-energy and interaction vertex on 
the imaginary axis, in the usual truncation where the 
six-point vertex is neglected. The Katanin modification 
(iKatanid (|20Q4| ). see also Section II.C.2) of the flow hier- 
archy turns out to be crucial here. If it is employed, the 
auxiliary-fermion functional RG describes the transitions 
from Neel order to collinear order through an intermedi- 
ate paramagnetic phase in the J1-J2 spin-1/2 model on 
the square lattice as function of J1/J2 in good agree- 
ment with numerical approaches . Furthermore, similar 



syste ms on the triangular lattice (jReuther and Thomale 



2011 ) and with longer-ranged couplings (jReuther et al. 



2011f ) were studied. The success of a relatively simple 



truncation in such an intrinsically strongly coupled sys- 
tem is explained by these authors in that the diagrams 
summed in this flow contain the leading contributions in 
both 1/iV- and l/S'-expansions plus particle-particle dia- 
grams, hence those contributions that are believed to be 
most important. 



IV. SPONTANEOUS SYMMETRY BREAKING 

In many interacting Fermi systems a symmetry of the 
bare action is spontaneously broken at sufficiently low 
temperatures and, in particular, in the ground state. In 
the fermionic flow equations, the common types of spon- 
taneous symmetry breaking such as magnetic order or 
superconductivity are associated with a divergence of the 
effective two-particle interaction at a finite scale Ac > 0, 
in a specific momentum channel. In Sec. HI we discussed 
several examples for such divergences. The truncation for 
the effective two-particle vertex leading to the A^-patch 
scheme described and used in Sec. Ill is insufficient to 
describe the symmetry-broken phase. To continue the 
flow below the scale Ac, an appropriate order parameter 
has to be introduced. 

There are two distinct ways of implementing sponta- 
neous symmetry breaking in the functional RG. In one 
approach the fermionic flow is computed in presence of 
a small (ideally inflnitesimal) symmetry breaking term 
added to the bare action, which is p romoted to a fl- 
nite o rder parameter below the scale Ac (|Salmhofer et aZl 
I2OO4I ). A relatively simple truncation of the exact flow 
equation captures spontaneous symmetry breaking in 
mcan-fleld models such as the reduced BCS model ex- 
actly, although the effective two-particle interactions di- 
verge. Another possibility is to decouple the inter ac- 
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tion by a bosonic order parameter field, via a Hubbard- 
Stratonovich transformation, and to study the cou- 
pled flow of the f ermionic and order parameter fields 
(|Baier et all \2004^. 

In case of competing instabilities a reliable calculation 
based on either of the above-mentioned routes to symme- 
try breaking is quite involved. For a rough estimate of 
order parameters and phase diagrams, one may also ne- 
glect low energy fluctuations and combine flow equations 
at high scales with a mean-field treatment at low scales. 
In this functional RG -I- mean-field approach, one stops 
the flow of the effective two-particle interaction at a scale 
Amf > Ac, that is, before it diverges. The remaining 
low energy degrees of freedom are treated in mean-field 
approximation, with a reduced effective interaction ex- 
tracted from the effective two-particle vertex r'^^^'^"'' . In 
a first application of this "poor man's" approach to sym- 
metry breaking the interplay and possible coexistence of 
antiferromagnetism and d-wave superconductivity in the 
(repulsive) two-dim ensional Hubbard model were studied 
( Reiss et al As in any hybrid method, the results 

depend quantitatively on the choice of the intermediate 
scale Amf (except for mean-field models), and there is 
no unique criterion for this choice. 

We now review the purely fermionic and the Hubbard- 
Stratonovich approaches to spontaneous symmetry 
breaking in the functional RG. The methods will be ex- 
plained for the case of a superconductor as a prototype 
for continuous symmetry breaking, and the reader is re- 
ferred to the literature on applications involving other 
order parameters. 



A. Fermionic flows 

The effective action obtained from the exact 
flow equation or from symmetry-conserving truncations 
thereof exhibits the same symmetries as the bare ac- 
tion S. To analyze spontaneous symmetry breaking, one 
therefore has to add a symmetry-breaking term 5S to 
the bare action and compute the flow of in presence 
of this term. In case of spontaneous symmetry breaking 
an arbitrarily small symmetry breaking term is promoted 
to a flnite order parameter at a scale Ac, which survives 
until the end of the flow. 

A crucial issue is then to find a managable trunca- 
tion of the exact flow equation which captures the essen- 
tial features of the flow into the symmetry broken phase. 
This is non-trivial since the effective two-particle inter- 
actions driving the symmetry breaking become large. In- 
deed, truncations based on neglecting vertices r*^^™^"^ 
with m > 2 in the hierarchy of flow equations fail mis- 
erably. A benchmark for truncations is the requirement 



that they should at least provide a decent solution for 
mean-field models. This r equirement is m et by an ap- 
proximation introduced bv iKataninI ( 2004[ ) to implement 
Ward identities in truncated fiow equations. Katanin's 
truncation, which was described already in Sec. II. C, con- 
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FIG. 13 Coupled flow equations for the self-energy and the 
two-particle vertex determining the fermionic flow with sym- 
metry breaking. 



sists of two coupled flow equations for the self-energy Ti^ 
and the two-particle vertex r'^''^^^, see Fig. [131 They are 
almost identical to the first two equations in the hierar- 
chy described in Sec. II. C, with F^^'''^ = 0, but in the flow 
equation for F*^^)^ the single-scale propagator is re- 
placed by 9a = S^ + G^^aS'^G^. This modification 
takes tadpole contributions obtained from contractions 
of the three-particle vertex F^''^^ into account. 

It is easy to see that the Katanin truncation solves 
mean-field models for symmetry breaking such as the 
Stoner model for ferromagnetism or the reduced BCS 
model exactly ( Salmhofer et al\ . 12004 ) . The exact self- 
energy in such models is given by the Hartree-Fock term 
E = VG (schematically), where V is the bare interac- 
tion, and the two-particle vertex by a ladder sum of the 
form F(4) = V{1 - GGV)-^. These equations hold also 
in presence of a cutoff A. Applying A-derivatives one 
finds immediately that E'^ and F^^^^ obey flow equa- 
tions of the ('schematic) form OaTi'^ = t'^^'^^S'^ and 
9aF(4)A = r(4)A9A(G^G^)F(4)A^ which corresponds ex- 
actly to Katanin's truncation. 

To be more specific, we now consider the case of singlet 
superconductivity, where the continuous U{\) symmetry 
associated with charge conservation is spontaneously bro- 
ken, while spin-rotation invariance remains conserved. 
Superconductivity can be induced by adding a term of 
the form 



5S = 



k 



[Ao(fc)V;t(fc)V"i(-fc) + AS(fc)V';(-fc)V'T(fc)] 



(116) 

with a (generally complex) external pairing field Ao(fc), 
to the bare action. It is convenient to use Nambu spinors 
"^aik) and ^a{k) with ^+{k) = tp^ik), *+(fc) = tp^{k), 
^'_(fc) = -0;(-fc), ^'_(fc) = i>i{-k). The effective action 
as a functional of the Nambu fields, truncated beyond 
two-particle terms, has the form 



k Ctl,Ct2 
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+ i E E ri''^)4^^^^(fci,fc2,fc3,fc4)*a,(fcl)«'c..(fc2)*a3(fc3)*ajfc4) . (117) 
fci ,...,A:4 ai ,...,CK4 



Due to spin-rotation invariancc only terms with an equal 
number of '5 and 4' fields contribute. The Nambu prop- 
agator = (r'^)"^)"^ can be written as a 2 x 2 matrix 
of the form 



G^{k) F^{k) 
F*^ik) -G^{-k) 
(118) 

It is instructive to discuss the flow of the supercon- 
ducting gap and the two-par t icle v ertex for the reduced 
BCS model (ISalmhofer et al\ . \2004) . which is defined by 
an action of the form 



k.a 



{-ika + ^k) ipa{k)%l;„{k) 



+ V{k, k') ^^(fc)^^(-fc)^^(-fc')V'T(fc')-(119) 

k,k' 

Note that the interaction is restricted to particles with 
strictly opposite momenta and spins. It is well-known 
that mean-field theory solves this mod el exactly in th e 
thermodynamical limit ( HaasJ . [l"962; ,Miihlschlegelll962[) . 
The restricted momentum dependence of the bare inter- 
action carries over to similar restrictions for the effective 
two-particle vertex r*-^-*^ in Eq. (|117p . Only two inde- 
pendent components appear, namely 



V^{k,k') = TfJ\_{k,k',k',k) , 
W^{k,k') = Tf_^__{k,k\k\k) . 



(120) 
(121) 



The first component is a normal interaction between two 
particles, and its initial value V^^{k, k') is the bare inter- 
action. The second component is an anomalous term de- 
scribing the creation of four particles. It is initially zero, 
but is generated by charge symmetry breaking terms in 
the course of the flow. Another anomalous term describ- 
ing the destruction of four particles is given by the com- 
plex conjugate of W^{k, k'). The diagonal element of the 
Nambu self-energy vanishes for the reduced BCS model, 
while the off-diagonal element is given by the gap func- 
tion A^(fc). 

For the special case of a momentum-independent s- 
wave interaction V ^ the flow equations obtained from 
the procedure described above are particularly simple. 
Choosing a momentum-independent and real bare gap 
Aq > 0, the flowing quantities A"^, and W'^ are real 
and momentum-independent, too. Their (exact) flow is 
given by 



A^ = -(T/^ + VF^)^ 



I A'^ fixed 



(122) 



propagator, and 



A\2 



-(F^±iy^) 

J2^[\G\krT\F\k)\'].{U3) 



A typical flow for an attractive bare interaction V < 
is shown in Fig. [T31 for two different choices of the bare 
gap Aq. The gap increases monotonically from the ini- 
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where -^F^\/^a Axed is the anomalous Nambu single-scale 



FIG. 14 (Color online) Flow for a reduced BCS model with a 
constant density of states at zero temperature; the band width 
is one and the bare interaction V — —0.3. Left: Flow of the 
gap A'^; the thick line is for a bare gap Ao = 2.4- 10~* and the 
thin line for Aq = 6 • 10~*, in units of the band width. Right: 
Flow of the linear combinations + VF* (solid lines) and 
V* — (dashed lines) of normal and anomalous vertices. 
Thick lines are again for Aq = 2.4 ■ 10^* and thin lines for 
Ao = 6- 10"^ 

tial value Aq upon lowering A, and reaches a finite value 
A ^ Aq for A — 0. A finite Aq regularizes the square- 
root singularity in the gap flow at A = Ac. The normal 
vertex reaches a large negative value at the critical 
scale Ac, while the anomalous vertex becomes large 
and positive. The linear combination + W^, which 
drives the gap flow, is also strongly negative at Ac, but 
it saturates at a moderately negative value for A below 
Ac. By contrast, — decreases monotonically and 
reaches a final value of order I/Aq for A — )- 0, which 
diverges for Aq — > 0. This divergence is the mean-field 
remnant of the Goldstone mode associated with the bro- 
ken continuous symmetry. 

In the case of a discrete broken symmetry, the effective 
interaction becomes large only at the critical scale, while 
no large components remain for A — > 0. This has been ex- 
emplified in a study of the RG fiow of a me an-field mode l 
for a commensurate charge density wave ( Gersch et all 
I2005D . 

The performance of the Katanin truncation for mod- 
els with full (not reduced) interactions has not yet been 
fully explored, since an accurate paranictrization of the 
flowing vertex is quite demanding. However, the re- 
sults obtained so far are encouraging. Staying with su- 
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pcrconductivity as a an example, the (Nambu) vertex 
contains 16 components, most of them corresponding 
to anomalous interactions. In addition to the anoma- 
lous terms appearing already in the reduced BCS model, 
there are anomalous interactions corresponding to the 
creation of three p articles and destruction of one parti- 
cle, a nd vice versa (|Gersch et all 120081: ISalmhofer et all 



Making full use of spin-rotation invariance, all 
the Nambu components can be actually expressed by 
only thre e independent functions of m omenta and fre- 
quencies ( Eberlein and Metzneil . |2010[ ) . The main chal- 
lenge is an adequate parametrization of the (three-fold) 
momentum and frequency dependence, since singulari- 
ties associated with symmetry-breaking and the Gold- 
stone mode appear in the course of the flow. Surpris- 
ingly, in a test case study for the weakly attractive Hub- 
bard model, a rather crude parametrization using the A^- 
patch discretization described in Sec. Ill turned out to 
yield a reasonable flow into the superconducting phase, 
with results for the gap in good a greement with result s 
obtained earlier by other means (jGersch et al I I2008D . 
This is encouraging, but the low energy fluctuations are 
clearly not well described in such a parametrization. To 
deal with the singular momentum and frequency depen- 
dence in the Cooper channel (and possibly also in the 
forward sc attering channel), the channel de composition 
devised by iHusemann and Salmhofeil (I2009D seems very 
useful, since it allows one to isolate singular dependences 
in functions of only one momentum and frequency vari- 
able, similar to a description of singular interactions by 
exchange bosons. The channel decomposition has been 
formally extended already to the superconducting state 
([Eberlein and Metzneil . l2010t ). but a concrete calculation 
beyond mean-field models has not yet been performed. 

In systems with a first order phase transition one may 
miss the symmetry broken phase if one tests only for local 
stability of the symmetric phase by offering a small sym- 
metry breaking field, since the latter may be metastable. 
However, one can escape from the metastable state by 
adding a scale dependent symmetry breaking countert- 
erm to the effective action, which has to be choosen 
sufficiently large at the beginning of the flow and fades 
out for A — 7> 0, such that the system is ultimately 
not modified. Formally this is just another choice of 
regularization within the general framework described 
in Sec. II. B. The counterterm method has been imple- 
mented for the exactly soluble t est case of a charge- 
density wave mean-field model bv lGersch et all ( 2006 



Popular approximations also beyond mean-field the- 
ory can be retrieved from the functional RG by keep- 
ing a suitable subset of contributions. In particular, 
the Eliashberg theory for frequency-dependent (usually 
phonon-induced) pairing interactions can be obtained as 
an approximati on to the e xact flow equations both in 



the svmmetri c (iTsai et al. j_ 20051) and in the sy mmetry- 
broken state ( Honerkamp and Salmhofeil . |2005| ) . This is 



channel for zero transfer in the flow of the interaction, 
and the Fock term for the self-energy. 



B. Flows with Hubbard-Stratonovich fields 

Collective order parameter fluctuations associated 
with spontaneous symmetry breaking in interacting 
many-body systems are often treated by introducing 
an auxiliary order parame ter held via a Hubbard- 
Stratonovich transformation ( Popovl [l987t) . A combina- 
tion of the functional RG with the Hubbard-Stratonovich 
route to spontaneous symmetry bre aking in an i n teract - 
ing Fermi system was flrst used by iBaier et~al. I (I2004D . 
They studied the formation of an antiferromagnetic state 
in the repulsive two-dimensional Hubbard model at half- 
filling and managed to recover the low-energy collective 
behavior (described by a non-linear sigma model) from a 
truncated set of coupled flow equations for the fermions 
and the order parameter field. In the following we de- 
scribe the method for the case of a superfiuid phase, sum- 
marizing the work of several groups. 

We consider an interacting continuum or lattice Fermi 
system with a local attraction V < 0. For continuum 
systems a suitable ultraviolet regularization is necessary. 
A local attraction can act only between particles with op- 
posite spin and leads to singlet pairing. It is thus natural 
to decouple this interaction by a Hubbard-Stratonovich 
transformation with a complex bosonic field (j)(q) cor- 
responding to the bilinear composite of fermionic fields 
VJ2k + This leads to an action of the 

form 



i'a{k) (ikf) - ^k) V'o-(fc) 



rrib 



achieved by keeping the Cooper channel for zero total 
momentum and frequency and the crossed particle-hole 



+ ^[v;t(fc + g)V^;(-fc)0(9)+h.c.] ,(124) 

k,q 

where (f)* is the complex conjugate of and = 
~1/V > 0. 

Spontaneous symmetry breaking can now be stud- 
ied by using the flow equation for the effective action 
T-^[(j),tl},'tp] for coupled bosonic and fermionic flelds de- 
rived in Sec. II. G. Relatively simple truncations capture 
several non-trivial fluctuation effects. Effective interac- 
tions beyond quartic order in the fields are generally ne- 
glected. Also boson-fermion vertices beyond the order 
appearing already in the bare action are discarded. The 
truncations are usually formulated as an ansatz for the 
effective average action 

■0] = r'^[0, ijj, -0] — regulator term 
which obeys the initial condition F^" = S, see Sec. II. G. 
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The ansatz used for the bosonic part is guided by the 
usual strategy of a double expansion in and g radients 
(see, for example. iTetradis and WetterichI ( 19941 )): 



Tb'M = ^ Ui^ciH^)) + gradient terms, 



(126) 



where x = {xo,xi, . . . ,Xd) collects imaginary time and 
real space coordinates. Note that we use the same letter 
(j) for the real space and momentum space representations 
of the bosonic field. The shape of the local potential 
U^^{(t>) depends on the scale. For A above a critical scale 
Ac it has the convex form 



Uti4>)-^mt\<f>\'+u'^\<f>\' 



(127) 



with a minimum at = 0. For A < Ac the potential 
assumes a mexican hat shape 



, A|21 



(128) 



with a circle of minima at \(f>\ = \a^\, where is the 
(flowing) bosonic order parameter. The regime A > Ac is 
called the symmetric regime. At A = Ac the bosonic mass 
TOft vanishes. In the symmetry-broken regime, for A < Ac, 
the order parameter rises continuously from zero to 
a finite value. Its flow can be computed by tracing the 
minimum of the flowing potential or, equivalently, 

by the condition that the bosonic one-point vertex F^,^^'^ 
vanishes. 

For the gradient terms in F^[(/)] various choices have 



been made . Th e simp l est on e ()Birse et all l2005t 

iDiehl et all . l2007l : iKrippal . l2007l ) compatible with the 
U{1) symmetry has the form of an inverse bare prop- 
agator for free bosons, 

J2 [zU*ix)d,.o<l^ix)~At^*ix)y^^ix)] = 

X 

J2 (?) [-^Ztqo + A^q^] ^{q) , (129) 



where V — (dxi dx^)- For lattice fermions one may 
replace by a periodic dispe rsion uig whi c h is p ropor- 
tional to only at small q ( Strack et all |2008| ). The 
term lin ear in qq is a bsent in particle-hole symmetric 
systems ( Strack et all l2008l ) . such that contributions of 
order Qq become important. Additional gradient terms 
have to be taken into account to fully capture the effects 
of the Goldstone mode, as discussed below. 

The normal fcrmionic part of the effective action is 
usually kept in its bare form, sometimes adjusted by 
renormalization factors for the frequency and momen- 
tum dependences. In the symmetry broken regime, an 
anomalous term is generated, such that r4 becomes 



J2 [A^{k)Mk)M-k) + h.c.] . (130) 



For a local interaction the fc-depcndence of the gap func- 
tion A^(/c) is very weak and in simple truncations fully 
absent. Quartic terms corresponding to effective two- 
fermion interactions are absent in the bare action by 
virtue of the Hubbard-Stratonovich decoupling, but are 
generated again in the course of the flow. These gener- 
ated terms are neglected in lowest order truncations, and 
sometimes they are treated by a dynamical decoupling 
procedure called dynamical bosonization, see below. 

For the effective boson-fermion interaction one also 
maintains the bare form of a local 3-point function, 

k.q 

+ anomalous terms , (131) 

with anomalous terms of the form ipipcj) and iptptj)* con- 
tributing only in the symmetry-broken regime. The cou- 
pling is frequently referred to as "Yukawa coupling" . 
The anomalous terms in the boson-fermion interaction 
are usually neglecte d. If taken int o acco unt, they remain 
indeed rather small (jStrack et al I I2OO8D . 

Instead of using a ?7(l)-symmetric ansatz for the ef- 
fective action, one may also start from the hierarchy of 
flow equations for the vertex functions and implement 
the [/ (l)-symmetry by Ward identities ( Bartosch et all 
l2009bD . 

Even with the simple ansatz (|129l) for the bosonic gra- 
dient terms, the effective action described above yields 
sensible results not only at weak coupling, but actually 
in the entire regime from BCS superfluidity to Bo se Ein- 
stein condensation of tightly bound pairs (Die h l et all 
I2OO7I ). In particular, the transition temperature Tc in- 
creases exponentially with the interaction in the weak 
coupling regime, reaches a maximum, and finally satu- 
rates in the strong coupling limit, as it should. 

The bosonic interaction and also the bosonic rcnor- 
malization f a ctors A^ and Zt van ish in the limit A ^ 
(|Birse et all 120051 : iKrippal . |2007[ ). This fluctuation ef- 
fect reflects the drastic renormalization of longitudinal 
order parameter correlations, which are well-known from 
the inter acting Bose ga s in di mensions c? < 3 (sec, for 
example. iPistolesi et al\ ( 20041 )). Note that for a node- 
less gap function the low-energy behavior of a fermionic 
superfluid is equivalent to that of an interacting Bose 
gas, since fermionic excitations are fully gapped. How- 
ever, with the simple ansatz (|129p the transverse order 
parameter fluctuations corresponding to the Goldstone 
mode are also strongly renormalized, which is not cor- 
rect. To distinguish between longitudinal and trans- 
verse fluctuations, one may flx the phase of the or- 
der parameter such that is real, decompose the 
complex order parameter fleld in real and imaginary 
parts (j){q) ~ a{q) + ii^{q) with cr(— g) = o'*('7) and 
7r(— q) = T:*[q), and intro duce different r e norm alization 
factors for a and tt fields (jPistolesi et al I I2OO4D . Using 
this deco mposition, the correc t infrared behavior was ob- 
tained by IStrack et al. I (|2008D where, however, the can- 
cellation of singular contributions to the renormalization 
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factors for the transverse tt fields was implemented by 
hand. To capture this cancellation intrinsically, one has 
to include an additional U{1) symmetric gradien t terrn 
of the form [crO^„,V)g + tt O.^,,, V)7r]^ (|straclj . [2009t 
iTetradis and Wetterichl . Il994[) . 

The fermionic flow based on the Katanin truncation 
described in Sec. IV. A reproduces the exact solution of 
the reduced BCS model (and other mean- field models). 
Within the truncation described above, the bosonized 
flow yields a reasonable solution without artificial fea- 
tures, but the gap comes out a bit too small. The reason 
for this is the truncation of U^^{(f>) at quartic order. To 
re cover the exa. c t solu tion, one has to keep all orders in 
(/) (jStrack et a/.l . l2008( ) . 

The ansatz (|130p for r^[-0, ^p] neglects the generation of 
fermionic interactions by the flow. In particular, quartic 
(two-fermion) interactions are generated by box diagrams 
with four boson-fermion vertices. These terms con- 
tain contributions from particle-hole fluctuations which, 
among other effects, lead to a significant reduction of 
the transition temperature. The (rc-)generated two- 
fermion interaction can be decoupled at each step in 
the fiow by a procedure cahed dynamical bosonization 
(jFlocrchi nger and Wetterichl . 120091 : iGies and Wetterichl . 
|2Q02j 2 0041 ). A general two-fermion interaction cannot 
be decoupled (exactly) by a single Hubbard- Stratonovich 
field, such that several fields may be needed to obtain 
accurate results. Dynamical bosonization was used to 
include effects from particle-hole fluc tuations in attrac- 
tively interacting Fermi systems by iFloerchinger et~al\ 
(|2008l) . 



Following the work of iBaier "el Ml (l2004) on the repul- 
sive Hubbard model at half-filling, functional RG flow 
equations with Hubbard-Stratonovich fields were also ap- 
plied to the Hubbard model away from half-filling. Com- 
mensurate and incommens urate antiferromagne tic fiuc- 
tuations were investigated ( Krahl et all l2009al ). More 
importantly, it was clarified how the generation of d-wave 
pairing from antiferromag netic fluctuations ca n be cap- 
tured by a bosonized flow ( Krahl et a/.l . l2009bh . and the 
fiow was continued into the symmetry broken phase, with 
coupled order parameter fields describi ng antiferromag- 
netism and d-wave superconductivity (jFriederich et all 
I20lll201lh . 

Compared to the purely fermionic RG described in 
Sec. IV. A, the treatment of order parameter fiuctuations 
is facilitated considerably by the Hubbard-Stratonovich 
field. On the other hand, fluctuation effects associ- 
ated with other channels (the particle-hole channel in 
case of superfluidity) look more complicated. For sys- 
tems with competing instabilities the choice of an ade- 
quate Hubbard-Stratonovich field becomes problematic, 
since the fermionic interaction can be decoupled in differ- 
ent ways, which, in combination with truncations, may 
lead to ambiguities in the results. In general, several 
Hubbard-Stratonovich fields must be used, and the anal- 
ysis done in Section IIIII indicates which ones are the 
most important. The decomposition of the interaction 



in ( Husemann and Salmhofeil . l2009l ) allows to switch to 
Hubbard-Stratonovich fields after the fermionic flow has 
been performed down to a certain scale, and may thus 
be used to combine the two flow representations. 



V. QUANTUM CRITICALITY 

Instabilities of the normal metallic state lead to a r ich 
variety of quantum phase transitions ( Sachdevl Il999l ) in 
the ground state of interacting electron systems, which 
can be tuned by a control parameter such as pressure, 
doping, or a magnetic field. Most interesting are contin- 
uous transitions which le ad to quantum critical fluctua- 
tions (jBelitz et al I I2005D . Near a quantum critical point 
(QCP) electronic excitations are strongly scattered by or- 
der paramete r fluctuations such that Ferm i liquid theory 
breaks down (von Loehnevsen et ad l2007t IVoitall2003l) . 
Quantum critical fiuctuations are therefore frequently in- 
voked as a mechanism for non-Fermi liquid behavior ob- 
served in strongly correlated electron compounds. 

Quantum phase transitions in interacting Fermi sys- 
tems are traditionally describe d by a n effe ctive order 
param eter theory pioneered by iHert j ( 1976f ) and iMillisI 
( 19931 ). An order parameter field (j) is introduced by a 
Hubbard-Stratonovich decoupling of the fermionic inter- 
action, and the fermionic fields are subsequently inte- 
grated out. The resulting effective action for the order 
parameter is truncated at quartic order and analyzed by 
standard scaling and RG techniques. However, more re- 
cent studies revealed that the Hertz-Millis approach is 
not al ways applica b le, especially in low-dimensional sys - 
tems ijBelitz et all l2005l : Ivon Loehnevsen et "oil . l2007t) . 



Since electronic excitations in a metal are gapless, in- 
tegrating out the electrons can lead to singular interac- 
tions in the effective order parameter action which can- 
not be approximated by a local quartic term. The na- 
ture of the problem was identified and essential aspects 
of its solution were presented first f o r diso rdered fer- 
romagnets by iKirkpatrick and Belitz and later 

elaborated on by iBelitz et al\ (l2001allbl)~ For clean fer- 
romagnets, IBelitz et al\ (|l997t) showed that Hertz-Millis 
theory breaks down, and no continuous quantum phase 
transition can exist, in any dim ension d < 3 ; the t ransi- 
tion is generically of first order ( Belitz et all\l99^ . The 
Hertz-Millis approach was also shown to be invalid for 
the q uantum antiferromagnet i c transition in two dimen- 
sions (| Abanov and Chubukov . 120041 : 1 Abanov et adl2003l : 
iMetlitski and Sachdevl . l20I0al ). In that case a continuous 
transition survives, but the QCP becomes non-Gaussian. 

Applications of the functional RG to quantum phase 
transitions and quantum criticality in interacting Fermi 
systems have appeared only recently. In Sec. V.A we ex- 
plain how the Hertz-Millis theory fits into the functional 
RG framework and we review some extensions relying on 
an effective order parameter action truncated at quartic 
or hexatic order. An application of a non-pcrturbative 
truncation, where all orders in (j) are (and must be) kept. 
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is discussed in Sec. V.B. Finally, in Sec. V.C we briefly 
discuss the possibility to study coupled flow equations 
for fermions and their critical order parameter fluctua- 
tions in the functional RG framework, and we refer to 
first steps in this direction. 



where 



(134) 



(135) 



A. Hertz-Millis theory 

In his seminal work on quantum ph ase transitions in 
metallic electron systems. iHert j ( 1976h proposed to de- 
couple the electron-electron interaction by introducing 
an order parameter field via a Hubbard- Stratonovich 
transformation. The resulting action is quadratic in the 
fermionic variables ip and ?/;, which can therefore be in- 
tegrated out. One thus obtains an effective action which 
depends only on (j). Truncating at quartic order in </>, and 
discarding irrelevant momentum and frequency depen- 
dences (in the sense of standard power counting), leads 
to the Hertz action. 
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where S^^^ is a field-independent term, and 

UioM = r<j)^ + uq^\ (133) 

We write our equations for the case of a real scalar order 
parameter for simplicity, using (again) the same letter 
for the real and momentum space representations of the 
field. Except for the frequency dependence, the action 
has the form of a 0^ theory for thermal phase transitions. 
The frequency dependent term stems from low-energy 
particle-hole excitations. Here the dynamical exponent 
z is an integer number > 2 depending on the type of 
transition. Tuning the parameter r one can approach 
the phase transition, in particular the quantum phase 
transition at T = 0. 

The action (|132|) has been analyzed by standard scaling 
and RG techniques. Due to the frequency dependence, 
the scaling behavior at the QCP corresponds to a system 
with an effective dime nsional i ty dpf i = d + z, where d is 
the spatial dimension ( Hertzl . [l976h . As a consequence, 
the QCP appears to be Gaussian in three- and e ven in 
two-d imensional systems. An important insig htbylMillisI 
(|1993D was that the 0^ term in the action is nevertheless 
crucial to obtain the correct temperature dependences 
near the QCP. He derived the temperature dependence 
of the correlation length ^ and other quantities by using a 
perturbative RG with a mixed momentum and frequency 
cutoff. From a functional RG perspective, Millis' scaling 
theory can be viewed as a simple truncation of the effec- 
tive average action r^[(/)] evolving from S[(f>], namely 
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and A parametrizes a mixed momentum and frequency 
cutoff. The fiow equations for the parameters in Eq. (|134p 
are obtained by inserting r;^[0] in the exact functional 
fiow equation (|9ip and comparing coefficients. Due to 
the local form of the interaction, the self-energy is 
momentum and frequency independent such that the pa- 
rameters A^ and remain invariant. The flow of u^, 
which is driven by a contribution of order (u^)^, is im- 
portant only in the marginal case d + z = A and near 
the thermal phase transition at Tc > 0. Hence, most of 
Millis' results on the region around the QCP in the phase 
diagram arc based on an analysis of the fiow of and 
the thermodynamic poten t ial V,^ = TP'"^^ (for a review, 
see lvon Loehnevsen et al\ ( 2007[ )1. 

Various extensions of Millis' analysis were derived 
within the functional RG framework. In particular, an 
extension to the symmetry-broken phase was formulated, 
for cases where the broken symmetry is discrete and does 
not g ap out the fermionic excitations ( Jakubczvk et all . 
120081) . One such case is a nematic transition driven by a 
Pomeranchuk instability of interacting electrons on a lat- 
tice, where the discrete point- group symmet r y of t he lat- 
tice is spontaneously broken ( Fradkin et all 2010t ). The 
symmetry-broken regime was described by the ansatz 
(|134p for r;^ with a quartic local potential which has 
a minimum away from zero: 
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The resulting flow equations were used to compute Tc 
and the Ginzburg temperature Tq below Tc as a function 
of the control parameter r. To access the non-Gaussian 
thermal critical regime near Tc, it is crucial to take the 
flow of the quartic coupling into account. The param- 
eters A^ and Z^ are now scale dependent, too. While 
Z^ remains almost invariant, the flow of A^ is important 
near Tc and gives rise to an anomalous scaling dimension. 
A main result of the calculation was that the leading r- 
dependence of Tc is the same as that of the Ginzburg tem- 
peratures below and above Tc (the latter was calculated 
bv lMillisl(|l993h ). b ut a fairly large Gin z burg region opens 
in two dimensions ( Bauer et al. . 1201 lU Jakubczvk et all 
I2008D . 

In another extension a ^^-interaction was included in 
^loc *° study a possi ble change of the o rder of the transi- 
tion by fluctuations ( JakubczvM 20091) . as well a s quan - 
tum tricritical points in metals (jjakubczvk et aZ.1 . l201Clt ). 

Note that the extensions mentioned above are based on 
perturbative truncations resulting in flow equations with 
few running couplings, which could have been obtained 
also by more conventional RG methods. 
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B. Full potential flow 

We now review an application to a problem where the 
effective action cannot be truncated at any finite order 
in (j), such that the possibilit y to make non-perturba tive 
truncations becomes crucial ( Jakubczvk et adl2009l ). 

The problem arises when asking how a nematic transi- 
tion caused by a d-wave Pomcranchuk instability in two 
dimensions is affected by fluctuations. Such a transition 
can be modelled by tight-binding electrons on a square 
lattice wit h an attractive d-wa ve forward scattering in- 
teraction ( Metzner et a/.l . [20031 ) : 
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/kk'(q)rik(q)nk'(-q) , (137) 
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where rik(q) = 
of lattice sites. The interaction has the form 



a ''k-q/2,cr k+q/2,cr 



and L is the number 



ansatz (|134p for r;^[(/)] without expanding the local po- 
tential U^^{(j)). The fiow of C^ioc('/') is then determined 
by a partial differential equation which contains a second 
derivative of the potential with respect to (j). 

In Fig. 15 an exemplary plot of the evolution of the 
flowing effective potential U(^^[(j)) is shown for A rang- 
ing from the ultraviolet cutoff Aq = ~ 0.37 to the 
final value A = 0. The fiow has been computed for 
electrons on a square lattice with nearest neighbor hop- 
ping i = 1, next-to-nearest neighbor hopping t' = — i, 
and a coupling strength g = 0.8. The initial (mean- 
field) potential has a minimum at 4>q = 0.112. The final 
potential exhibits spontaneous symmetry breaking with 
an order parameter 0o = 0.102. Fluctations shift 00 
toward a slightly smaller value compared to the mean- 
field solutions. The fiat shape of U^^{(j)) for (j) < (j)o at 
A = is imposed by the convexity property of the grand 
canonical potential. The transition line between normal 



/kk'(q) = -g(q)dkdk' 



(138) 



where dk = cosfc^; — cos ky is a form factor with (1^2 _y2 
symmetry. The coupling function giq) > has a max- 
imum at q = and is restricted to small momentum 
transfers by a cutoff Aq. For sufficiently large g = g{0) 
the interaction drives a d-wave Pomeranchuk instability 
leading to a nematic state with broken orientation sym- 
metry, which can be described by the order parameter 



(139) 
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In the plane spanned by the chemical potential and tem- 
perature a nematic phase is formed below a dome-shaped 
transition line Tc(/i) with a maximal transition temper- 
ature near Van Hove filling. In mean-field theory, the 
phase transition is usually first order near the edges of 
the transition line, that is, where Tc is r elatively low, and 
second order at the roof of the dome (iKee et al. 1. 120031: 
iKhavkine etaH 12004 lYamase et oil . |2005[ ). 

Introducing an order parameter field via a Hubbard- 
Stratonovich transformation, integrating out the 
fermions, and keeping only the leading momentum and 
frequency dependences for small q and small go/|q| leads 
to a Hertz- type action S[(j)\ of the form (|132p . with z = 3 
and a local potential given by the mean-field potential 



FIG. 15 (Color online) Flowing effective potential Ui^c(.'t>) for 
various values of A between Aq = and 0, at /i = —0.78 
and T = 0.05 (| Jakubczvk et d.Ll2009l ). 



and symmetry-broken phases is shown in Fig. 16 for two 
choices of Aq. Compared to the corresponding mean-field 
result, the transition temperature is suppressed, with a 
larger reduction for larger Aq (corresponding to a larger 
phase space for fiuctuations). For Aq = 1 the transi- 
tion is continuous down to T = 0, leading to quantum 
critical points at the edges of the nematic dome. Increas- 
ing Aq further (or reducing g), one may even eliminate 
the nematic phase c ompletely from the phase diagram 
([Yamase et aLl . l201l[ ). 



At low temperatures, the coefficients of a Landau ex- 
pansion of [/ioc(0) in powers of the field, U{4)) = ao + 
a2<^^ + a4(/)^ -l- . . ., are typically negative for all exponents 
2m > 4. Hence, S[4)\ and consequently also the effec- 
tive action r^[(/)] cannot be truncated at any finite order 
in (j). Fortunately, for bosonic fields the functional RG 
allows also for non-perturbative approximations, where 
on e expands on l y in gradients, and not in powers of 
4> ( Berges et al\ . 20021 ). In particular, one can use the 



C. Coupled flow of fermions and order parameter 
fluctuations 

There are various systems where integrating out the 
electrons leads to an effective order parameter ac- 
tion with singular interactions wh ich cannot b e ap- 
proximated by a loca.1 cou pling ( Belitz et al\ . 120051 : 
Ivon Loehnevsen et al\ . |2007[ ). In such cases it can 
be advantageous to keep the electrons in the action. 
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FIG. 16 (Color online) Critical temperatures versus chemical 
potential for Ao — (larger dome with dots and crosses) 
and Ao = 1 (smaller dome with dots). Th e mean-field tran- 
sition line is also shown for comparison (|Jakubczvk et al\ . 
120091 '). 



treating the coupled system consisting of electrons 
and their order parameter fluctuations. Several cou- 
pled boson-fermion systems exhibiting quantum crit- 
icality have already been analyzed by v a rious meth- 
ods; se e, for example, lAbanov et aZj ( 2003 ): Belitz et al. 
2001ai rbl): Huh and Sachdev ( 200S1 : iKaul and Sachd~ 



(2008 1) : iMetlitski and SachdevI (|2010ai bf) : iRech et al. 



2006r i: IVoita et aZ.I (|2000[) . The interplay of bosonic and 



VI. CORRELATION EFFECTS IN QUANTUM WIRES 
AND QUANTUM DOTS 

As our last application of the functional RG to cor- 
related fermion systems we discuss many-body effects 
in quantum wires and dots. The focus is on transport 
through such systems which arc coupled to two or more 
semi-infinite leads. While in many of the above applica- 
tions it was crucial to devise an approximation scheme 
in which the flow of the two-particle vertex (the effective 
two-particle interaction) was properly described, in the 
ones reviewed in this section the physics is dominated 
by the flow of the self-energy. We start out with a brief 
discussion of quantum transport through a region con- 
taining correlations in Sec. IVI.AI To study transport be- 
yond the linear response regime the functional RG was 
recently extended to Fermi systems out of equilibrium. 
In Sec. IVI.BI we review the main steps of this gener- 
alization. After discussing the most elementary exam- 
ple of linear transport through an inhomogcneous cor- 
related quantum wire - a chain with a single local im- 
purity - in Sees. IVI.C.HlVl.C.41 we show how the func- 
tional RG can be used (i) to describe transport on all 
energy scales for more complex setups (Sec. IVLC.Sp . (ii) 
to identify unconventional low-energy fixed points of such 
systems (Sec. IVI.C.Gp . and (iii) to study finite bias non- 
equilibrium transport (Sec. IVI.C.7[) . As an example of 
the application of the functional RG to quantum dots, 
in Sec. IVI.D.2I we consider an interacting chain of only 
three lattice sites corresponding to a dot. 



fermionic infrared singularities at the quantum critical 
point poses an interesting problem. 



The functional RG for coupled bosons and fermions 
described in Sec. II. G provides a suitable framework 
to study such problems. So far, it has not been ap- 
plied to quantum phase transitions in metallic elec- 
tron systems. However, encouraging works on rela- 
tivistic field-theoretic models with gapless bosons and 
fermions have already appeared. For example, func- 
tional RG flow equ ations have been used to study the 
Gross-Nev eu model (iRosa et all 200 1[) , quantum electro- 
dynamics (jGies and Jaeckel . 20041), and supersymmetric 



Wess-Zumino models ( Gies et al . I2009D . In the context 



of condensed matter physics, a toy model for a semimetal- 
to-superfluid quantum phase transition has been studied 
by coupled flow equatio ns for the electrons and the super- 
fluid order parameter ( Obert et all l201ll : IStrack et all 
I2OIOD . In dimensions d < 3 the fermions and the order 
parameter fluctuations acquire anomalous scaling dimen- 
sions at the QCP of that model, leading to non-Fermi 
liquid behavior and non-Gaussian criticality. 

It will be interesting to devise suitable truncations 
of the coupled boson-fermion flow equations for mag- 
netic and nematic quantum phase transitions in low- 
dimensional metallic systems, where many open ques- 
tions need to be clarified. 



A. Quantum transport 

Experimental progress has made it possible to mea- 
sure transport through mesoscopic regions like one- 
dimensional (ID) quantum wires of up to microme- 
ter length and "zero-dimensional" quantum dots. The 
experiments prov i de evidence for correlati on effects 
(jPeshpande et all I2OIOI : iHanson et ail . 120071) .'^ It is a 
great theoretical challenge to describe transport when 
correlations in the mesoscopic system are important. 
Usually the leads to which this correlated region is con- 
nected are modeled as non-interacting. A general formal 
expression in terms of Keldysh Green functions for the 
current / through an interacting system coupled to two 
le ads (indices L, R) in the sta tionary state was presented 
eir and WingreenI (|1992| ). Either for specific models 
or applying certain approximations to the two-particle 
interaction for each channel C it can be brought into a 
Landauer-Biittiker type form (jBiittikerl . Il986l : iLandauerl . 



^ Correlation effects in effectively ID electronic systems are 
also studied using pliotoemission. For a recent review, see 
iGrioni et aJ\ | |2009|) . 
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Tcie,T,V,)[fL{e)-fRie)]de, (141) 



where we set e = 1 = /i such that the conductance quan- 
tum per channel is given by l/(27r). Here 7^ is an ef- 
fective transmission probability, Vf, = /.jl — l^ift. the bias 
voltage, and f^ /r are Fermi functions with the chemical 
potentials of the left and right leads fJ-L/R- For the trans- 
port through a non-interacting system 7^ is the single- 
particle transmission probability. The goal is to compute 
T( in the presence of correlations. Here truncations of the 
functional RG flow equations which lead to frequency in- 
dependent self-energies are considered. In this case using 
Eq. (|14ip . with being proportional to the "contact-to- 
contact" matrix element of the (retarded) one-particle 
Green function (see below) docs not present an addi- 
tion al approxiination as cur rent vertex corrections van- 
ish I Enssl . [2005t lOguri l200l[ ). In this approximation the 



two-particle interaction affects the transport only via the 
renormalized self-energy which acts as an additional, T 
and Vf, de pendent scattering potential on non-interacting 
electrons (jLanger and Ambegao kar. 1961; Oguri, 2001). 
For the linear conductance (T) the transmission prob- 
ability enters only at zero bias Vb = 0, 



ZTT 



Tc{e,T,0)[-^)de, (142) 



i.e. 7^(e,T, 0) is an equilibrium property. At zero 
temperature Eq. (|142|) simplifies further to 5^(0) = 
rc(M,0,0)/(27r). 



B. Functional RG in non-equilibrium 

Recently, the functional RG approach was extended 
to study steady state non-equilibrium transport through 
quantum wires and dots in the presence of a finite bias 



of the left and right leads Vt = — ^ 


B (iGezzi et al.. 


2007: Jakobsl 2003. 20101: Jakobs et all l2007al 2010a 


Karraschl 20101: Karrasch et al. 2010c 


) . The basic 



idea behind this extension is the use of real time or 
real frequency Green funct ions on the Keldysh contour 
( Rammer and Smithl . I1986D instead of Matsubara Green 
functions. As usual in diagrammatic approaches based 
on Keldysh Green functions one assumes that the ini- 
tial st atistical operator d o es no t contain any correla- 
tions ( Rammer and Smithl . [l986l) . One can then either 



use a functional integral formulation (iKamenevl. 12004 ) 



of the non-equilibrium many-body problem (jGezzi et al 



20071) or a purely dia grammatic approach (jJakobs . .2003 ; 
Jakobs et all l2007a[) to derive the (same) flow equations 
for the self-energy and higher order vertex functions. Al- 
though the method allows to work with two-time Green 
functions and to study transient dynamics, in the current 
implementation of the functional RG in non-equilibrium 



for fermions the system is assumed to be in the steady 
state. For interacting bosons functional RG was also 
used to study dynamics ( Gasenzer and Pawlowskil 120081 : 
iKloss and Kopietzl . [20TTI ) . On a technical level and com- 
pared to the equilibrium Matsubara functional RG in the 
steady state the Keldysh structure only leads to an addi- 
tional index (the so-called Keldysh index ± referring to 
the upper and lower branch of the Keldysh contour) to 
be added to the set of quantum numbers. 



One of the main challenges of the non-equilibrium func- 
tional RG is to devise cutoff schemes which do not violate 
causalit y and general Kubo-Martin-Schwinge r (KMS) re- 
lations (jJakobk l2010t iJakobs et all . l2010bl) after trun- 
cation of the infinite hierarchy of fiow equations. For 
a general cutoff fulfilling the requirements discussed in 
Sect. III.Dl it is only guaranteed that causality and KMS 
relations hold up to the truncation order, e.g. first or- 
der for the level- 1 truncation. In fact, an infrared (real) 
frequency cutoff sim ilar to Eg. (11431) violates causal- 
ity in second order ( Gezzi et all 120071 : IJakobsl l2003l) . 



Its breaking constitutes a severe problem as relations 
connecting the Keldysh contour matrix elements of the 
Green function cannot be used. In particular, one can- 
not rotate from the G°'^ Green function, with Keldysh in- 
dices a, 6 = ±, to retarded, advance d, and Keldysh Green 
functions as it is usually done ( Rammer and Smithl . 
I1986D . The momentum cutoff scheme used in other sec- 
tions of this review avoids this problem but is not suit- 
able for models with broken translational invariance. We 
will return to this issue and discuss appropriate cutoff 
schemes when reviewing applications of non-equilibrium 
functional RG to transport through wires in Sec. IVI.G.7I 
and dots in Sec. IVI.D.21 



Besides its ability to be applied to non-equilibrium 
problems the real time (or real frequency) Keldysh func- 
tional RG has a distinct advantage even in equilibrium: 
the analytic continuation from Matsubara to real fre- 
quencies can be avoided. In fact, computing the real fre- 
quency dependence of observables obtained by a numer- 
ical solution of the Matsubara RG flow equations trun- 
cated on a level which includes a flowing two-particle ver- 
tex and self-energy of initially un known frequency struc- 
ture p resents a formidable task ( Karrasch et all l2008l 
l2010bD as it is known from other imaginary time quantum 
many-body methods such as quantum Monte-Carlo. This 
advantage was utilized in the real time functi onal RG 
study of the single-impur ity Anderson model ( Jakob^ 
120101 : iJakobs et all [2010al) in which a frequency depen- 
dent two-particle vertex and self-energy (complete level-2 
truncation) were kept (for th e imaginary time ana log of 
this study, see Sec. IVLUll iHedden et all (|2004l ). and 
iKarrasch et al. M2008D ). and can be expected to be useful 
also in functional RG studies of other models. 
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C. Impurities in Luttinger liquids 

The metallic state of correlated ferniions in one 
dimension is a non-Fermi liquid. It fa ll s int o 
the Luttinger liquid (LL) class ( Haldan j . Il980[) . 
This state of matter is characterized by a power- 
law decay of space-time correlati on function s with 
interaction depend e nt ex p onents ( Giamarchil 2004 ; 



Luther and Peschell \l97^ iLuttingeil. Il963l: iMattis 



19741: iMattis and Liebl Il965l: ISc honha mmed. 12005) and 
spin-charge separation (jDzyalosh inskii and Larki 
iMeden and Schonhammerl . Il992t IVoitl. Il993l) . In par- 
ticular, after Fourier transforming, the 2A:i;^-component 
of the density-d ensity response function show s a po wer- 
law divergence (|Apel and Ricd . Il982l: iMattisl . Il974[ ) for 
repulsive interactions, instead of a logarithmic one for 
non-interacting electrons (Lindhard function in ID). This 
indicates that any local inhomogeneity with a non- 
vanishing 2fci? component strongly affects the low-energy 
physics and thus the transport characteristics of a LL at 
low temperatures. 

Perturbation theory is insufBcient to describe the in- 
teraction effects as it fails to capture the RG flow of the 
inhomogeneity appearing even to first order in the inter- 
action and higher order diagrams diverge. As will become 
clear below the functional RG approach captures the im- 
purity flow and does not require a simplified modeling of 
the inhomogeneity. 



1. A single local impurity-the local sine-Gordon model 

To understand the effect of a single impurity on the 
low-energy physics of a spinless, infinite LL (absence of 
non-interacting leads), bosoniza tion was used first (for 
a revi e w of this me t hod s ee e.g. Ivon Delft and Schoelleil 
(I1998D . iGiamarchil (l2004l) . and ISchonhammeil tm^ . 
Within this approach the Fourier components of 
fermionic density operators are split into their chiral 
parts which obey Bose comm u tation relations in the low 
energy subspace (jTomonagal . Il950l ). In this limit the 
kinetic energy and the two-particle i nteraction can b e 
written as bilinears in these operators (iTomonagal Il950l) , 
while a single-particle scattering term (impurity) gener- 
ally takes a complicated form in the bosonic degrees of 
freedom. It simplifies if only the pure forward scatter- 
ing V{0) and backward scattering V{2kF) contributions 
are kept. Using bosonization to obtain results beyond 
this approximate modelling of the impurity is rather in- 
volved. 

The forward scattering term is linear in the bosons 
and can easily be treated leading to a phase shift. This 
is irrelevant for the physics described in the following. 
The backward scattering term consists of the cosine of 
a local bosonic field and the resulting Hamiltonian is 
known from field theory as the local sine-Gordon (LSG) 
mod el. One can either us e the exact Bethe ansatz solu- 
tion ( Fendlev et a/.l . [19951 ) of this model or a bosonic RG 



which i s perturbative \D.y{2kp ) (jFurusaki and Nagaosal 
Il993al: [k ane and Fished 119921) to obtain analytical re- 
sults. Alternativ ely numerical methods can be applied to 
the LSG model (|Egger and Grabertl . Il995l : iMoon et all 
I1993D . This led to a complete picture of the RG scaling 
of V{2kF) which has direct consequences for the linear 
conductance. The RG flow connects two fixed points, the 
perfect chain fixed point with conductance g = K/{2-k) 
and the open chain fixed point with g = 0. In ID spinless 
fcrmion systems only a single transport channel exists 
and we thus suppress the channel index C from now on. 
Here K > Q denotes the so-called LL parameter which 
depends on the underlying model of the quantum wire 
and its parameters such as the strength and range of the 
two-particle interaction as well as the band filling. Inde- 
pendently of the model considered K < \ for repulsive 
interactions and K > \ for attractive ones, while the 
non-interacting case corresponds io K = 1. The correc- 
tions to the fixed point conductances are power laws s'^p'" 
with the infrared energy scale s (e.g. the temperature T 
or the energy cutoff A in a RG procedure). The expo- 
nents are independent of the bare impurity strength and 
given by 7p = 2[K — 1) and 7o = 2[\/K — 1), respec- 
tively. The sign of the scaling exponents implies that the 
open chain fixed point is stable for repulsive interactions 
and unstable for attractive ones. The opposite holds for 
the perfect chain fixed point. These insights confirmed 
earlier indications that impurities with a backscattering 
component strongly alter t he low-energy physics of LLs 
with repulsive interactions ( Apel and Ricd . [19821 : iMattid 
I1974I ). The behavior can be summarized by saying that 
even a weak single impurity grows and eventually cuts 
the chain into two parts with open boundary conditions 
at the end points. 

The exponent 70 characterizing (for repulsive two- 
particle interactions) the suppression of g on small scales 
is twice the scaling exponent of the local single-particle 
spectral function of a LL close to an open boundary. This 
can be understood by viewing transport across the impu- 
rity as an en d-to-end tunneling betwe en two semi-infinite 
LLs (see e.g. iKane and Fishe il (ll992l) l. 

Bosonization was not only used for an infinite LL wire 
but also for the experimentally more relevant case in 
which an interacting wire containing a single impurity 
is contacted to two semi- infinite non-interacting leads. 
Contacts generically lead to single-particle backscatter- 
ing and thus have an effect similar to the impurity 
(see below). To disentangle the effect of the con- 
tacts and the impurity one often models the contacts 
such that they do not lead to any backscattering. In 
this case and in the absence of the impurity the con- 
ductance takes the maximal value l/(27r) (instead of 
K/(2tt) for an infinite LL). Using bosonization this can 
either be achie ved within th e so-ca l led local Luttinger 
liquid picture Oanzen et all 20061: Masloy and Stond . 



I1995I : iPonomarenkol 119951 ISafi and Schuld . 11991 " or by 

appropriate boundary conditions for the bo sonic fields 
(|Egger and Grabertl . Il997l : lEgger et all I2OOOI ) . The fixed 



39 



points and scaling exponents turned out to be the same 
as in the LSG model. 

For weak two-particle interactions the problem of 
a single impur ity in a LL was also studied u sing a 
fermionic RG (jMatveev et all Il993l: lYue et all Il994[) . 
In this approach a flow equation for the transmission 
coefficient at kp is derived using poor man's RG. An 
extension of this metho d to interactions o f arbi trary 
strength was presen t ed by l Aristov and W61fl"3 ( 20081) and 
lAristov and Wdlfld (I2009D . 

Remarkably, no intermediate fixed points appear 
within the LSG model and the crossover between the per- 
fect and open ch ain is charac terized by a one-parameter 
scaling function fegger et all. 200 0; Fc ndlev et a/.l . ll995t 
iKane and Fisheil . 119921: lMoon~ al\ . Il9"93h . In one spa- 
tial dimension a general impurity can be described by 
matrix elements Vk^k' (only Vkp^-kp = V{2kp) is kept in 
the LSG model). The RG analysis for such an impurity 
involves the coupling of all matrix elements Vk^k' m the 
flow and one might wonder if this leads to an intermedi- 
ate fixed point absent in the LSG model. To lowest order 
in the impurity strength the flow of Vkp,-kp upon lower- 
ing the cutoff A is driven by Vkp.-kp itself (see the next 
section). For repulsive two-particle interactions Vkp^-kp 
increases and the pcrturbative RG breaks down. Now the 
other couplings, absent in the LSG model, might become 
relevant and eventually cut off the flow of Vkp,-kp to 
large values, that is towards the open chain fixed point. 
That this docs not happen can nicely be revealed within a 
functional RG approach. Before reviewing how this ques- 
tion was approached (see Scc. lVI.C^ we first present the 
most elementary functional RG flow equation to tackle 
inhomogeneous LLs and analytically show that it leads 
to the scaling behavior known from bosonization in the 
limit of weak impurities. 



obtain the fixed point value of the effective interaction to 
leading order one can therefore replace the flowing two- 
particle vertex by the antisymmetrized bare interaction 
corresponding to the level-1 truncation introduced in Sec. 
III. C. 21 An improvement which includes the flow of the 
static two-particle vertex is reviewed in Sec. I VI. G. 41 Af- 
ter presenting the most elementary functional RG flow 
equation for the self-energy we show that it leads to the 
correct scaling properties for a weak impurity. 

As one deals with systems in which translational invari- 
ance is broken it is advantageous to introduce the infrared 
cutoff A in frequency space. To set up the functional RG 
flow equations for the self-energy the propagator Go of 
the non-interacting Hamiltonian Hq containing only the 
kinetic energy is replaced by G^(ia;„) = 0^(a;„)Go(ii^n) 
with a function which is unity for |a;„| 3> A and van- 
ishes for |aj„| ^ A. More specifically 
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> A + ttT 
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was used ( Enss et~al\ . |2005| ). where A starts at oo and 
goes down to 0.^ For T — >■ 0, becomes a sharp 0- 
function and Matsubara frequencies with |a;| < A are 
suppressed. 

In this subsection a general continuum or lattice model 
of spinless fermions is considered with one-particle states 
\a), which in the following will be either local states 
where a; = j is the site index for the lattice model, with 
lattice constant a = 1, or momentum states |fc„) with 
kn = 2im/L. A general two-particle interaction Eq. ([2]) is 
assumed and an impurity term Mmp = p ^a,P i'lt'^p- 
The flow for the self-energy Eg. ([50]) in the level-1 trun- 
cation reads ( Meden et al. . I2OO2D 



2. The functional RG approach to the single impurity problem 

As in the previous applications because of the neces- 
sary truncations the functional RG can (presently) only 
be used for small to intermediate two-particle interac- 
tions. For the application to inhomogeneous LLs it is 
crucial that it is non-perturbative in the single-particle 
inliomogeneity. Here the focus is on a description in 
which the RG flow and the interaction dependent expo- 
nents characterizing the physics close to fixed points are 
kept at least to leading order in the interaction. Follow- 
ing the discussion in the last paragraph of the preceeding 
subsection the feedback of the impurity into the flow of 
the self-energy dominates the physics to be studied. It 
is thus advantageous to use the one-particle irreducible 
functional RG scheme in which the full propagator in- 
cluding the self-energy appears on the right hand side of 
the flow equations. For spinless fermions in a homoge- 
neous wire the electron-electron interaction is renormal- 
i zed only by a finite amount of order interaction squared 
( S61voml . ll979l ). A single impurity does not alter this. To 
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dA 



1 - ^"^G^iiui) 



Ua,^,l3,S , (144) 



(5,7 



where Gq and T,^ are matrices. The initial condition is 
given by = V^,^ and Ua,p.^,s denotes the antisym- 

metrized bare two-particle vertex. 

At temperature T — Q and applying Eq. (|60|) to prod- 
ucts of Q- and 5-functions, Eq. (|144[) simplifies to 



dA 



^a,f} — 



a;— ±A 7,5 



^ A significant speed-up of the numerical solution of the difli'eren- 
tial flow equations can be achieved using the alternative cuto ff 
scheme discussed in the Appendix of lAndereassen et al. I ll2006bf) . 
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where 



For |fci| ^ kp one can linearize the dispersion ^_ 



(146) 



is the full propagator for the cutoff dependent self-energy 
and Qq = (Go)^^. The convergence factor e*"° in Eq. 
p45p is relevant only for determining the flow from A — 
oo down to some arbitrarily large Aq . For Aq much larger 
than the band width this high energy part of the flow can 
be computed analytically leading to the initial condition 

Within the present scheme ^ is frequency indepen- 
dent and can be considered as a flowing effective impu- 
rity potential. To obtain an approximation for the Green 
function of the original cutoff free problem for arbitrary 
impurity parameters one has to determine the self-energy 
^ at A = by numerically solving the set of differen- 
tial equations (|144p or (|145p . To compute the right hand 
side of the flow equations one has to invert the matrix 
Eq. (|146p . i.e., to solve the problem of a single particle 
moving in the effective scattering potential E^ ^ . 

Because of the matrix inversion involved in calculating 
the right hand side of Eq. (|145p the flow equations can 
be solved analytically only in limiting cases, one being 
the situation of a weak impurity. One then works in 
momentum space and considers E^j., for k ^ k' . The 

term linear in E^ presents the leading approximation in 
the expansion of on the right ha nd side of Eq. (|145p 
and one obtains ( Meden et aL . [2003 ) 



ki ,k2 
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*A-6: 



-E^ 



, (147) 



where = — ^^ with the one-particle dispersion 
efc. The antisymmetrized two-body matrix element 
Uk,ki,k'.k2 contains a momentum conserving Kronecker 
delta, where k + ki ^ k' + k2 modulo the reciprocal lat- 
tice vector 27rn for the lattice model with n = 0,±1, 
when the four momenta are in the first Brillouin zone. 
The umklapp processes n = ±1 involve low energy exci- 
tations only for special flllings-half filling for the nearest 
neighbor hopping model discussed later. 

First models are considered for which umklapp pro- 
cesses are absent. To determine the backscattering prop- 
erties of the self-energy one can put k = kp and k' ~ 
—kp — q with |g| <C kp. In order to read off the dominant 
behavior for small A the remaining summation variable 
is shifted ki = —kp + ki 

— E^ =-—S^U 

kF,-kF-q 27J- kF,-kF+ki,-kF-q,kF+q+ki 

ki 



1 



1 



kF+q+ki, — kF+ki _ ^ 



'kF+q+ki "^^ S-fc^+fei 

-K(A^-A)j. (148) 



kF 



-vpki and ^ 



vp{ki +q), where vp denotes the 



ykF+q+ki \---L ■ lyi .-J, 

Fermi velocity. In the thermodynamic limit the two Go 
factors for 5 = and A are proportional to S{ki). If 
only the singular contributions are kept the differential 
equation for E^ reads 



_yl' 



UkF 



27rvF 



A^* 



-kf 



(149) 



where UkF,-kF,kF,-kF = LUkF,-kF,kF,-kF is indepen- 
dent of the system size. For the continuum model 
UkF,-kF,kF,-kF = UiO) - U{2kp), where U{k) is the 
Fourier transform of the two-particle potential U{x — x'). 
This leads to the scaling relation 



^kp 



-kp/C^-^VF) 



(150) 



As Eq. (|149p was derived by expanding the Green func- 
tion in powers of the self-energy, the scaling behavior 
Eq. (|150p can be trusted only as long as E^ stays small. 
Thus a single-particle backscattering term is a relevant 
perturbation for repulsive interactions and an irrelevant 
one for attractive ones consistent with the bosonization 
result. 

Equation (|149p holds even in the half-filled band case 
for the nearest neighbor hopping lattice model. The ad- 
ditional singular contribution due to umklapp scattering 
is proportional to UkF.kF--kF,-kF which vanishes because 



of the antisymmetry of the matrix element (jMeden et al\ . 
120021 ). 

Neglecting all interaction effects beyond the renormal- 
ization of the impurity potential and using the Born ap- 
proximation the correction to the perfect chain conduc- 
tance is given by the self-energy squared. The corre- 
sponding exponent UkF.-kF.k ^.-k^ I (''^vf) agrees to lead- 
ing order in the interaction ( Schonhammeii l2005l ) with 
7p = 2{K — 1) obtained within bosonization. 

The opposite limit of a weak link can be treated an- 
alytically as well leading to results consistent to those 
of the bosonization approach with an exponent charac- 
terizing the deviation from the open chain fixed point 
conductance g — Q which agrees t o leading order in t he 
interaction with 7^ = 2(l//s: - 1) (|Meden et a/.l . l2002t) . 

Next the numerical solution of the RG flow for a spe- 
cific lattice model with arbitrary impurity strength is dis- 
cussed. This allows us to address the question of addi- 
tional fixed points. 



3. Basic wire model 

In the following the tight-binding model of spinless 
fermions with nearest neighbor interaction supplemented 
by an impurity is considered. The Hamiltonian is given 
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FIG. 17 Schematic plot of the quantum dot situ ation, where 
the b arriers are modeled by two site impurities (|Enss et all 
120051 ). 



hy H = i/kin + -ffint + Himp with kinetic energy 



Hkin — — 



E 



4+1 + 4 s+ 



i) 



(151) 



where and Cj are the creation and annihilation op- 
erators on site j, respectively. The corresponding non- 
interacting dispersion is —2 cos A:. The interac- 
tion is restricted to electrons on N neighboring sites 
(|Enss et al\ . 120051) 



i/(n,f/)] [n,+i - K", t/)] , (152) 



Cj Cj . The two re- 



with the local density operator nj 

gions of the lattice with j < 1 and j > N constitute 
the semi-infinite non-interacting leads. To model con- 
tacts which do not lead to single-particle backscattering 
(see above) the interaction Ujj+i between electrons on 
sites j and j -I- 1 is allowed to depend on the position. A 
conductance g = l/(27r) in the absence of single-particle 
impurities is only achieved if f/j.j+i is taken as a smoothly 
increasing function of j starting form zero at the bond 
(1,2) and approaching a constant bulk value U over a 
sufficiently large number of bonds. Equally, the C/j j+i 
are switched off close to the bond {N — 1,N). The re- 
sults are independent of the shape of the envelope func- 
tion as long as it is sufficiently smooth. An abrupt two- 
particle inhomogeneity acts similarly to a single-particle 
impurity. A detailed disussion on the effect of the spatial 
variation of the two-particle int erac tion was presented b y 
iMeden and Schollwockl ( 2003a[) and lJanzen et "aZI (I2006D 



In Eq. (|152p the density nj is shifted by a parameter 
v{n^U)^ which depends on the filling factor n and the 
bulk interaction U . This is equivalent to introducing an 
additional one-particle potential which can compensate 
the Hartree potential in the bulk of the inter acting wire - 
In th e half-filled band case 1^(1/2, [/) = 1/2 (jEnss et all 
I2005D . 

The general form of the impurity part of the Hamilto- 
nian is 



imp 



3,3' 



(153) 



where Vj^ji is a static potential. Site impurities are 



given by Vj 



3,3' 



V,5 



jOjji. For a single site impurity 
Vj = V Sjjg jo is chosen to be far away from both 
lea ds. Impurities close t o the conta cts were di s cusse d 
bv iFurusaki and Nagaosal ( 19961 ) and lEnss et al. 1 (l2005l ). 
Resonant tunneling can be studied considering two site 
impurities of strengths Vi and Vr on the sites ji ~ ji ~ 1 
and jr = jr- + 1- The N]j sites between ji and jr define 
a quantum dot. The effect of a gate voltage restricted to 
the dot region is described by a constant Vg on sites ji to 
jr- This situation is sketched in Fig. 1171 Hopping impu- 



rities are achieved setting Vjj' 



V 



3 ,3 



^3,3 + i ^3', 3 



+ 1- 



For the special case of a single hopping impurity, j+i = 
(t' — 1) Sjjg, the unit hopping amplitude is replaced by t' 
on the bond linking the sites jo and jo -f 1 . In the double- 
barrier problem a hopping ti across the bond {juji) and 
tr across {jnjr) is considered. 

The homogeneous model H = i?kin + -ffint with a con- 
stant interaction U across all bonds~not only the ones 
within [1, A^]-can be solved exactly by the Bethe ansatz 
and K is determine d by a system of coupled integral 
equations (Haldanc, Il980l ). In the half- filled case they 
can be solved analytically leading to 



K = 



2 f U 
— arccos 

TT V 2 



1 -1 



(154) 



for |[/| < 2. At other fillings the integral equations can be 
solved numerically. The model shows LL behavior for all 
particle densities n and any interaction strength except 
at half filling for \U\ > 2 were either phase separation 
sets in (for U < —2) or the system orders into a charge- 
density- wave state (for U > 2). 

Due to the presence of the semi-infinite leads the direct 
calculation of the non-interacting propagator related to 
^kin + -ffimp requires the inversion of a n infinite matr ix. 



Using a standard projection technique ( Tavloii |2000() it 
can be reduced to the inversion of a.n N x N matrix. The 
leads then provide an additional tj„-depende nt diagona l 
one-particle potential on sites 1 and N (Enss et all |2005| ) 



Vl'^fi^ur.) = !^^^|1-W1 



2 \ y {iujn + m)^ j 

xSj,j, (Sij + 6n,j) ■ (155) 

Since the interaction is only non-vanishing on the bonds 
between the sites 1 to the problem including the semi- 
infinite leads is this way reduced to the problem of an 
A^-site chain. In the functional RG it is then advanta- 
geous to replace the projected non-interacting propaga- 
tor Go including t he impurity by a cutoff dependent one 
(lEnss aLl . l2005[) . 



4. Numerical solution of improved flow equations 

With a minor increase in the numerical effort one 
can go beyond Eq. (|144p for the flow of the self- 
energy and include a A-dependent static interaction 
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( Andergassen et all . |2004|) . Its flow equation is derived 
from the general one for the two-particle vertex Eq. ([5^ 
applying the following approximations: (i) The three- 
particle vertex is set to zero, (ii) All frequencies are set 
to zero, (iii) The feedback of the inhomogeneity on the 
flow of the interaction is neglected, (iv) The interaction 
is assumed to remain of nearest-neighbor for m. Then 
obeys a simple differential equation ((Andergassen et 'all . 



(156) 



where the function h{A) depends only on the cutoff A and 
the Fermi momentum kp. The soluti on of the flow equa- 
tion is lengthy for arbitrary fillings ( Andergassen et oil . 
1200411 but has a simple form for half-filling 
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A 



C//(27r) 



(157) 



These approximations are sufficient to correctly describe 
the RG flow of the two-particle v ertex on the Fermi 
surfa ce of the homogeneous system ( Andergassen et oil . 
I2004D to second o rder, as i t is us ually done in the so-called 
g-ology method ( Solvoml . Il979( ). For the inhomogeneous 
LLs the flow of the effective interaction leads to improved 
results for the scaling exponents (see below). 

Within these approximations and using the projection 
of the leads, the self-energy at T > is a frequency- 
independent tridiagonal matrix in real space determined 
by the flow equation j ± 1 G [1, A'^]) 
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Qo(«^n) - eA(c^„)I]A 
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(158) 

where the matrix Qq = (Gq)^^ is the inverse of the pro- 
jected non-interacting propagator with impurity. The 
symbol |w„| sa A stands for taking the positive as well as 
negative frequency with absolute value closest to A.^ The 
initial conditions for E at A = Aq — > oo are independent 
of the precise realization of the inhomogeneity and read 



v^Ao 

^3,3 



[l/2-K",t/)](C/,-i,,+f/w+i) andE 



0. 



The frequency dependence of the self-energy which ap- 
pears in the exact solution in order is not captured 



* To achie ve this result the cutoff sc heme discussed in the Ap- 
pendix of [Andergassen et al. 1 booebl) was used. 
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FIG. 18 Self-energy near a site impurity of strength V = 
1.5 filling n = 1/4, and interaction (7=1; the impurity is 
located at .70 = 51 3 with A'' = 1025 sites. (Data taken from 
[Andergassen et al\ l|2004D .') 



by this scheme. Thus only the leading order is completely 
kept in the flow of E. 

The coupled flow equations can be solved numeri- 
cally by an algorithm wh ich approximately scales as N 
( Andergassen et oZI , |2004| ) . Typically systems of 10^ lat- 
tice sites were considered, roughly corresponding to the 
length of quantum wires accessible to transport experi- 
ments. For the interacting wire of flnite length the energy 
scale (^Tv = vp/N forms a cutoff for any RG flow. The 
flowing self-energy Eq. ()158p depends on the three scales 
T, Sn, and A. Saturation of E'^ for A < T or A < 6n 
sets in "automatically" in contrast to more intuitive RG 
schemes in which the flowing couplings depend on A only 
and the flow is stoppe d "by hand" by replacirig A — » T or 
A Sn, respectively ( Kane and Fished . Il992l : mie et all 
I1994D . 




10° 10^ 10^ 10^ 1000 2000 3000 
j-io 

FIG. 19 Decay of the oscillatory part of the off-diagonal 
matrix element of the self-energy away from a single hopping 
impurity at bond jo, jo + 1. Results for t' — 0.1, jo — 5000, 
— 10*, U = I, n = 1/2, and different temperatures T = 
10"^ (solid line), T = 10"^ (dotted line), T = 10"^ (dashed 
line), and T — 10~* (dashed-dotted line) are presented. The 
left panel shows the data on a log- log scale, the right panel 
on a linear-log scale. For comparison t he left panel conta ins 
a power- law (j — jo)~^ (thin solid line) (|Enss et am2005D . 

Figure [15] shows the self-energy E at the end of the 
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RG flow, for T = in the vicinity of a site impurity 
of intermediate strength. Both the onsite energy j 
as well as the hopping Sjj+i become oscillatory func- 
tions with wave number 2k p and a decaying amplitude. 
The asymptotic value of S^^+i away from the impu- 
rity leads to a broadening of the band due to the in- 
teraction. A more detailed analysis of the oscillatory 



part lAE.j+il 



Eoff , with the spatial av- 



erage Soff, for different T > is presented in Fig. [TH 
The left panel shows that for |j — jo I ~ 10 it decays as 
1/1 j — jo I up to a thermal length scale ^ 1/T (provided 
T > (5jv) beyond which the decay becomes exponential; 
see the right panel. For U > this is the generic behavior 
for large bare impurities or on asymptotical large length 
scales. It is the scattering off such a long-ranged oscilla- 
tory potential-so-called Wigner-von Neumann potential 
( Reed and Simonl . Il975() -which leads to the power-law 
suppression of the conductance and the local spectral 
weight. One can analytically show that the amplitude of 
the asymptotic 1/1?'— ?o I decay determines the exponent 
(jBarnabe-Theriault et aU |2005[) . By virtue of the RG 
flow this amplitude-and thus the exponent-becomes in- 
dependent of the impurity strength. For this reason first 
order perturbation theory fails. It also leads to an oscil- 
latory self-energy which decays as l/|j — jo| but with an 
amplitude wh i ch dep ends on the bare impurity strength 
( Meden et all . |2002|) incorrectly leading to a power-law 



with an impurity dependent exponent. The idea of an 
oscillatory decaying potential is simil arly inherent to a 
poor man's fermionic RG approach ( Yue et all . Il994[) . 
Often these oscillations of the effective renormalized po- 
tential are referred to as Friedel oscillations. This is mis- 
leading as this term is reserved to the spatial oscillations 
of the electron density. In particular, in an inhomoge- 
neous LL the effective potential decays as \j — jo\~^, 
while the d ensity oscillations asympt otically decay as 
\j - jo\~^ (|Egger and Grabertl . Il995[ ). The latter can 
also be sho wn within the functional R G formalism pre- 
sented here ( Andergassen et aZI . |2004| ). The application 
of the self- consistent Hartree-Fock approximation leads 
to an oscillatory self-energy with a const ant amplitude 
and t hus to a charge density wave state ( Meden et all . 
I2OO2D . This is an unphysical artifact of the approxima- 
tion. 



Using scattering theory ( Enss et~ail . |2005[ ) one can 
show that the effective transmission T(e, T) is given by 
the {1,N) matrix element of the single-particle Green 
function T{ek,T) ^ 4sin2 k \{N\G{ek + iO)|l)|2. Via the 
T-depcndcnt self-energy (see Fig. [TO]) G and thus T car- 
ries a temperature dependence. A typical example for 
the T-dependence of the linear conductance g{T) for a 
strong local impurity is shown as the solid line in Fig. 
[5T] It clearly follows the expected power-law behavior 
for (Jat < T < B with the band width 5 = 4. The T^^ 
scaling at larger T is a band effect. For —0.5 < U < 1.5 
and fillings n = 1/2 as well as 1/4 the exponent extracted 
(see lower panel of Fig. [^T|) agrees well with the one of 
the LSG model 70 — 2{1/K — 1), with K taken from the 



Bethe ansatz. Even for U = 1.5 the rela t ive er ror is less 
than 5 percent (see Fig. 5 of lEnss et ^ (|2005D ). Higher 
order corrections in U present in the numerical solution 
of the flow equations (|157p and (|158|) clearly improve the 
result over the one of the perturbative (in the impurity 
strength) analytical solution of Scc. lVI.C?2] which yields a 
purely linear exponent. We emphasize that this improve- 
ment is not systematic as second and higher order terms 
are only partly kept in the RG. A similar agreement can 
be found for 7^ (|Enss et aZ.I . |2005[) . 




x=[T/TJU,n,V)]' 



FIG. 20 One-parameter scaling plot of the conductance. 
Open symbols represent results obtained for U = 0.5, n = 
1/2, and different T and V, while filled symbols were calcu- 
lated for U — 0.851, n = 1/4. Both pairs of U and n lead 
to the same K = 0.85 (within the present approximation). 
The solid line indicates the non-interacting scaling function 
{l+x^)~^ and the dashed one the LSG m odel power- l aw de cay 
with exponent —2/K. (Data taken from lEnss et all (|2005l ).) 

Within the LSG model no intermediate fixed points 
appear which is reflected by one-parameter scaling 
g = gK{x)/2'K wi th x = \T I Tq\^ ^ and a non- 
universal scale Tn (|Egg er et al. \ '2000": Fcn dlev et all . 



ll995l:lKane and Fisheij .[l992: Moon et a/. , 1993|) For ap- 
propriately chosen To data for different T and V{2kF) 
but flxed K can be collapsed onto the iiT-dependent 
scaling function gK{x). It has the limiting behavior 
gK{x) (X 1 — for 2; — >■ and gxix) cx x""^!^ for a; — >■ 00. 
One can perform a similar scaling with data from the 
numerical solution of the flow equations for the micro- 
scopic lattice model considering different V and T as well 
as two sets of [U . n) leading to the same LL parameter 
( Enss et al\ . |2005[ ) . The perfect coUaps of the data of Fig. 
1201 shows that the improved description of the impurity 
flow beyond the single amplitude approximation inherent 
to the LSG model does not lead to additional fixed points. 
The functional RG scaling function shows a sensible K 
dependence. The exponent of the large x power-law de- 
cay is smaller than the non-interacting one —2 (solid line 
at large x) and very close to the LSG model exponent 
~2/K shown as the dashed line in Fig. [201 This has 
to be contrasted to the K independent (non-interacting) 
scaling function g = (1 -\- x'^)~^ resulting from the poor 
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man's fermionic RG ( Yue et all . Il994l) : solid line in Fig. 



"""lOO ^10000 



The functional RG results for a single local impurity 
in a LL show that the LSG model describes the physics 
(two fixed points, exponents, one-parameter scaling) of 
a broader class of models. The same approach was 
also used to study the persistent current through a LL 
ring with a loca.1 impurity prierced by a mag i ietic fl ux 
(iGendiar et all 120091: iMeden and SchoUwockl . r2003allbh 



Aspects resulting from t he spin degree of fre e dom of elec- 
trons were discussed by Andergassen et al\ ( 2006a[ l and 
lAndergassen et al\ ^0^. 



5. Resonant tunneling 

We next review the results on resonant trans- 
port through a double barrier-defin ing an inte r acting 
quantum do t emb edded in a LL (jEnss et all 120051: 
iMeden et all . |2005|) . The setup is sketched in Fig. [T71 
The linear conductance g is characterized by a hierar- 
chy of energy scales. The functional RG is a unique tool 
to access this problem as it provides reliable results on 
all scales. For a fixed dot size Nd and fixed barriers 
Vi/r (or ti/j.) the dot can be tuned to resonance varying 
Vg. Only for symmetric dots with Vi ~ Vr (or ti ~ tr) 
the peak conductance becomes "perfect" gp = l/(27r). 
For asymmetric barriers a backscattering component of 
the single-particle inhomogeneity remains, leading to a 
reduced gp. Due to the interaction backscattering grows 
during a RG procedure and on asymptotic scales the con- 
ductance vanishes with scaling exponent 7o. The same 
holds away from resonance regardless of the ratio Vi/Vr 
(or ti/tr). Thus the non-interacting resonance of finite 
width cither disappears (asymmetric barriers) or turns 
into a resonance of zero width (symmetric barriers). A 
rich T-dependence is found on resonance and for sym- 
metric barriers on which we now focus. Without loss of 
generality we only consider site impurities as barriers. 

The functional RG procedure can directly be applied 
to the double barrier problem. The dot parameters 
only enter via the non-interacting propagator. Figure 
[211 shows the peak conductance gp{T) for a dot with 
high barriers and two different dot sizes. The relevant 
energy scales B, 5no = vp/No, (see below), and 

Sn are indicated by the arrows. For 5no ~ ^ the 
two barriers behave as independent impurities. Using 
scattering theory one can show that in this case gp is 
obtained by adding the resistances of the two barriers 
(|Enss et all . 120051: Ijakobs all . HqOTH). This explains 



why for Nz) = 100, for which this temperature regime 
is clearly developed, gp{T) agrees to the solid line ob- 
tained by taking g(T)/2 of a single site impurity of equal 
heigth as used for the double barrier. Note that it is 
a non-trivial fact that in this temperature regime the 
individual resistances can be added to give the total 
resistance. In the presence of inelastic processes one 
would of course expect this result (resistors in series) 




FIG. 21 Upper panel: gj,{T) for U = 0.5, iV = 10^ Vi/r = 10, 
n — 1/2, and N = 6 (squares), 100 (diamonds). The arrows 
indicate the relevant energy scales B, JjVj^, and 5jv. The 
solid curve shows g{T)/2 for a single barrier with = 10 and 
U = 0.5, N = lO''. Lower panel: Logarithmic derivative of the 
conductance. Da shed line: 7 o : dash ed-dotted line: 7o/2 — 1. 
(Data taken from lEnss et al\ (|2005l ').') 



but they arc absent in the mcsoscopic sctup-the ones 
resulting from the electron-electron interaction are sup- 
pressed by the approximations. In fact, the case of three 
barriers constitutes an exam ple for which adding resis- 
tances does no longer hold (jjakobs et al I l2007bli . For 
T ^ Sno the width of —df/de is smaller than Sno and 
only the resonance peak around e = of T(e, T) con- 
tributes to the integral in Eq. ( 11421) . T he width w of 
this peak vanishes as T^"/"^ /Nu ([Enss et a l.. 2005) lead- 
ing to gp{T) cx T'^"/^^^, The lower bound of this scal- 
ing regime, first discussed using bos onization ( Furusakil . 
1 19981 : iFurusaki and Nagaosai Il993bl ). is reached when T 



equals w, i.e. at T;^^ cx iV^^^ ^ For T < T^^,27rgp 
approaches 1. For T reaching Sn any power-law scaling 
in T with an interaction dependent exponent is cut off 
by the finite size of the interacting part of the quantum 
wire. In addition of identifying the different tempera- 
ture regimes the functional RG approach allows to (i) 
quantify the size of the crossover regime-typically half an 
order of magnitude-and to (ii) obtain results for "non- 
universal" regimes as e.g. realized for No = 6 and Sn^^ < 
T <^ B. For dots with weak barriers and sufficiently 
large No only the regime with scaling exponent 7o/2 — 1 
is realized and for weak barriers and sma ll Nn none 
of th e above powe r- law r egimes emerges ( Enss et all 
l2005l : IMeden et all . l2005f ). Resona nt transport in 
LLs was also studied b y boso nization ( Furusakil . Il998l : 
iFurusaki and Nagaosai 1993b), poor man's fermionic 



RG (Nazarov and Glazman. |2003|; Polvakov and Gornvil . 



120031 ). and numerically (mgle and Eggerl . l20oi r 

The temperature dependence of the peak conductance 
of resonant tunneling nicely exemplifies that the func- 
tional RG approach provides sensible results on all en- 
ergy scales even for problems with a hierachy of scales. 
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Other examples from this class are situations in which 
the wi re-lead conta. c ts are not modeled as being "per- 
fect" ( Jakobs et and models in which the 
leads and contacts ar e described in a more realistic way 
(IWachter et aLl . l2009l) . 



6. Y-junctions 

The power of the functional RG approach to uncover 
unconventional fixed points and the related interesting 
physics was exemplified by discussing a specific junction 
of three ID wires, a so-called Y-junction. The three LL 
wires (index ly = 1,2,3) each of length N and coupled to a 
non-interacting semi-infinite lead via a "perfect" contact 
are described by the basic model discussed in Sec. IVI.C73] 
The symmetric junction pierced by a magnetic flux (f>, is 
sketched in Fig. [22] and given by 



V \^ 



FIG. 2 2 Sketch of the symmetri c Y-ju nction of three quantum 
wires l|Barnabe-Theriault et al\ . l2005l l . 



Hy = -tY^(^c\,^co.i, + h.c.^ +V^no,i, 

3 

X(e^'^/'4,.co,.+i+h.c.) , (159) 



where the wire indices 4 and 1 are identified. The junc- 
tion is characterized by the thr ee parameters ty^V , and 



t^. Using scattering theory (jBarnabc -Theriault et all 
12005]: iBarnabe-Theriault et alUlQQhi .Enss et all |2005[) . 



the [/ = conductance from wire v to wire v' can be 
written as 



4(ImK)^ je"*"* - k\ 
— 3k + 2 cosoil^ 



(160) 



with a single complex parameter n ~ {—V~t\ G\ i)/|iA|- 
The Green function & is obtained for one of the equiva- 
lent disconnected (iy = 0) wires and i S C denotes its 
diagonal matrix element taken at the first site j = 1. It 
is evaluated at energy e 4- iO with e — ?> 0. Equation (|160p 
holds if V and v' are in cyclic order and is independent of 
the pair considered; g^i^i, follows by replacing </> — > —(f>. If 
(j) does not correspond to an integer multiple of tt and for 
generic junction parameters the conductance from v to v' 
differs from the one with reversed indices indicating the 
breaking of time-reversal symmetry. This constitues the 



most interesting situation and we focus on such fluxes. 
In Fig. [23] the conductance from wire v to v' (cyclic) for 
(f) — QAtt is shown for the upper half of the complex k- 
plane. For restored time-reversal symmetry the largest 
conductance allowed by the unitarity of the scattering 
matrix is 2TTgtj,^i =4/9 (denoted the "perfect junction 
value" in the following); even for optimized parameters 
a reflection of 1/9 is unavoidable. 
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FIG. 23 (Color online) The non-interacting conductance 
2Tig^^^i (cyclic indices) as a function of the complex parame- 
ter Hi which in turn is a function of the junction parameters 
tY, V, and tA- The flux is = QAtt. 

For [/ 7^ the Y-junction can straightforwardly be 
treated within the functional RG base d approximation 
scheme ( Barnabe-Theriault et a/.l . l2005l ) . We here focus 
on T = 0. To compute the conductance from Eq. (jl60p . 

must be replaced by the auxiliary Green function G 
obtained by considering S (at the end of the RG flow for 
the full system) as an effective potential for a single dis- 
connected wire setting = (jBarnabe-Theriault et all . 
I2OO5] ). Via the RG flow of S, G develops a dependence 
on {tY,ti^,V), [/, and 5m = vp/N. The latter energy 
scale is a natural infrared cutoff-in contrast to the flow 
parameter A which is artiflcial and sent to 0. A com- 
prehencive picture of the low-energy physics is obtained 
from the dependence of k on Jat. In Fig.[24]each line is for 
a fixed set of junction parameters and 5^ as a variable. 
The flux is chosen as (j) = QAtt and the arrows indicate 
the direction of decreasing 5iq- As ImK has the opposite 
sign of ImGi4 < it is restricted to positive values. 

Equation (|160p allows for four distinguished conduc- 
tance situations (see Fig. [23]) : (i) on the line Im«; = 0, 
9u,v' = gv',u = for almost all Rck; (ii) it is interrupted 
by three points having flux-dependent positions with the 
conductance 2T:g^^^i = 2TTg^in = 4/9; (iii) for a spcciflc 
flux-dependent k one finds 2'Kg,^ ,jt — 1 and 2TTgi^t^ — 0; 
(iv) = = is also reached for |k| — > 00. These 
are the fixed points of the RG flow as is evident from 
Fig. [M] (i) is an interupted line of decoupled chain fixed 
points with vanishing conductances which is stable for 
U > and unstable in the opposite case. Analyzing 
the dependence of on Sn in its vicinity for different 
U one finds that the scaling exponent is independent of 
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0.8, 



FIG. 24 (Color online) Flow of for f/ = -1, n = 1/2, 
and (j) = 0.47r. Arrows indicate the direction for U < 
0. For U > it is reversed . For details see the text 
(jBarnabe-Theriault ei "oil 120051 '). 



(j) given by 70 as obtained for the single impurity, (ii) 
constitute three perfect junction fixed points (circles in 
Fig. [M)) . For U > each of the three fixed points has a 
basin of attraction given by one of the three parts of the 
curve C{4>) (curved thick line in Fig. 1241 interrupted by 
the square) on which the reflection 1 — — '^'^gu',v 

takes a local minimum. The U dependence of the scal- 
ing exponent when approaching one of the fixed points 
along its corresponding line is shown in Fig. 1251 (circles). 
It is independent of </> and for small \U\ it can be fitted 
by U/{Ztt). These fixed points have not been found by 
any method which is based on bosonization and the exact 
dependence of their scaling exponent on K is presently 
unknown. Because of the factor 1/3 in the leading or- 
der it must be different from the if-dcpendence of any of 
the exponents discussed so far. (iii) The basins of attrac- 
tion are separated by the maximal asymmetry fixed point 
(maximal breaking of time-reversal symmetry; square in 
Fig. [M)) . For = 7r/2 this fix ed point was identified by a 
bosonization based approach ( Chamon et aU[2003t ). and 
it was conjectured that the behavior found holds for all 
fiuxes different from integer multiples of tt. The func- 
tional RG results indeed confirm this-at least for small 
to intermediate |/7|-as one finds this fixed point for all 
such 4> and obtains a flux-independent scaling exponent 
which to leading order agrees to the bosonization result 
7Y = 2(A - 1) with A = 4if/(3 + K^) (see Fig. 
The bosonization exponent shows a non-monotonic de- 
pendence on K and thus C/, which the approximate func- 
tional RG approach does not capture. This implies that 
the maximal asymmetry fixed point is unstable for repul- 
sive interactions, and stable for sufficiently small attrac- 
tive ones but turns unstable again for larger attractive 
interactions, (iv) In the mapping of the complex plane 
onto the Riemann sphere the g = 00 fixed point (north 
pole) is part of the projected line of decoupled chain fixed 
points and shows the same stability properties and scal- 
ing dimension. 
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FIG. 25 Scahng expo nents of the Y-jun ction close to the fixed 
points (jBarnabc-Thcriault et a/] . 120051 ). 



The most interesting physics is associated with the per- 
fect junction fixed points which for U > each have one 
stable direction. If the junction parameters of a non- 
interacting system at fixed </> 7^ rmr, m £ No are chosen 
such that the resulting k lies on C{(j)), but not on one of 
the three special points (ii), gi,,,^' ^ and the conduc- 
tance indicates the breaking of time-reversal symmetry 
as expected. Turning on an interaction [/ > the "fine- 
tuned" system flows to one of the perfect chain fixed 
points with equal "perfect" conductances 2Tig^,^i —4/9 
and 27rg^/.i, = 4/9. Therefore, at small energy scales the 
junction conductance does no longer indicate the explicit 
breaking of time-reversal symmetry. For generic junction 
parameters away from C{(j)) one flnds related behavior. 
Close to the line of decoupled chain flxed points the rela- 
tive difference \gi^,v>—gv>,i^\/ {gv,v'+gv' ,v) scales as a power 
law in (5jv with an exponent given by 7o/2 and thus van- 
ishes if [/ > 0. This implies that g^ jji and g^',u become 
equal faster than they go to zero. In that sense for [/ > 
and up to the unstable maximal asymmetry fixed point, 
on small scales the conductance does not show the break- 
ing of time-reversal symmetry-time reversal symmetry is 
"restored" by the interaction. 

Other types of junctions of an arbitrary num- 
ber of LL wires were s tudie d using functional RG 
( Barnabe-Theriault et ai j_ 20051) as w e ll as by the poor 



man s 



2002 



fermionic RG ( Aristov et all 12010: Lai et al 



2002 ) and bosonization based approaches 



iNavak et a[| . ll999l ). 



Chen et al. 



7. Non-equilibrium transport through a contacted wire 

Non-equilibrium functional RG was used to study a fi- 
nite bias transport geometry with an impurity-frec N site 
interacting wire contacted to two non-interacting semi- 
infinite leads by tunnel barriers modeled by reduced hop- 
ping matrix elements as introduced in subs ection I VI . C . 3 1 
^n.i = {tL — 1) and ti\r,N+i = {tn — 1) (jjakobs et al 



l2007af ). In equilibrium this model features a local single- 
particle spectral function Pj{uj) which close to the chemi- 
cal potential, in the vicinity of the contacts, and for repul- 
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sive i nteractions is suppressed: Pj{io) ~ lo~'° (jEnss et al\ . 
I2005D . The hnear conductance behaves as g{T) ^ T'^° 
which can be understood from viewing transport as an 
end-to-end tunnehng between a LL and a Fermi hquid 
lead and using the sum of two resistances as discussed in 
subsection IVI.C.5I 

A cutoff scheme which conserves causality to any trun- 
cation order (Jakobs, 2010) is given by an imaginary fre- 
quency cutoff. The Fermi function of the two leads which 
can be written as a Matsubara sum 



(161) 



is replaced by 



e(K|-A)e'' 



(162) 



Details of this procedure including a discussion of the 
initial condi tions and its relation to the tem perature 
flow sche me |Honerkamp and Salmhofed . l2001al) are pre- 
sented bv I Jakobs ( 2010f ). Within the lowest-order trun- 
cation and after taking the equilibrium limit this cutoff 
implemented for Keldysh Green functions yields the same 
flow equations as the Matsu bara functional RG with the 
frequency cut off Eq. (|143|) (jjakobd . l2010l: IJakobs et 'all . 
I2007al . l2010bh . 

In the presence of a finite bias voltage the level- 1 trun- 
cation scheme (bare two-particle vertex) with the cutoff 
procedure (|162p was applied. As discussed in Sec. lVI.(T2l 
in equilibrium this is sufficient to obtain scaling expo- 
nents correctly to leading order in U. For weak tunnel- 
ing Fl/a = ^^'l/rPo ^ li with pq the density of states 
of the disconnected, non-interacting leads taken at the 
last lattice site, the flow of the retarded non-equilibrium 
self-energy matrix 5]'°*'^ is given by a weighted sum of 
two equilibrium flows 



^ yvret.A 



E 



^ V.cq.A 
dA 



(163) 



where the terms inside the brackets on the right hand 
side are given by Eq. ()158p with the chemical poten- 
tial set to PL or pji respectively (and U). As 
discussed in subsection IVLC.4I each such term leads to 
an oscillatory slowly decaying self-energy originating at 
the inhomogeneity-the tunnel barriers in the present 
case-and extending into the interacting part of the wire. 
The two chemical potentials Pl/r imply two different 

wave numbers 2k^p^^\ Because of the weighting factor 
^x/i^L + ^r) the amplitudes of the two superimposed 
decaying oscillations are generically different and depend 
on the strength of the bare inhomogeneity. The resulting 
non-equilibrium effect of two different and FA-depcndent 



exponents characterizing the scaling of the spectral func- 
tion close to pl Enid pji goes beyond the naive expecta- 
tion that the bias voltage plays the role o f an infrared 
cutoff scale only (see e.g. ISchoelleil (|2009l )). In Fig. [H] 
the local spectral function near the left contact and for 
a restricted energy range around pl/r is shown. Due 
to the finite temperature (T = 10~^) and the finite size 
of the interacting wire {N = 2 x lO''^) the suppression 
at Pl/r is incomplete (cut off by max{r, (5^}), but the 
difference in the exponent is still apparent. A detailed 
analysis shows that the exponents at Pl/r are given by 
IL/R = '^l/r1o{pl/r)/{^l + Tfl) where the argument 
in the open boundary exponent 7o indicates that it de- 
pends on the band filling and thus the chemical potential. 
After adding a third probe lead these "non-universal" 
exponents can be measured in a t ransport experiment 
(lJakobsl . [2010t IJakobs et al. I. l2007at) . 




FIG. 26 Suppression of tlie local one-particle spectral weight 
as a function of energy near the left contact (at site 5) of 
an interacting wire driven out of equihbrium by a finite bias 
current. The parameters are T = 10~*, A'' = 2*, U = 0.5, 
tL = 0.075, tfl = 0.1 5, and p-l/r = ±0.05. (Data taken from 
IJakobs et all l|2007al ').) 



D. Quantum dots 

A spatially confined system featuring a few energy lev- 
els is called a quantum dot. In a transport geometry the 
dot is coupled to at least two leads. Quantum dots show 
interesting physics if all relevant energy scales (e.g. level- 
lead couplings and T) are smaller than the level spacing 
of the isolated system. Due to the strong confinement the 
two-particle interaction on the dot cannot be neglected 
and leads to phenomena such as Coulomb blockade and 
the Kondo effect. 



1. Spin fluctuations 

In the Kondo regime the physics is dominated 
by spin fluctuations. The virtues and limitations 
of the functional RG approach to describe aspects 
of Kondo physics in and out off equilibrium were 
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extensively studied within the single-impurity Ander- l2008l : IPovonl . [20071 ) . Field theoretical methods were ap- 
plied in the scaling limit in which all energy scales are 
much smaller than the band width B. In the follow- 
ing the focus is on this limit. Functional RG results 
for the equilibrium and non-equilibrium properties be- 
yond the scaling limit including a favorable compari- 
son with recent numerical ti me-dependent density- matrix 
renormali zation group data (iBoulat et al. . 2p08[) ar e pre- 
sented by Karrasch et all ( 2010cf ). Karraschl (j2010l ). and 
iKarrasch et al\ (|2010aD . 



(lAndergassen et all 2006a: iBartosch et all l2009a; 


C. Karrasch and Mcden, 120081: Eichlei et al.. |200£: 


Gezzi et al.. 2007; Goldstein et al. 


. 20091 Hedden et al. 


2004: Isidori et al. 2010: Jakobs 


, 20101: Jakobs et al. 


2010al: Karrasch, 20101 Karrasch et all '20061 2007allb, 


2008 


: Karrasch and Meden. 2009: Kashchevevs et al.. 


2009 


:lMeden and Marauardt. 2006; Schmidt and Wolfle, 


2O10 


Wevrauch and Sibold. i2008l et al. 2005, 


2O10 


). As the number of correlated degrees of freedom 



in quantum dots is small the static truncation used for 
LLs was extended to contain all second order processes 
including a frequency dependent two-particle vertex 
and self-energy, ca pturing the f u ll rea l -space as well 
as spin stru c ture (IHedden et al 
Jakobs et aj\ . 12010. 



20041: IJakobsL [201C : 
Karraschl l2010l: I Karrasch et al. 



20081 l2010d) . In fact, the studies of the single-impurity 



Anderson model constitute one of the rare examples of 
the functional RG approach to correlated Fermi systems 
using a complete level-2 truncation (even supplemented 
by parts of the six-point vertex through the replacement 
discussed in the second part of Sec. III.C.2p . Including 
the frequency dependence clearly improves the results 
beyond bare perturbation theory of the same (that is 
second) order but it is presently not possible to reach the 
strong coupling regime in a controlled way. A discussion 
of the problems yet to be solved was pre s ented by 
IKarrasch et al\ (l20 08V'Karrasch et "aH (l2010(J).lKarraschl 
(|2010[ ). I Jakobs et a l (2010a[), and iJakobsl (|2010[ ). An 
alternative way to include the full frequency dependence 
was recently intro d uced for another impurity model by 
ISchmidt and Ens^ ( 201l[ ). In the following a simple 
quantum dot model dominated by charge fluctuations is 
discussed. 



2. Charge fluctuations in non-equilibrium 

The quantum dot model belongs to the class of spin- 
less models introduced in subsection I VI. C . 3l For a three 
site interacting chain (TV = 3) with Ui^2 = Ul and 
C^2,3 = Ur two hopping impurities are located at bonds 
(1,2) and (2,3): ti^2 = — 1, ^2,3 = tR — 1- Lattice 
site 2 constitutes a single-level dot which can be ad- 
justed in energy by a gate voltage Vg (see Fig. [T7]). An 
electron on this site interacts with lead electrons via a 
nearest-neighbor coupling Ul/r which are otherwise non- 
interacting. Choosing v = 1/2 in Eq. (jl52p . Vg=Q cor- 
responds to the particle-hole symmetric point with dot 
occupation (n2) = 1/2. This model is a lattice real- 
ization of the interacting resonant level model (IRLM). 
The use of a variety of analytical as well as numerical 
methods led to a rather complete understa nding of the 
physic s of this model in equilibrium (see e.g. iBorda et al\ 
( 20071 ) and references therein). In addition, the cur- 
rent under a finite bia s voltage ^il = Vt /2 and hb = 
—Vb/2 was investigated (jBorda et a/.l . 120071 : iBoulat et all 



First order perturbation theory in Ui^/n leads to 
logarithmic terms in the self-energy of the form 
Ul /r In {ih /rI B) which in the scaling limit become large. 
They indicate the appearance of power laws in tL/fi with 
U]^/ji dependent exponents. To uncover them requires 
a treatment which goes beyond perturbation theory. In 
the limit of weak to intermediate two-particle interac- 
tions a Kcldysh functional RG approach to the IRLM in 
the level- 1 truncation leads to a comprehensive picture 
of the physics in and out of equilibrium. In particular, it 
allows to identify the relevant energy scales. 

For the present model instead of Eq. (I162D an- 



other cutoff scheme suitabl e for non-equilibrium ( Jakobj , 
I2OIOI: iJakobs et all l2010aD was implemented and tested 
(|Karraschl I2OIOI : IKarrasch et all l2010cD . In this ap- 
proach each of the three interacting sites is coupled to 
its own auxiliary lead-in addition to the coupling of sites 
1 and 3 to the physical leads. The local density of states 
at the contact points of the auxiliary leads is assumed to 
be energy independent (wide band limit) such that the 
hybridization is energy independent and forms an addi- 
tional onsite "energy" iA on each of the three sites. The 
auxiliary couplings are then considered as the cutoff and 
flow from A = (X), at which regularization is achieved, 
down to A = 0, at which the auxiliary leads are decou- 
pled and the original problem is restored. One can show 
that in the lowest order truncation and in the equilib- 
rium limit the Keldysh contour flow equations become 
equal to the equilibrium ones obtained using the Matsub- 
ara formalism with the (at T = 0) sharp energy cutoff 
Eq. (|143|) . Similarly to the imaginary frequency cutofi^ of 
subsection lVI.C.Tl it conserves causality even after trun- 
cation of the functional RG flow equation hierarchy. In 
addition, in the equilibrium limit this so-called reservoir 
cutoff scheme obeys the KMS relation in any truncation 
order (| Jakobsl [2010l : I Jakobs et al. I. l2010ai rbh. 



In the scaling limit and to lowest order in U only flow 
equations for the hydridizations F^ with initial values 
r™' = iipQt\ appear (A = L/R); the flow of the level 
energies of the sites 1 to 3 is of order J7^. The renor- 
malized hybridizations set the width of the resonance at 
= 0. For A being smaller than the band width the 
flow equations for the rates read (F"^ = L^ + F^) 



A 



- VgY + (A + FA)2 



(164) 
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They have the approximate solutions 



^ A 



An 



max{|AiA -K,|,r/2} 



(165) 



The scale Aq is of the order of the band width. Within 
the static approximation the current takes the form of the 
non-interacting expression with the bare hybridizations 



The bias voltage dependence is clearly more involved 
than in Eq. (jl68p . In particular, simple power- law scaling 
with exponent —2poUii is only recov ered in the extreme 
limits of either Vf, 3^ T„ or c 1 ( Andergassen et all . 
I2OIII ) as the exponent of the second term in the denomi- 
nator 2pQUii is small. Off resonance (e.g. at Vg = 0) and 
for Vb ^ r the current is given by 



r'^' replaced by the renormalized ones 



1 TrT 



Li R 



arctan 



Vb/2 - Vg 



arctan 



Vb/2 + Vg 

r 



\Vb/2-V,\ 



2c 



'2poUl / rp \2poUr I ^2 



\VJ2-V,\ 



c \\VJ2+V,\ 



(166 

It turns out to be useful ( Andergassen et all . 1201 ll 
iKarrasch et al. l. l2010aD to introduce the two scales = 
r™Ao/r„, with Tu = T^ + T^, and the asymmetry pa- 
rameter = T^/T^. The same flow equation can be 
derive d using the so-called rea l -time RG in frequency 
space (lAndergassen et al. I, I2OIII: IKarrasch et al. I, l2010al: 
ISchoelleil . I2009D . Within this approach also the re- 
laxation into the steady state was analyzed in detail 
( Andergassen et oZI . l201lUKarrasch et al. I. l2010al) . 

We first review the results obtained for left-right sym- 
metric model with t^ = tn = t' and Ul = Ub. = U 'as wel l 
as particle-hole symmetry Vq = ( Andergassen et all . 



(170) 

The more involved role of Vb is again apparent. A power 
law is obtained in the above studied left-right symmetric 
case or for very strong left-right asymmtery (c <C 1 or 
c 1) (jAndergassen et al\ . I2O1I . 

This concludes the analysis of the IRLM which shows 
that the functional RG can be a tool to obtain a compre- 
hensive picture of the equilibrium and steady-state non- 
equilibrium physics of a dot model dominated by charge 
fluctuations. The approach allows for an unbiased anal- 
ysis of the non-equilibrium rates and cutoff scales. 



I2OIII : iKarraschl |201GI : IKarrasch et~al 



P65p it follows that in this case the maximum of cither 
l^r^;^! = \Vb\/2 or r itself cuts off the RG flow. The charge 
susceptibility x = ~ d (712) /dVg\y is directly given by 

the renormalized width ~ ^^^^ which at V;, = and 
to leading order in U (Fa — F"" on the right hand side 
of Eq. (|165p ) gives the scaling relation 



2010aD . From Eq. VII. CONCLUSION 



X ^ (F-y 



ay 



2poU + 0{U^) 



(167) 



In the non-interacting case x ^ (|p"")-i ^s expected. 
From Eq. (|166p it follows that the current for = Fa <C 
Vf, ^ B is given by 



F - 



ai = 2poU + 0{U^) 



(168) 



Equations p67p and (11681) were also obtained using 



other approa ches ( Borda et am2007l : iBoulat et a/.l . l200S 
iDovonl . [2007h and suggest that the bias voltage merely 
plays the role of an ad ditional infrared cutoff, beside s 
e.g. F or temperature ( Borda and Zawadowskil l2010t) . 
That this is in general not the case is nicely shown by a 
functional RG treat ment away from particle-hole and/or 
left-r i ght symmetry (Andergassen et al. l.[201ll lKarrasdl 



I2OIGI : IKarrasch et al. 



2010 



The differential conductance g = dl/dVb has a max- 
imum when Vn crosses the chemical potential at Vn = 
±Vb/2 (jAndergassen et aLl . [201ll; IKarrasch et am2010a[ ). 
In the on-resonance case the current for V ^ T reads 



(%) 



2poUr 



2F 



j^.2poUl 1 (t^Y'"'^" 1 + c2 ■ 

r / + c{vbj 

(169) 



A. Summary 

The functional RG has proven to be a valuable source 
of new approximation schemes for interacting fermion 
systems. The heart of the method is an exact flow equa- 
tion, which describes the flow of the effective action as a 
function of a suitable flow parameter. The flow provides 
a smooth evolution from the bare action to the final effec- 
tive action from which all properties of the systems can 
be obtained. Approximations are obtained by truncating 
the effective action. In many cases, rather simple trunca- 
tions turned out to capture rather complex many-body 
phenomena. Compared to the traditional resummations 
of perturbation theory these approximations have the ad- 
vantage that infrared singularities arc treated properly 
due to the built-in RG structure. Approximations de- 
rived in the functional RG framework can be applied di- 
rectly to microscopic models, not only to renormalizable 
effective field theories. Remarkably, the functional RG 
reviewed here as a computational tool is very similar to 
RG approaches used by mathematicians to derive general 
rigorous results for interacting fermion systems. 

We have dedicated a large portion of this review to 
general features of the functional RG method for inter- 
acting Fermi systems (Sec. II). After defining the rel- 
evant generating functionals, we have presented a self- 
contained derivation of the exact functional fiow equa- 
tion and its expansion leading to an exact hierarchy of 
flow equations for vertex functions. We have reviewed 
the different choices of flow parameters used so far, along 
with their advantages and disadvantages. Truncations 
and their justification by power-counting have been dis- 
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cussed in detail for translation invariant bulk systems, 
with links to the closely related mathematical literature. 

In Sees. Ill- VI we have reviewed applications of the 
functional RG to specific systems. Most of the approx- 
imations used in these sections are based on relatively 
simple truncations involving only the flow of the two- 
particle vertex and/or the self-energy. Nevertheless a rich 
variety of phenomena associated with low energy singu- 
larities and instabilities is captured. Instead of summa- 
rizing the content of each section, let us merely highlight 
some distinctive features. Sec. Ill reviews functional RG 
work on the stability analysis of two-dimensional electron 
systems with competing instabilities. The main advan- 
tage of the functional RG based one-loop computation of 
the two-particle vertex, compared to other weak-coupling 
approximations, is that particle-particle and particle-hole 
channels are treated on equal footing, such that there is 
no artificial bias due to a selection or a different treat- 
ment of channels. In the conventional many-body frame- 
work a summation of all parquet diagrams would be 
required to achieve this, but a solution of the parquet 
equations is extremely difficult. Spontaneous symmetry 
breaking, the topic of Sec. IV, can be treated either by 
a purely fermionic flow, or by coupled flow equations for 
the fermions and a Hubbard-Stratonovich field for the 
order parameter. It seems that a comprehensive treat- 
ment of all relevant fluctuation effects related to symme- 
try breaking can be achieved. Applications of the fRG to 
quantum criticality, reviewed in Sec. V, have begun only 
recently. Approximations beyond Hertz-Millis theory can 
be obtained from non-perturbativc truncations of the ef- 
fective order parameter action, or by treating fermions 
and order parameter fluctuations in a coupled flow in- 
stead of integrating the fermionic degrees of freedom at 
once. While the applications reviewed in Sees. III-V ad- 
dress translation invariant bulk systems, the purpose of 
Sec. VI is to show how the functional RG can be fruit- 
fully applied to inhomogencous systems such as quan- 
tum wires and quantum dots - in thermal equilibrium 
and also in a non-equilibrium steady state. A strikingly 
simple truncation of the flow equation hierarchy turned 
out to describe a wealth of non-trivial quantum trans- 
port properties characterized by low-energy power-laws 
and complex crossover phenomena. 



B. Future directions 

The number of functional RG based works on inter- 
acting Fermi systems has increased steadily over the last 
decade, but the possibilities opened by this approach are 
far from being exhausted. There are many opportuni- 
ties and challenges concerning both fundamental devel- 
opments of the method and the extension to a broader 
range of systems. 

On the methodological side there are a number of open 
issues. In systems with an instability of the normal 
metallic state, the flow of the effective interactions is not 



yet fully understood, even on the level of truncations in- 
volving only the two-particle vertex and the self-energy, 
since a faithful parametrization of singular momentum 
and frequency dependences of the vertex is not easy. 

The most outstanding challenge is probably to iden- 
tify accurate and computable truncations of the exact 
flow equation for strongly interacting systems such as 
systems close to a Mott metal-insulator transition. It is 
clear that three-particle and higher order vertices cannot 
be discarded in a strongly interacting system. However, 
they will usually not lead to qualitative changes such as 
new singularities. Hence, there is hope that the contri- 
bution from many-body vertices can be absorbed in the 
structure appearing already on the two-particle level. Af- 
ter all, many strong coupling phenomena, including the 
Mott transition, consist essentially in the formation of 
two-particle bound states. To capture effects related to 
strong local correlations, such as the Mott transition, one 
may also try to treat higher order vertices in a local ap- 
proximation. This would make a link to the dynamical 
mean- field theory (DMFT), where all vertices, includ- 
i ng the self-energy, a re approximated by local functions 
(|Georges et adll996l) . 

For systems with strongly interacting order param- 
eter fluctuations there are already a number of non- 
perturbative approximations for bosonic actions on the 
market. The local potential approximation presented in 
Sec. V is only the simplest one. It can be extended by 
taking non-local cont ributions into accou nt, either in a 
derivative expansion ( Berges et al\ . |2002|) , or by includ- 
ing the full momentum or fre quency dependence u p to a 
certain level in the hierarchy (jBlaizot et al\ . |2005() . Such 
approximations may be very useful for studying incom- 
mensurate density wave instabilities in cases where the 
modulation vector of the density wave can be determined 
only after taking fluctuations into account. 

Recently, the functional RG was extended to a real 
time (or real frequency) Keldysh functional RG which 
can be used for studying correlated Ferrn i systems in 
non-equil ibrium (|Gezzi et all l2007t IJakobsl l2003l I2OIOI: 
iKarraschi [jiTflh. First applications, partly reviewed in 
Sec. VI, indicate that also for these type of problems 
the functional RG constitutes a useful tool of outstand- 
ing flexibility. So far only non-equilibrium steady states 
were studied. To investigate a time evolution is tech- 
nically straightforward but requires a significantly in- 
creased computational effort or additional approxima- 
tions. 

With few exceptions, applications of the functional RG 
to interacting Fermi systems have so far been limited to 
purely fermionic one-band systems. There are many ex- 
tensions of this restricted class of systems, where the fiex- 
ibility of the functional RG can be fruitfully used in the 
future. Multi-band models have been studied already for 
the pnictide superconductors, but there are many more 
and qualitatively different models for transition metal ox- 
ides with orbital degrees of freedom. One may include 
phonons and analyze the electron-phonon interaction ef- 
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fects beyond the EUashberg theory. Allowing for disorder 
one may study the complex interplay of interaction and 
disorder effects. It is not hard to generalize the exact flow 
equations for the extensions listed above. The interesting 
task is then to devise suitable truncations. 

Last but not least, the functional RG is an ideal many- 
body tool to be combined with ab initio band structure 
calculations. A lot of work in the last 15 years has been 
dedicated to the ab initio calcu lation of correlated elec- 
tron materials with the D MFT ( Anisimov and Izvumovl . 
I2OIOI: iKotliar ef all . |2006| ). As in DMFT, an arbitrary 
band structure can be used as input for a functional RG 
calculation. Furthermore, one can easily implement non- 
local potentials and non-local two-particle interactions. 
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Appendix A: Wick-ordered flow equations 

In this appendix we present a derivation of Wick or- 
dered flow equations for fermions, which have been used 
for calculations of instabilities and symmetry-breaking in 
the two-dimensional Hubbard model. 

Wick ordered m-particle functions W^^ ar e generated 
from the Wick-ordered effective interaction ( Salmhofeil . 
I1999D 



2 "^0 V 



[x,x] 



(Al) 



The exponent in the Wick-ordering factor is the func- 
tional Laplacian Aqa = (dy^jG^dj^) with Gg = Go — Gq . 
The Wick-ordered interaction converges to V for A — >■ 0, 
since Gq vanishes in that limit. However, the flow equa- 
tions for and the corresponding m-particle functions 
differ from those for V^. The flow equation fo r the gen- 
erating functional reads ( Salmhoferl . Il999l ) 



dA 



'^0 



(A2) 



where the superscript " diff" indicates that the Laplacian 
takes one derivate on the first, and the other on the sec- 
ond factor W"^ on the right-hand side. This equation is 



obtained as follows. Using the definition of and the 
flow equation for V^, one can write 



dA 



dA' 



A AV^ + -A<i'?V^V^) 

C'O 2 Gr, 



Go 



Using the decomposition Aqa 



A% 



factor 1 



A*^ 



factor2 



A^^A (when acting on a product), this yields 9a W"^ 



"=0 lA'i'rfe^^oV 
^ Go 



'^aV^) and thus Eq. dM 



Expanding in powers of Grassmann fields and compar- 
ing coefficients, one obtains a hierarchy of flow equations 
for the m-particle functions T4^(2m)A^ which is illustrated 
diagrammatically in Fig. [271 The line with the dash 



W 



(2m)A 



r(2n)A 



(2m-2n+2j)A 



d 

dA 




FIG. 27 Diagrammatic representation of the flow equations 
for the effective m-particle interactions W^^"^^^ in the Wick- 
ordered version of the functional RG; the internal line with 
a dash corresponds to OaGq, the others to Gq ; all possible 
pairings leaving m ingoing and m outgoing external legs have 
to be summed. 

is due to contractions generated by AI'a in Eq. (|A2[) . 

Go 

the other lines are generated by the exponential of A^'a ■ 

Gq 

Note that the right-hand side of the Wick-ordered flow 
equations is bilinear in the effective interactions, and no 
tadpole terms appear. Note also that the propagators 
connecting the vertices have support for energies at and 
below the cutoff scale A, such that the integration re- 
gion shrinks as A decreases. One might worry that the 
low-energy propagators lead to infrared divergences even 
for A > 0. This is not the case, as can be seen from 
th e general infrared power-counting analysis presented 
bv lSalmhofeJ (|1999D . 



Appendix B: Details of power counting 

1. Propagator bounds 

Here we show, using properties of the dispersion func- 
tion and the Fermi surface, that 

SA<a + blog^ and j|G'^||<cA-i (Bl) 
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where a, b and c are constants that do not depend on A. 
In absence of Van Hove singularities, 6 = 0. We first con- 
sider the case where the self-energy effects are not taken 
into account (they are discussed in Subsection III.E.4P . 

Then one can simply take x^(^) = X> ( ) where 

X>iv) is a fixed (i.e. A-independent) increasing function 
that vanishes at least linearly at r/ ~ 0, tends to 1 as 



J dl = ^ X] / {2tt)'' contains both frequency and momen- 

turn summations. Taking the absolute values inside the 
sum and estimating the factor by its maximum , 
we obtain 



itr (S^ P^)|| <sa||P^ 



(B6) 



V 



oo and satisfies x>{v) < V ^ ^oi large rj. We can with ,sa = max^^ / dk |5'^„,(fc)|. The second sim- 



then verify (|Bip by scaling, as follows. The full propaga- 
tor is G^ik) = {iko-Ckr^x''{k), so \G^{k)\ < f where 
c = ma.x{ri^^x>{v^) : ?7 > 0} is finite. The single-scale 
propagator is 



S\k) 



2 

A3 



(B2) 



Using that the Matsubara sum is a Riemann sum for the 
convergent integral of over fco and introducing the 
density of states N{E) = J d'^kS{^]^ — E), we get 



SA < / dka 



dEN{E) 



/r2 



E^X'> 



,(B3) 



where the 4 instead of 2 gives a (crude) bound for the 
change from the Matsubara sum to the integral for large 
enough /3. Changing variables to p = (fcg 4- i?^) 2 and a 
polar angle Lp, we obtain 



.SA<^y^ P^dpx'>[^)J^ d^Nipcos^). (B4) 

If the density of states N is bounded, using N{E) < 
Nq and scaling out A implies sa < a, with a ~ 
SttNq J p^x'>{p^)dp < 00. In presence of a Van Hove 
point on the Fermi surface, N stays bounded in dimen- 
sions d > 3, but diverges logarithmically for d = 2. In 
this case, the (^-integral contributes an additional factor 
log A. 

Thus (jBip holds. The hypotheses on x> are satisfied 
in particular for the standard strict cutoff functions that 
vanish identically near rj = 0, and which are identically 
1 for ?7 > 1. For such cutoffs, the single-scale propagator 
is nonvanishing only in a "momentum shell" of thickness 
A around the Fermi surface, and the above bounds can 
also be obtained by estimating the k-space volume of this 
shell (see also Section fB.Sp . 



2. Power counting 

Here we prove (|70p to all orders in the running coupling 
expansion. All terms on the right hand side of the flow 
equation for Fr^™'^ are of the form 

itr(s^p^) iLs.)^lJ diY,siAi)PtAks:,i,-i). 

a, a' 

(B5) 

where — p(,2m+2)A ^-^^ ^^^^ term of ([S7)) and given 
by the other summands in (j67p for the other terms, and 



Pill . . . ||Pn||. It implies bounds for all £p -contributions 



pic inequality that we shall use is that ||Pi 

-A 



in terms of \\G'^\ 
bounded by 
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(2m,)A| 



SO that ||^f(2'"'^| 
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SA 



[l|r( 



(2m+2)A| 



IV^^II IIG^ 



p>2 



-.(2m)A|| 
p(2mp)A| 



(B7) 



The power counting is now determined by sa and HG^H. 
Given ((B7| and (|BT|) . the proof of is an effortless 
induction argument. The inductive sc heme proceeds in 
the standard way of ( Polchinskt Il984l ) , namely upwards 
in r > 1 and at fixed r, downwards in to, starting at 
TO = r, where p(,2'"+2)^ _ q_ ^pj-^g induction start r = 1 is 
trivial. Let r > 2 and assume (|70p to hold for all {r',m') 
with r' < r and for r' = r, m! > m. The right hand 
side of (jB7p contains only terms to which the inductive 
hypothesis ((70)) applies. Inserting it, using (|Bip . and 
collecting powers in the form 1 — p ^ (2 — niq) = 1 — m 



and YZ=i{rq -ruq + l) 



TO, we obtain 



I d p(2m)A| 



< t^^'^^sa^-^+Va^'AI 



where the constant jr"^"'' is a weighted sum of products of 
the 7^f' We now use the initial condition vi^"^''^ — 



(B8) 



to write F 



(2m)A 



, take the norm of 



this equation, and use (jBSp . This gives 



||p(2m)A|| < ~(2m) / ^r~,n+ljr^l-,n 



(B9) 



By definition, /a > /a' if A < A', so < /a for all 
£ in the integration interval. Thus the last integral is 

bounded by 7^^"Va'' /a° d^- 5/""+^^^"". Because si is 
at most logarithmic in i, and to > 3, d£ sg^l^'™- < 
i^A^~™s^ with a constant K that depends on a and 

771. This, together with an appropriate choice of 
completes the induction step. 

For TO = 2, doing the last integral increases the power 
of the logarithm by one. This case is discussed in more 
detail in Subsection IB. 31 

For TO = 1, the self-energy term, the same simple 
bound gives H^E'^H < sa /a, so the integral gives a 
contribution of order f\. When a countcrterm is used to 
keep the Fermi surface fixed, the initial condition for E'^ 
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at A = Aq is given by the counterterm, which needs to 
be adjusted such that at low scales A, E'^(0,k) = 0(A) 
whenever = 0. This leads to the self-consistency rela- 
tion mentioned in section Hi. E. 41 

A sim ilar proof can be given in the Wick ordered 
scheme (ISalmhofeil . ll998bD : it is even simpler because 
the double induction used here is replaced by single in- 
duction on r. 

A crucial point in obtaining (|B8p is that all the depen- 
dence on p and on the mq drops out when the power of A 
is collected. It is this property that makes many-fcrmion 
models with short-range interactions renormalizable in 
the strict quantum- field-theoretical sense. The classifica- 
tion in relevant, marginal and irrelevant terms now also 
becomes apparent because for m > 3, the r^^'"-*^ grow 
as A decreases: suppose we add an additional (2m > 6)- 
point interaction vertex 

y(2m)Ao of Qj.jgj. I ^Yie initial 
interaction at Aq. Its insertion at a lower scale A is a fac- 
tor (;^)'"~^ smaller than that of the effective 2m-point- 
vertex created by the two-particle interaction. Thus the 
influence of wanes at low scales - it is irrelevant. 

A simple adaptation of the above inductive argument in- 
deed shows that the inclusion of such additional terms 
with m > 3 in the interaction at Aq changes only prefac- 
tors in the power counting bounds. For m = 2, this sup- 
pressing factor is absent, so that these terms are marginal 
(the more detailed analysis of Section IB. 31 shows how to 
separate the marginally relevant from the marginally ir- 
relevant terms). Moreover, it is clear that this power 
counting breaks down when F^*^^ develops singularities 
as a function of k and w, because then f\ = oo. Finally, 
for m ~ 1, the scale derivative of the self-energy obtained 
by the above argument is of order A^~™, as in (jB8[) . but 
since m = 1, this integrates to 0(1) instead of 0(A^^™) 
- this term is relevant. To get its size back to 0(A^^™) in 
the momentum shell where |^(k)| ^ A, one needs a can- 
cellation by a counterterm, as described briefly in Section 
III. E. 41 In the Taylor expansion required to do the can- 
cellation, the derivative of the self-energy appears. By 
the above power counting, this is a marginal term. In 
the Luttinger model, it is really marginal and causes the 
anomalous exponents. For curved Fermi surfaces in c? > 2 
dimensions, it is seen to be irrelevant by the arguments 
discussed in Section [B. 3. bl 

There is a hard problem hidden in the recursion of 
the constants 7^^™'. In the recursion described above, 
the number of terms that gets added corresponds to 
the number of Feynman graphs with r vertices, which 
grows factorially in r, so that the bound obtained in 
this way for 7^^™^ and c,^.^"'^ is of order r!. If satu- 
rated, it would lead to a convergence problem. How- 
ever, due to the fcrmionic antisymmetry, sign cancella- 
tions in the sum over Feynman diagrams prevent this fac- 
torial from arising. For proofs, we refer interested read- 



ences therein). In their application to propagators with 
Fermi surfaces, these proofs also provide a rigorous 
bas is for the use o f Fer mi surface patches (first used 
m (|Feldman et all Il992l ) and there called "sectors"). 
Patching the Fermi surface has become an essential tool 
also in applications, see Section lllll 

In a typical lattice model, the kinetic energy per par- 
ticle is bounded, so that the fiow is usually started at 
the highest value (the bandwidth) of the kinetic energy, 
Aq. The '^6 used here also cuts off large frequencies. 
Thus the starting interaction is in this case one where 
the degrees of freedom with frequencies j/col above Aq 
have already been integrated over. This starting action 
can be obtained by converg ent perturbation theory, see 
iPedra and Salmhofed (|2008l ). 



pr 

ers to the literature (seelPisertori and Rivasseaul (I2OOOI): 
iFeldman et a/.l (ll998all2002lll992h:lPedra and Salmhofeil 



(|2008l ): ISalmhofer and Wieczerkowskil (|2000( ) and refer 



3. Improved power counting 

This is a refinement of power counting, valid in a large 
class of bulk ferm i on sy s tems in > 2 (iFeldman et al. . 
19961. Il998bl. Il999l l2000t iFeldman and Trubowitd . ll99C : 
Shankaj . ll99l ~ It is the deeper reason behind the emer- 
gence of Fermi liquid behaviour and of dominant Cooper 
pairing tencencies in weakly coupled standard fermion 
systems, and it provides a precise link between Fermi 
surface geometry and scaling properties of the effective 
771-particle vertices in general. 

We discuss this in the absence of self-energy effects, to 
bring out the main effects as clearly as possible. (The 
self-energy changes the Fermi surface; if the interacting 
Fermi surface is regular, the following analysis remains 
essentially unchanged.) We also assume a strict cutoff 
hmction, i.e. x>{v) = for 77 < (1 — S)^, where < i5 < 
1/2 is fixed, and x>{v) = 1 for V^^- Again, this choice 
is not essential; it just simplifies the discussion. 



a. Effects of curvature on power counting 

The integral lA{k) = J dpod'^p\S^{p)\ \G^{±p + fc)| 
arises from the trace on the right hand side of the RG 
equation when all effective vertices and all but one of 
the propagators have been estimated by their max- 
imal values. It thus determines the maximal possible 
value of a term on the right hand side of the RG equa- 
tion, where the dependence on one external momentum 
is kept. In particular, /a is directly relevant for the one- 
loop contributions to the flowing four-point vertex. The 
power counting done above corresponds to the estimate 
Ihik) < \\G^\\ s\ < A~^s\, so that Jj^I\dX grows loga- 
rithmically in A for small A. 

Since = 0, S^{k) = {iko - Ck)"^ dAX^{k) and 

G^{k) = -I^^M^ = cjAo (fc) _ / " dX S^{k). (BIO) 
iko - ?k J A 

The term G"^" is nonvanishing at large frequencies, 
but not important (HG"^"!] < Aq^, hence a factor 
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A/Aq smaller than \\G^\\ when A gets small. Thus 
/ dpod'^p\S^{p)\\G^°lp + k)\ < saA-\ hence its inte- 
gral over A is bounded by a constant). The derivative 
of the strict cutoff function vanishes unless A{1 — S) < 
K^o — Ck| < A, so S^{p) vanishes unless |po| < A and 
p is in the momentum space shell J-' a = {k : |^k| < A}, 
and there, |S'^(p)| < The po-sum in /a gives at most 
2A, and the p-intcgral gives the d-dimcnsional volume of 
the intersection J- a H (k ± of two momentum space 
shells, where one is shifted by k. It follows that 



/A(fc) < O(Ao-i) + - 



Ao 



^voU (J-An(k± j-a)). 



(Bll) 

This links the scaling behaviour of terms in the RG equa- 
tion to the geometric properties of the Fermi surface. 

Obviously, the volume of the intersection is at most as 
large as the volume of Ja itself: vo\d{J^A H (k ± J^a)) < 
vold^A < const. A. Using this in (jBll[) gives the general 
power counting bound mentioned at the beginning of this 
section, lA{k) < const. A~^. Assuming that = ^k, 
this bound is always saturated for k = 0, and also for 
those k for which the shift by k makes the two shells 
overlap over a significant region of the Fermi surface, that 
is, when k is a nesting vector of the Fermi surface. 

For other values of k, the intersection volume can be 
much smaller than that of J^a- A general definition of 
non-nesting was given, and powe r counting bounds wer e 
derived when it is satisfied, by iFeldman et al. 
an d extended to the case with V an Hove singularities 
in ( Feldman and Salmhoferl . l2008al fbl) . Here we only cite 
the result for the case of a strictly convex and positively 
curved Fermi surface without Van Hove singular ities, dis- 
cussed also in the Appendix of ( S almhof eri . 1 1 9 9 9() . In that 
case, and for A < A < wp min|k|, one can show that the 



volume ratio ^'"^ '''Toirf ^^'^^''^ proportional to 



|k|i;F,mi 



if k^J"; 



(2) 



if kGJ" 



(2) 



(B12) 
(B13) 



Here WF,min is the smallest value of |Ve| on the 

Fermi surface, J-)^ is a 0(A)-ncighbourhood of the set 
{2k : ^k = 0} (note that 2k is taken modulo reciprocal 
lattice vectors), and k denotes the minimal curvature on 
the Fermi surface. This is illustrated for A = A in Fig. [511 
In the first case, the intersection is transversal, which de- 
creases the intersection volume by the factor in (|B12p . 
The second case corresponds to a 2fci?-intersection, where 
the curvature in a region of size ^/X determines the in- 
tersection volume, corresponding to (|B13p . In the third 
case, |k| is so small that the volume of the intersection is 
essentially equal to that of Fa- 

The scale in the flow where the improvements set in 
is determined by the curvature of the Fermi surface, be- 
cause there is really only an improvement if the addi- 
tional factors are smaller than one. In cases where the 
curvature is small on large parts of the Fermi surface, as 




FIG. 28 Intersections of a momentum shell around the Fermi 
surface with its translate, as arising in loop integrals on the 
right hand side of the RG equation. When the Fermi surface 
is curved, the intersection volume decreases strongly unless 
the translating momentum is small. 



in the Hubbard model near to half-filling and at small 
ncxt-to-nearest hopping, one thus has an effective nest- 
ing at those scales and at those k where the quotients in 
(|B12[) and (|B13p arc so big that they give a bound that 
is bigger than the trivial bound 1 for the volume ratio. 
Eqs. (|B12j) and (jBlSp imply that for small |k|. 



Ao 



lA(k)dA < const, log 
Ik^F 



(B14) 



(where the constant depends on the curvature of the 
Fermi surface) and that the function re mains bounded 
for |k | not close to zero (f o r det ails, see iFeldman et al\ 
(I1996D : ISalmhofeil (ll998aL I1999D '). Thus, for convex 
curved Fermi surfaces, the four-point function can di- 
verge only at k = and there, only logarithmically (by 
a similar argument, one can see that it can diverge only 
at ko = 0). The particle-particle correction to the vertex 
function has exactly this behaviour. In the particle-hole 
term, there is an additional sign cancellation that re- 
moves the logarithm. The same argument shows that in 
general, divergences can occur only at nesting vectors of 
the Fermi surface. 



b. Uniform improvement from overlapping loops 

An extension of these geometric estimates to two-loop 
integrals of the type J dp J dq S^jp) S^' ( q)S^" ip±q±k) 
is ver y useful for d > 2: it is shown in ()Feldman et all . 
Il996l ) that in absence of nesting and Van Hove singu- 
larities, such integrals contain a scaling improvement 
independently of k. Such two-loop integrals associ- 
ated to graphs with overlapping loops arise when the 
RG equation gets iterated ; the graph classification of 



(Feldman et 



1 gets Iterated ; 
al\ . I1996LI1999D 



shows that in a precise sense, 
the overwhelming majority of graphs in the Feynman 
graph expansion contains one or even two such subin- 
tegrals, hence becomes subleading at low scales. As is 
explai ned in detail in ( Feldman et all Il996l : I S almhof eil . 
ll998aD . in absence of nesting and Van Hove singularities, 
this justifies the particle-particle-ladder aproximation, it 
singles out the Hartree-Fock type contributions to the 
self-energy by scaling arguments, and it allows to show 
that the derivative of the self-energy is RG-irrelevant in- 
stead of marginal. 
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